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1 Introduction 



Deep inelastic lepton-hadron scattering (DIS) has played a seminal role in the de- 
velopment of our present understanding of the sub-structure of elementary particles. The 
discovery of Bjorken scaling in the late nineteen-sixties provided the critical impetus for 
the idea that elementary particles contain almost pointlike constituents and for the sub- 
sequent invention of the Parton Model. DIS continued to play an essential role in the long 
period of consolidation that followed, in the gradual linking of partons and quarks, in the 
discovery of the existence of missing constituents, later identified as gluons, and in the 
wonderful confluence of all the different parts of the picture into a coherent dynamical 
theory of quarks and gluons - Quantum Chromodynamics (QCD). 

In more recent times the emphasis has shifted to the detailed study of the x- 
dependence of the parton distribution functions and to the study of their Q^-evolution, 
probably the most direct test of the perturbative aspects of QCD. 

Polarized DIS, involving the collision of a longitudinally polarized lepton beam on 
a polarized target (either longitudinally or transversely polarized) provides a different, 
complementary and equally important insight into the structure of the nucleon. Whereas 
ordinary DIS probes simply the number density of partons with a fraction x of the mo- 
mentum of the parent hadron, polarized DIS can partly answer the more sophisticated 
question as to the number density of partons with given x and given spin polarization, in 
a hadron of deflnite polarization, either along or transverse to the motion of the hadron. 

But what is quite extraordinary and unexpected ab initio is the richness and subtlety 
of the dynamical effects associated with the polarized case. Whereas the unpolarized scal- 
ing functions Fi_2(3:^) have a simple interpretation in the Naive Parton Model (where the 
nucleon is considered as an ensemble of essentially free massless partons) and a straight- 
forward generalisation in the framework of perturbative QCD, the spin dependent scaling 
functions gi^2{x) are much more subtle, each fascinating in its own way. The function 
gi{x) which, at first sight, seems trivial to deal with in the Naive Parton Model, turns 
out, within perturbative QCD, to have an anomalous gluon contribution associated with 
it. In addition the first moment of g\{x) turns out to be connected with an essentially 
non-perturbative aspect of QCD, the axial ghost which is invoked to resolve the U{1) 
problem of the mass of the r]' . And g2{x) turns out not to have any interpretation at all 
in purely partonic language. 

What is also fascinating is the extraordinary interplay of theory and experiment in 
the study of gi{x). For a long time the theory of gi{x) remained comfortably at the level of 
the Naive Parton Model. Then, in 1988, came the disturbing results of the European Muon 
Collaboration (EMC) [ASH 88, 89], which differed significantly from the naive theoretical 
predictions. These results could be argued to imply that the sum of the spins carried by 
the quarks in a proton (S'^;'"^'''^^) was consistent with zero, rather than with 1/2 as given 
in the quark model, suggesting a "spin crisis in the parton model" [LEA 88]. This led to 
an intense scrutiny of the basis of the theoretical calculation of gi{x) and the discovery of 
the anomalous gluon contribution [EFR 88]. (As so often happens in theoretical physics 
it turns out that such an effect had already been studied to some extent in a largely 
overlooked paper of 1982 [LAM 82]). So surprising was this discovery that the calculation 
was immediately checked by three different groups [ALT 88; CAR 88; LEA 88a] who all 
arrived at the same result. (Somewhat fortuitously, as it turns out, as was demonstrated 
in [CAR 88]). 

In the modified theoretical picture, the quantity AE = 2Sf^"-'^^, whose value had to 



be consistent with zero as a consequence of the EMC experiment, is replaced by the hnear 
combination (for 3 flavours) AS — (3Q;s/27r) Ag^ where /S.g — Jq dx Ag(x) and /^g{x) is 
the polarized gluon number density. 

It has to be stressed that as a consequence of QCD a measurement of the first 
moment of gi{x) does not measure the 2;-component of the sum of the quark spins. It 
measures only the superposition 

1 r . ^ 3as . 
- AE -Ag 

and this linear combination can be made small by a cancellation between quark and gluon 
contributions. Thus the EMC results ceases to imply that AE is small. 

The function g2{x), on the other hand, docs not have any simple interpretation in 
the Naive Parton Model and it is a triumph of pcrturbativc QCD that one can derive a 
sensible, gauge invariant result for it in the QCD Field Theoretic Model [EFR 84]. 

In this review we concentrate almost exclusively on the polarized case. A good 
survey of the unpolarized case can be found in [ALT 82] or, at a more introductory level, 
in [LEA 94] . Our treatment attempts to cover both the theory and the phenomenology of 
the subject and is meant to be reasonably self-contained. 

On the phenomenological side we outline the formalism for discussing DIS in terms of 
one photon exchange. We briefly mention the extension to weak interactions (both neutral 
and charged current) initiated by charged leptons. We introduce the scaling functions 
gi,2ix) and consider how data on them can be extracted from experiments on cross-section 
asymmetries using both longitudinal and transversely polarized targets. We consider both 
hydrogen and nuclear targets and we also address the problem of how to extract data at 
fixed over the whole range of x, as needed for the testing of sum rules. 

Since the data taken in different x-bins correspond to different ranges of Q^, the 
production of a data set for all x at some fixed Q'^ is a non trivial "experimental" problem. 
There is some evidence that the Compton asymmetry Ai{x,Q^) is almost independent 
of and the extraction of gi{x,Q^) at the desired is carried out by assuming Ai 
independent of Q"^. Now it is reasonable to believe that Ai{x,Q^) will vary only slowly 
with Q^, but there is no justification for believing it to be exactly constant. Nonetheless 
it is experimental practice NOT to include an error in gi{x^ Q"^) to refiect the uncertainty 
in the variation of Ai with Q^. This is, to say the least, quite unjustifiable. A possible 
strategy for carrying out the extrapolation in is discussed in Section 2.1.3 d). 

An extensive phenomenological analysis of all the data on gi(x,Q'^) and its first 
moment ri(Q^) is given including the very recent data from experiments on deuterium 
and ^He targets. The implications of the data are considered in the framework of the Naive 
Parton Model and of the Operator Product Expansion (OPE). Aspects discussed include: 
information from hyperon /9-decay; the extrapolation to a; = 0; higher twist effects. 

We argue, in agreement with several recent analyses, that there is absolutely no 
evidence for a failure of the very fundamental Bjorken sum rule. 

On the theoretical side we distinguish between three approaches of varying degrees 
of sophistication: 

i) In the Naive Parton Model the nucleon is viewed as an assembly of free, parallel 
moving quarks. The treatment is probabilistic and its essence is summarized by the relation 
connecting the virtual photon-hadron cross-section a'^*^ to a convolution of the 7*-quark 



cross-sections and the number density of quarks in the hadron, fq/H{x)' 

<7^*^ = J dx'a^\x/x')f,/H{x'). 

ii) In the QCD Improved Parton Model one supplements the Naive Parton Model 
with QCD-controUed, probabihstic, Altarelli-Parisi splitting functions. The description 
remains probabilistic and the main effect is the replacement of number densities fq/H{x) 
by Q^-dependent densities fg/H{x,Q'^). 

iii) In the QCD Field Theoretic Model one starts with genuine QCD Feynman 
diagrams. We call it a "Model" because the diagrams have to be divided into a hard 
scattering part, which is calculated strictly according to perturbative QCD, and a soft part 
which we cannot calculate and which has to be modelled. One makes certain reasonable 
sounding assumptions about the behaviour of this non-perturbative part. The Model is, 
in general, not probabilistic and involves hadronic matrix elements of products of quark 
and gluon interacting fields at different space-time points. In those cases where the result 
is probabilistic one can usually recover the Naive Parton Model upon treating the fields 
as free fields, or, with a less drastic approximation, the QCD Improved Parton Model. 

The Operator Product Expansion (which deals with matrix elements of local op- 
erators), where it can be justified leads to results perfectly compatible with the Field 
Theoretic Model, but the latter has a much wider field of applicability and can be utilised 
in situations where use of the OPE cannot be justified. 

The calculation of gi,2{x) in the Naive Parton Model is covered in some detail and 
the source of the difficulties with g2{x) is, we hope, very clearly delineated. We introduce 
and discuss the general structure of the Operator Product Expansion and give expressions 
for the n-th moments of gi,2{x, Q"^). We are careful to state the range of n for which the 
results can be proved and this leads to an examination of various claims in the literature 
about relationships between gi and (^2) about the first moment of g2{x) (/ dx g2{x) = 0: 
the Burkhardt-Cottingham sum rule) and about the second moment of gi{x) + 2g2{x) 
(/ dx x[gi{x) + 2g2{x)] ~ 0: the Efremov-Leader-Teryaev sum rule). 

We explain that neither of these sum rules follows directly from the OPE and 
that their validity rests upon an assumption about the invertability of certain Fourier 
transforms (cquivalently, about the behaviour as a; 0) and that this invertability does 
not occur if g2{x) has the small x behaviour l/x"^ suggested by Regge Theory, so that the 
above integrals diverge. 

The axial anomaly and the gluon axial current are discussed in some detail. We 
study the connection between K^^ and the gluon spin operator and derive the anomalous 
gluon contribution to gi{x). There is some disagreement in the literature regarding the 
kernel relating the polarized gluon number density ^g{x) to gi{x), and we give arguments 
in favour of a particular choice. 

There are also arguments in the literature against the whole concept of an anoma- 
lous gluon contribution, on the grounds that K^^ is not a gauge invariant operator. Here 
we would hke to remind readers that in gauge theories (even in the Abehan case of QED) 
it often happens that the operators representing fundamental physical quantities like mo- 
mentum are not gauge-invariant. This lack of gauge invariance is innocuous and irrelevant 
provided that the expectation values of those operators between physical states are gauge 
invariant. Ultimately, the gauge invariance of /S.g[x) is guaranted by showing how it is 
related to the physically measurable cross-section for 2-jet production. 

There are also claims in the literature that the anomalous contribution to gi vanishes 



if the calculation is carried out with strictly massless gluons. We point out that this is 
irrelevant. In the QCD Field Theoretic Model there is, in principle, an integral over the fc^ 
of the gluon and the point A;^ = is of zero measure. Only in a Parton Model calculation 
with gluon mass fixed at A;^ = would there be any problem. 

The whole question of the calculation of gi^2 and of the anomaly is finally approached 
from the point of view of the QCD Field Theoretic Model. Contrary to some assertions 
in the literature we show that there is perfect consistency between this approach and the 
Operator Product Expansion. 

We contrast the calculation of g2{x) and gi{x) in QCD Field Theoretic Model, 
showing that the Feynman diagrams which have a simple quark-parton interpretation lead 
to an unacceptable result for g2{x), i.e. a result which does not respect gauge invariance. A 
detailed demonstration is given of how the unwanted gauge dependent terms are cancelled 
when one includes diagrams corresponding to quark-gluon correlations in the nucleon, a 
concept which is outside the framework of the Parton Model. 

We provide a brief introduction to non-perturbativc and topological effects in QCD. 
to the U{1) problem, to the axial ghost and to the isosinglet. Generalized Goldberger- 
Treiman relation (GGT). A connection is established between the total quark spin and 
the coupling of the rj' to nucleons and this suggests that the sum of the quark spins is, in 
fact, close to the canonical value of 1/2, in agreement with the interpretation of the EMC 
data utilising an anonalous gluon contribution to gi{x). 

A dominant feature throughout is the polarized gluon distribution Ag{x). It is 
essential that this important element in the internal structure of the nucleon should be 
measured independently, in as many other reactions as possible. We end with a brief 
survey of reactions in which this might be feasible. 

Finally, a word about notation. It may seem pedantic to squabble about choices 
of notation. But when a subject threatens to become incomprehensible as a result of an 
injection of illogical, confusing notation, it is necessary to draw attention to the danger. 
Our Section on notation, consequently, will comment on this issue. 

1.1 Notation 

The Field Theoretic approach to hard scattering was developed in a pioneering 
paper by Ralston and Soper in 1979 [RAL 79]. They introduced a precise, logical, clear 
and informative notation for the hadronic matrix elements of the bifocal operators that 
appear. 

On the other hand the spin dependent scaling functions gi,2{x) were introduced by 
Bjorken in 1966 [BJO 66]. They too have a precise and clear meaning in the expression 
for the measured spin-dependent cross-section. 

A complete relabelling of the Ralston-Soper matrix elements [JAF 91] ought to 
require serious justification (especially the use of gi for one of them when it is only the 
leading twist approximation to Bjorken gi) and such justification has never been offered. 
We therefore follow the Ralston-Soper notation. 

Concerning 7 matrices and spinors we adopt the conventions of Bjorken and Drell 
[BJO 65] except that our spinors are normalised so that for a fermion with four-momentum 
P'' = {E,P) 

u^u = 2E, v^v = 2E, uu = 2M, vv = -2M 
for both massive and massless particles. 



One has then 

u{p)Yu{p) = 2jf v{p)Yv{p) = 2p^ . 

We use 

1 2 

75 7 7 7 7 

and our convariant spin vectors S'^ are normahsed to S^S^ = —1 . 
For a fermion of mass M, then, 

u{p, S)riMP, S) = -v{p, S)r^5v{p, S) = . 

For massless fermions of hehcity A = ±1/2, or of chirahty 2A = ± 1, 

u{p, A)7^5«(P, A) = -vip, X)Y'iMp, A) = hm 2M5^(A) = 4A 

All states are normalized so that 

{P'\P) = {2nf 2E S%p' -p). 

We use 

2^0123 — +1 ■ 

Finally, because 'transversity' has a quite specific meaning [DAL 56; KOT 70] in 
terms of a phase convention etc., we avoid its use and simply talk of longitudinal and 
transverse polarization states. 

2 Polarized lepton-nucleon Deep Inelastic Scattering - General aspect 

2.1 General formalism in one photon exchange approximation 

Consider the inelastic scattering of polarized leptons on polarized nucleons. We 
denote by m the lepton mass, k [k') the initial (final) lepton four-momentum and s [s') 
its covariant spin four- vector, such that s-k — Q {s' -k' — and s- s — —1 {s' ■ s' — —1); the 
nucleon mass is M and the nucleon four-momentum and spin four-vector are, respectively, 
P and S. Assuming one photon exchange, see Fig. 2.1, the differential cross-section for 
detecting the final polarized lepton in the solid angle dVl and in the final energy range 
(E', E' + dE') in the laboratory frame, P = (M,0), k = {E,k), k' = {E\k'), can be 
written as [LEA 85] 

^^ = ^-L W^'^ (2 11) 

dn dE' 2Mq^ E ' ^ ^ 

where q — k — k' and a is the fine structure constant. 
In Eq. (2.1.1) the leptonic tensor L^j^ is given by 

V(A;, s; k\ s') = [u{k\ s') 7^ u{k, s)]* [u{k' , s') 7, u{k, s)] (2.1.2) 

and can be split into symmetric {S) and antisymmetric (A) parts under /x, v interchange: 



L^,{k,s-k',s') = L^^ {k;k')-riLfJ {k,s;k') (2.1.3) 
+ L' f ) {k, s; k', s') + ihX^ (k; k', s') 



where 



LlV (fc; k') 


~ k^k^ + k^ky g^u {k ■ k 


- m^) 


(2.1.4) 


Lil) {k, s; k') 


= m Ef^^a/d -s" {k - k'Y 




(2.1.5) 




= {k ■ s') {k'^s^ + s^kl - g^^ 


k'-s) 






- {k-k' - w?) + ^'fji^y 


- g^iu s- s') 


(2.1.6) 




+ {k' ■ s){s'^K + k^s',) - {s ■ 


s ) {k^kjj + k^jky) 




(k;k',s') 


= m e^yai3 s"^{k-k'Y . 




(2.1.7) 



On summing Eq. (2.1.3) over s' and averaging over s one recovers the usual unpo- 
larized leptonic tensor, 2L^^J. On summing over s', one obtains 21^^^^ + 2iLfJ . 

The unknown hadronic tensor W^v is similarly defined in terms of four structure 
functions as [DRE 64; CAR 72; HEY 72]. 



with 



2M 



W,M. P. S) = Wll\q; P) + I Wlf{q; P, S) 



(2.1.8) 



+ 



P. / % 



Pf o 



W^iP-q.q^ 



M2 



:2.i.9) 



2M 



< MS^Gi{P ■ q, q^) 
+ [ip.q)S^-iS.q)P^]^^^^^y 



(2.1.10) 



Prom Eqs. (2.1.1, 3, 8) one has 



(Pa 
dQ dE' 



a 



E' 



2M54 E 



(2.1.11) 



The individual terms in square brackets can be separately studied by considering 
cross-sections or differences between cross-sections with particular initial and final polar- 
izations [ANS 79]. Each of these terms is, at least in principle, a measurable quantity 
which is either a function of the two spin-averaged structure functions W\ and 14^2 (terms 
containing W^^') or of the two spin-dependent structure functions Gi and G2 (terms 
containing W'^). For example, the usual unpolarized cross-section is proportional to 



^2 unp 1 



d^a 



dn dE' 



4 dn dE' 



- (k, s, P, S; k', s') 



a 



E' 



2Mg4 E 



2^(5) ^,.iS) ^ 



(2.1.12) 



while differences of cross-sections with opposite target spins single out the L^^^ 
term: 



E 



d'^cr dj^G 

- [k, s, P, —S; k , s) — s, P, S;k ,s) 



dn dE' 



E' 



(2.1.13) 



2.1.1 Structure functions and Bjorken scaling 

The cross-section for the inelastic scattering of unpolarizcd leptons on unpolarized 
nucleons, in the laboratory frame, can be written explicitely using Eqs. (2.1.4, 9, 12) and 
neglecting the lepton mass, as 



dn dE' 



9 6 
2Wi sin^ - + W2 cos^ - 



(2.1.14) 



where 6 is the LAB scattering angle of the lepton. Its measurement supplies information 
on the unpolarized structure functions Wi{P ■ q, q^) and W2{P • q, q^). 
In the Bjorken hmit, or Deep Inelastic Scattering (DIS) regime. 



00 



V — E — E' ^ 00 X — 
they are known to approximately scale: 



where Fi^2 vary very slowly with at fixed x. 

Similarly, from Eqs. (2.1.5, 10, 13), one has 



2P ■ q 2Mu 



, fixed (2.1.15) 



lim MPFi(P-g,g2) = Fi{x) 

Bj 



lim iyW2iP-q,Q^) = H^) , 

Bj 



E 



(2.1.16) 



d'a 



q^E 



dn dE' dn dE' 



{q ■ S){q ■ s) + Q\s ■ S) MG, + {s ■ S){P ■ q) - {q ■ S){P ■ s) 



(2.1.17) 

G2 
M 



which yields information on the polarized structure functions Gi{P-q, (f) and 6*2 (P-g, q^). 
They also, in the Bjorken limit, are expected to scale approximately: 



hm ^^^G,{P.q,Q') = g,{x) 

Bj V 

lim v{P-q) G2{P-q,q^) = g2{x) . 
Bj 



In terms of gi^2 the expression (2.1.10) for Wjj^^ becomes 



(2.1.18) 



Wlt\q:P.s) = e,,^p q4sPg^{x,Q') + 



, {S • q) P^ - 



g2{x,Q^)\. (2.1.19) 



2.1.2 Cross-section differences 

Aiming at gathering information on the polarized structure functions Gi and G2 
we look at Eq. (2.1.17) and specialize it to particular spin configurations of the incoming 
leptons and target nucleons. 

Let us consider first the case of longitudinally polarized leptons, that is initial 
leptons with spin along (^) or opposite (^) the direction of motion, while the nucleons 
at rest are polarized along {S) or opposite {—S) an arbitrary direction S. We then have 

st, = -st^ — (\k\, kE) fc=T^ (2.1.20) 
m \k\ 

= (0,5). 

We take the 2;-axis along the incoming lepton direction and define (see Fig. 2.2) 

k" = (E,0,0, |fc|) ~ ^(1,0,0,1) 

k'f" = {E',k')c^E'{l,k) 

— E'{l,sm9cos(p, sin^sint/?, cos^) (2.1.21) 

S — (sin q; cos /3, sin o; sin /3, cos a). 

Using Eqs. (2.1.20, 21) into Eq. (2.1.17) one obtains (at leading order in m/E) 

d^a^'^ d^a^'-^ Aa" E' 



dfl dE' dfl dE' E 

X \[Ecosa + E' cos O] MGi + 2EE' [cos O - cos a] i . 



(2.1.22) 



a is the polar angle of the nucleon spin direction, i.e. the angle between k and S, and © 
is the angle between the outgoing lepton direction, k , and S : 

cos G = sin 6^ cos ip sin a cos {3 

+ sin 6* sin V? sin a sin /9 + cos ^ cos a (2.1.23) 
= sin 9 sin a cos + cos 6 cos a 

where (f) — P — (pis the azimuthal angle between the {k,k) scattering plane and the {k, S) 
polarization plane (see Figs. 2.2 and 2.3). 

For particular values of a one recovers familiar results [CAR 72; HEY 72]. For 
nucleons polarized along {=>) or opposite (<^) the initial lepton direction of motion one 
has a = 0, e = 6, and Eq. (2.1.22) gives 



4^2 E' 



dfl dE' dfl dE' 



g2 E 



{E + E' cos e)MGi-Q^G2 



(2.1.24) 



If the nucleons are transversely polarized, that is the nucleon spin is perpendicular 
to the direction of the incoming lepton, then a — 7r/2 and Eqs. (2.1.22, 23) yield 



W E'^ 



dfl dE' dfl dE' 



E 



sin ^ cos {MG1 + 2EG2). 



(2.1.25) 



Notice that if the nucleon spin is perpendicular to the scattering plane {a — (f) — 
7r/2), then the cross-section difference in Eq. (2.1.25) is zero; such a difference has its max- 
imum absolute value when = or tt, that is when the nucleon spin vector, perpendicular 
to k, lies in the scattering plane. 

In the above we considered longitudinally polarized leptons. Note that it is im- 
practical to deal with transversely polarized leptons. Because, for transversely polarized 
incoming leptons, that is leptons with spin perpendicular to the direction of motion, we 
have 

s = {0,s), (2.1.26) 

with the unit vector s orthogonal to fc, s - k = 0. In such a case, contrary to what happens 
for longitudinally polarized leptons, Eq. (2.1.20), there is no factor E/m to cancel the 
factor m/E which appears in the cross-section differences (2.1.17), and the latter turn 
out to be vanishingly small in the large energy limit {m/E — > 0). 

Useful information on the unpolarized structure functions Wi and W2 could be 
obtained by looking at initial- final lepton spin asymmetries [ANS 79] ; such a technique, 
however, requires the measurement of the scattered lepton polarization, which is probably 
too difficult to contemplate at present. 

2.1.3 Experimental measurement of gi and g2 on nucleon targets 

Measurements of cross-section differences with particular lepton and nucleon spin 
configurations provide information on the polarized structure functions Gi and G2 or on 
the scaling functions gi and g2, Eq. (2.1.18). 

A single difference measurement, however, only yields information on a combination 
of Gi and G2, rather than on the separate structure functions. Extracting from the data 
values of Gi or G2 alone still requires some further approximate procedure. 

a) Longitudinally polarized target 

In the performed experiments [ALG 78; BAU 83, 88; ASH 88, 89; ADE 93, 94; ANT 
93] the quantity actually measured is the longitudinal spin-spin asymmetry in ip iX, 



da —da . _ 

All = ' (2.1.27) 

da^ + da^ 

where da is short for d^a/{dQ. dE') and the denominator is simply twice the unpolarized 
cross-sections. Prom Eqs. (2.1.14) and (2.1.24) we derive 

Q'[{E + E'cosd)MG^-Q^G2] ^2.1.28) 



2EE' [2Wi sm\e/2) + W2 cos^{e/2)] 



For reasons to be explained shortly the asymmetry Ay is usually expressed in terms 
of virtual Compton scattering asymmetries ^41^2 [LEA 85] , 

A\i^D{Ai+r]A2), (2.1.29) 

where D and rj are known coefficients given in Appendix A. The analysis of the data 
then proceeds through subsequent approximate steps, explained in detail in Appendix A, 
which lead to the expressions 

All ~ (2.1.30) 



and 

where ^2(2:) is the unpolarized scahng structure function, Eq. (2.1.16), and R is the ratio 
of the longitudinal to transverse cross-section. 

All the approximations involved in the above simplifications can be shown to be 
quite harmless when measuring gi [ASH 88, 89]. Indeed not only is rj small in (2.1.29) but 
one can show that 

IA2I < VR (2.1.33) 

and R is known to be small. 

If the goal is simply the knowledge of gi and g2 on an event-by-event basis then 
there is a somewhat more direct way of utilising the data. 

From (2.1.14, 18 and 28) we have 

MuQ^E , 2xM 

A|| =g,- — — g2 (2.1.34) 



2a^E'{E + E' cos 9) dQdE' " E + E' cos 6 

which can be written more concisely as 

g^ - Kg2 = 2K dfj""*' A|| (2.1.35) 

with 

2xM xM 

K — 
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where 

acos(^/2) 



CMott 



_2Esm^{e/2) 

The RHS of Eq. (2.1.35) is constructed directly from experiment, with no need of 
further data analyses in order to extract F2 and R, as required in Eq. (2.1.31). 

As we said previously the single measurement of A\\ (and dcr"'*^) gives us information 
on the combination gi — Kg2, rather than on gi or g2 alone. The usual argument [ALG 78; 
BAU 83, 88], however, is that the g2 term in Eq. (2.1.35) can be safely neglected because 
of the kinematical coefficient k, which, in the large energy limit, is very small, as can be 
seen from Eqs. (2.1.36). This is indeed confirmed by a more careful analysis of the g2 term 
[LEA 88]. We can thus conclude that the measurement of the quantities on the RHS of 
Eq. (2.1.35) provides essentially a direct measurement of the polarized scaling structure 
function gi. 

If one requires data on gi{x^Q'^) on an event-by-event basis, (2.1.35) is the most 
direct way to use the experimental data on A\\. However, it is of great interest to assemble 
results for gi{x, Q"^) over the entire x-range Q < x <1 at the same Q^. Experimentally this 
is impossible. The kinematics of the experiment dictates that smaller x will correspond 
to a smaller range of accessible Q^. 



The experimentalist is thus forced to extrapolate the data in at fixed x, and 
the question is then: which quantities vary most smoothly and slowly in Q^l It is an 
experimental fact that, where it has been studied, A\\{x^Q'^)/D varies only slowly with 
Q^. (This question is discussed further in Section 4.5). For this reason experimentalists 
prefer to express their measurements in terms of data on Ay (a;, Q"^) via (2.1.31). It is then 
assumed that A\\{x,Q'^)/D is essentially independent of and the value of gi{x, (Q^)) 
quoted, for an experiment with mean equal to (Q^) , is really 

It is important to bear this in mind when testing sum rules which hold at fixed Q^. 
A possible strategy for improving the extrapolation [EFR 94a] is discussed in Section 
2.1.3 d). 

As mentioned above, all the approximations leading to (2.1.31) are safe if one is 
trying to evaluate gi{x). In the next Section we consider the perpendicular asymmetry 
which is used to measure a combination of gi and g2 with the aim of extracting g2. 
Since g2 is expected to be much smaller than g\ some care may be necessary in utilising 
an approximate version of gi{x). 

b) Transversely polarized target 

If we want to have information on g2 we must then consider other spin-spin asym- 
metries. By scattering longitudinally polarized leptons on transversely polarized nucleons, 
one can measure the quantity 

A± = - — 2.1.38 

da^^ + da^^ ^ ' 

where, again, da is short for ^al{d£l dE') and the denominator is twice the unpolarized 
cross-section da'^'^^. 

Prom Eqs. (2.1.14) and (2.1.25) one obtains 

. Q'^sme{MGi + 2EG2) 



2E [2Wi sin2(^/2) + W2 cos'^{9/2)] 



cos(/) (2.1.39) 



where 4> is the difference between the azimuthal angles of S and , (j) = f3 — if (see 
Figs. 2.2, 3). For sake of clarity we show in Fig. 2.4 the spin configuration leading to Eq. 
(2.1.39), for = 0. In the following we shall always assume = 0; if necessary one could 
integrate Eq. (2.1.39) over the range covered by the experimental apparatus. 

One can repeat for g2 the same direct procedure followed for A\\, Eqs. (2.1.34-36), 
and can use Eqs. (2.1.14, 18 and 39) (with = 0) to write 

9^ + i^9i={^^K' da--^ A^ (2.1.40) 

where 

^, ^ Q^EMu ^ Ecos\0/2) ^ 
Aa^E''^ sin 9 2xaMott sin 61 
and (TMott was given in Eqs. (2.1.36). 

A measurement of A^, i.e. the RHS of (2.1.40), provides direct information on the 
structure function combination g2 + u/(2E) gi. 

In this case the coefficient of gi is not negligible and to isolate g2 one must feed in 
one's knowledge of g^ obtained from the A\\ measurement. 



The A± type experiments are much harder than the A\\ type so it is perhaps worth 
noting that one can obtain an estimate of g2{x) from the A\\ experiments if performed at 
several different beam energies [LEA 88]. 

At fixed values of and x the only dependence on the beam energy E in the LHS 
of (2.1.35) is in the coefficient of g2{x). 

A study of the energy variation of the RHS of (2.1.34) thus allows one, in principle, 
to learn about g2ix). Unfortunately, because the coefficient of g2{x) is so small, it is not 
clear whether this method is practicable or not. 

In a similar way the measurement of A± at different beam energies E, allows the 
isolation of g2{x) from (2.1.40). For example, measurements at E = Ei and E2 yield 

{El - E2) g2{x) = [uK' da^'^P Aj_]e=e, - Wk' da^''^ Aj_]e=e, ■ (2.1.42) 

Equations (2.1.40, 41) provide the most direct access to g2 on an event-by-event 
basis, provided one's knowledge of gi is accurate enough. 

However, as with gi, if data is desired at fixed over the entire range of x a 
different strategy is required which we now describe. 

c) Combined analysis using Ay and A±_ 

From the longitudinal polarization data one can extrapolate A\\/D smoothly in 
at fixed X. Combining this with the perpendicular polarization data where it is measured 
one can construct (see Appendix A for definitions of kinematical coefficients) 

A'-'^X- ^ fr^" I (2143) 

at the values of the A± experiment. But A' should itself vary slowly with Q"^ because one 
has 

= ,,.,f!^^g'<g.)| A'M^) (2.1.44) 

SO that extrapolation in at fixed x should be smooth. Via (2.1.44) one can thus estimate 
gi + g2 over the entire x, range. 

Again, with A' extrapolated to the relevant Q^, one can obtain an improved evalu- 
ation of gi{x, Q'^) via 

^'^^'^^-2x[l + R{x,Q^)] l + 4M^xyQ^ [d VQ'V VQ" J i ' ^ ^ ^ 

The formulae (2.1.44) and (2.1.45) are exact. No approximations have been made. 
And they are expressed in terms of the functions A'(x, Q^) and A\\/D which should both 
be slowly varying functions of at fixed x. 

d) The problem of extrapolating in 
We mentioned in a) above that 

A(:r,Q^) = ^M|^ (2.1.46) 

is taken to be independent of in the experimental evaluation of gi{x, {Q^)) via (2.1.37). 
There is no allowance made for the error in extrapolating from the measured region of 
for the X involved to the required value (Q^)- A simple linear parametrization a + bQ^ or 
a + blnQ^ for Ai will not work because a best fit will yield a very small value of 6, but 
with very large errors, leading to unrealistic error estimates on gi{x, (Q^))- 



There is no rigorous theoretical solution to this experimental problem, but the 
following approximate procedure should lead to an improved estimate of gi{x, {Q"^)) and 
its error [EFR 94a] . 

Let xi, X2, ■ ■ ■ be the values of x for the data bins and let Ql, Q2, ... be the mean 
values of for each x-hin. 

We define a zeroth order approximation to gi{x, (Q^)) for each Xi via (2.1.37), i.e. 

g^^\x,, {Q')) ^ A^{x,,Q',)F^{x,, {Q')) (2.1.47) 

where we have written 

Fi = — - (2.1.48) 
for brevity. We define an improved estimate of gi{x, {Q'^)) for each Xi via 

g,{x,, {Q')) ^ A,{x,, {Q'))F,{x,, {Q')) (2.1.49) 
where Ai{xi, (Q^)) is obtained using 

A,{x,,{Q'))^A,{xi,Q',) + b{xi,Q',) H{Q')/Ql) (2-1-50) 



and b{xi, Qj) is estimated from 

b{xi,Q',) 



dgi{xi,Q'^) 



Fi(x,,g|) 91ng2 
A^ixi,Ql) dFiixi,Q^) 



(2.1.51) 



Fi(x,,gf) 9ing2 



All the terms on the RHS of (2.1.51) are known from experiment except for the 
derivative of gi. The latter can be approximately calculated from the evolution equation 
provided we use the experimental fact that the flavour singlet part of gi is much smaller 
than the non-singlet part [see (4.2.22)]. Thus we treat the evolution as if it were purely 
non-singlet, i.e. we compute 



dgi{xi,Q^^ 



dlnQ^ 
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, f^gi{y,Qh^P,J-] , (2.1.52) 



where APqg is the non-singlet polarized splitting function. In fact, for the non-singlet case, 
polarized and unpolarized splitting functions are equal [ALT 77] , 

AP,,{x) = P,,{x) = \ (^)^ • (2-1-53) 

Finally, in the convolution integral in (2.1.52), we approximate gi by its known 
zeroth order estimate (2.1.47). Thus we suggest the approximate formula 



dgi{xi,Q'^ 



ainQ2 



^CjS?Hy,Qf)ApJ^). (2.1.54) 



This will, of course, provide an estimate for the value of b in (2.1.51) and for the 
error on it. 



2.1.4 Measurement of on nuclear targets 

We consider spin 1 targets, like deuterium, or spin 1/2 targets like ^He. The asym- 
metries considered are the analogues of (2.1.27 and 38) and are defined for nucleus A 

by 

Af = "^^^ ""^^-^ (2.1.55) 
, A da~7^ — da~7^ 



da^ + dcTj^ 
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(2.1.56) 



where a^, aX means = ±1/2 for a spin 1/2 target, but = ±1 for a spin one target 
longitudinally polarized; and similarly for transverse polarization. 

We deal only with the case where the constituents of the nucleons are assumed to 
contribute independently to the scattering. This means we neglect shadowing and Fermi 
motion. 

For unpolarized scattering on deuterium this is tantamount to taking 

CTd^CTp + CTn, (2.1.57) 

a perfectly reasonable approximation for but one which can be misleading when dif- 
ferences of cross-section are being studied. For example the Gottfried sum rule requires 
the combination ap — an for which the approximation 2ap — aa may be dangerous because 
the subtraction of comparable quantities magnifies the error [EPE 93] . 

In a similar way, for the testing of the Bjorken sum rule we shall require the quantity 
9i ~ o'^^ might worry about obtaining from nuclear data on the basis of 

independent scattering. Happily because gf and g^ are expected to be quite different in 
magnitude this should be quite reliable [EPE 92]. 

On the basis of independent scattering on the Z protons and N neutrons in the 
nucleus the nuclear cross-section difference is given by 

da^ - daX = 2 [Zdap Vp A\ + Ndan Vn A[f] , (2.1.58) 

where dap^n ai'c the unpolarized nucleon cross-sections, A^^''^ the nucleon longitudinal asym- 
metries and Vp^n the longitudinal polarization of the nucleons in the nuclear state with 
J, = 1/2 or 1. ' 

For the asymmetry defined in (2.1.55) one then has 

Af = fpVp + fnVn (2.1.59) 

where 

f = ^ f = ^ (2 1 60) 

Zdap + Ndan Zdap + Ndan v ■ ■ ; 

are the fractions of events originating on protons and neutrons respectively. 
An analogous result holds for Aj_. 

Equation (2.1.59) is the fundamental formula which should be used to extract 
from A^ and a knowledge of Ay. This is especially true when one does not have data on 
A-^ and Ajj at the same (Q^) because it appears that the asymmetries themselves show 
very little dependence. 



For deuterium, because of the D-state admixture one has 

3 



= = (1 - r (2.1.61) 



where ud — 0.058 is the D-state probabihty. 

For the fractions fp^n one has, approximately. 



Fi/{l + RP) F^/{l + R-) 



where F2 is the deuteron F2 per nucleon. 

Often one defines gf per nucleon via (2.1.30, 31) 



so that (2.1.59) becomes 



gt M') = ^^-^^|^?(^,Q^) + ■ (2.1.64) 

In the above (2.1.55, 56) we have defined asymmetries for 100% polarized targets. 
For spin 1/2 targets with degree of longitudinal polarization V the generalization 
of these equations is quite straightforward: 

_ 1 l da-{-V)-d<j-{V) \_ 1 da-{-V)-da-{V) 

- V\da-{-V) + da-{V)r V 2da ^^"^-^^^ 

where da is the Tinpolarized cross- sect ion. A similar result holds for A^. 

For spin 1 targets the result is a little more subtle. If p_, po are the probabilities 

of finding states with = 1, — 1, in the target, then the degree of polarization is 
[BOU 80] 

V^p+-p- (2.1.66) 

and the alignment is 

A^l-3po. (2.1.67) 

Then one has 

da-^{-V) - da^iV) = V{da^ - da^ } , (2.1.68) 

but 

A 

da^{-V) + da^{V) = 2da + —[da+ + da_ - 2 dao] (2.1.69) 

where da is the unpolarized cross- sect ion for the spin 1 target. 
It follows that for a spin 1 target 

, ^ 1 da-(-V) - da-jV) 

" V 2da+{A/^)[da+ + da^-2daQ] K ■ ■ ) 



and this is strictly only equal to 



1 da^{-V) -da^{V) 
V 2da 



(2.1.71) 



if the alignment A is known to be zero. 

Of course for the 100% polarized target corresponding to (2.1.55), Po — and the 
alignment has its maximum value A = 1. However, even for A ^ 0, the term in square 
brackets in (2.1.69, 70) vanishes in the approximation of independent scattering from the 
nucleons. Since the results (2.1.59, 60) are based on this approximation it is consistent in 
these equations to ignore the term proportional to A. 

It is worth noting that one could, in principle, check the assumption of independent 
scattering by experimentally testing (2.1.69). 

For ^He the wave-function is almost entirely an S'-state with the two protons having 
opposite spins. Thus all the spin is carried by the neutron but there is some mixing in the 
wave-function [WOL 89; CIO 93] and one estimates 

V'f^ = (87 ± 2)% V'/^ = (-2.5 ± 0.3)% . (2.1.72) 

For the fractions fp^n one takes the approximation 

where ^2*^^ is the helium-3 F2 per nucleon. 
Defining g^^^ per nucleon via (2.1.31), 



one obtains 



gl"" (x, g^) = \ [(0.87 ± 0.02) ^^(x, Q^) - (0.050 ± 0.006) gl{x, Q^)] . (2.1.75) 

Note that the above is again based on completely independent scattering from the 
constituent nucleons. This might be reasonable at high but it is not likely to be a 
good approximation in the E142 experiment which has (Q^) = 2 (GeV/c)^ and which 
will be discussed in Section 4.5. Unfortunately we do not have any simple prescription for 
improving the analysis. 



2.2 Neutral and chcirged current weak interactions initiated by chcirged 
leptons 

Consider now the deep inelastic scattering of longitudinally polarized leptons on 
polarized nucleons at very high energies, taking into account weak interactions of both the 
neutral and charged current type. We do not study neutrino initiated processes, although 
the extension of our results to such a case would be trivial, because of the technical 
difficulties in polarizing the large nucleon targets needed for neutrino scattering, which 
make such experiments quite unrealistic. Useful information can be found in [NAS 71; 
DER 73; NIK 73; AHM 76; KAU 77; JOS 77; CAH 78; BAR 79; CRA 83; LAM 89; VOG 
91; JEN 91; RAV 92; MAT 92, 92b; XIA 93; DEE 93]. 



Let us start from neutral current processes, IN — > IX, which, at lowest perturbative 
order, proceed via the exchange of one photon or one boson. Figs. 2.1, 5a. The cross- 
section then receives contributions from a purely electromagnetic, a purely weak and an 
interference term, so that Eq. (2.1.1) becomes 

(Pa E' 

- % E (2.2.1) 



dE> 2M54 ^ 



On summing over the final lepton spins and neglecting terms proportional to mjE 
or m/E', the leptonic tensors L^^ are given, for negatively charged leptons (e~,/x~), by 

LI. = E [^(^'' 7m <K A)]* [u{k\ s') 7. u{K A)] 

= 2 + - k ■ k'g^, - 2i\ e^.^p k'^k'"] (2.2.2) 
Ll^ = E [^(^'' 7m(5v - 5.75) u{k, A)]* [«(^', s') 7. «(^, A)] 

= (^,-2A^Jl;^, (2.2.3) 
^Ji. = E ['^(^'' 7/.(^^ - 9a15) u{k, A)]* ['iZ(A;', s') ^.{g^ - c/^75) w(A;, A)] 

s' 

= (^, - 2A^J2 L^;, (2.2.4) 

where 2A = ±1 denotes twice the helicity of the initial lepton. L^^ agrees with Eqs. (2.1.3- 
5), upon remembering that for a fast moving lepton with helicity A one has ~ 2\k^ /m. 
For positively charged leptons (e"*", /x"*") one should simply replace, in Eqs. (2.2.3, 4), qa 
by —Qa- In our notation 

gy = -\ + 2sin^ew 9a = -\- (2.2.5) 

The factors rf in Eq. (2.2.1) collect some kinematical factors, coupling constants 
and the relative weights of different propagators, namely: 

where G is the Fermi coupling constant and Mz is the Zq mass. Notice that GM|/(2V^7rQ;) 
= (4sin2^VKCos2^n.)-i ~ 4/3. 

Finally the hadronic tensor W^"^ defines the coupling of the electro-weak current 
to the nucleon; exploiting Lorentz and CP invariance it can be expressed in terms of 8 
independent structure functions as [BAR 79] 



El + 2?^^ S^^^ - ^xP'^S^ gi 
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(2.2.7) 



P/j-Su + 'S'^-Py i S ■ q j S ■ q ^ 

2P-q ^3+ (p.^)2 



Notice that terms proportional to or q'^ can be dropped in the definition of 

WJj^^ because they give no contribution (in the m/E — > hmit) when contracted with 
L^'^ . Other definitions of the hadronic tensor appearing in the hterature differ from ours 
because of these terms. In particular the coefficient of g2 could be written in a more 
familiar way using the identity [HEl 73] 

W^"^^ = :^^[iQ-S)q''P^-(P-q)q''S^] 
IxF ■ q 

— {Qfi^ua/Sj — Qu£nal3-y) ^ (2.2.8) 

and dropping the last two terms. 

The structure functions Fj{P ■ q, q^) and g'^ {P ■ q, q^) are expected to scale approx- 
imately in the Bjorken hmit and to depend only on x = Q'^/{2P ■ q). The are the 
unpolarized structure functions and the the polarized ones: when averaging over the 
nucleon spin one sums lV^i/(P, q, S) to W^i^{P, q, —S) and all terms proportional to gj can- 
cel out. In Eq. (2.2.7) we have allowed for parity violation and indeed IV^^, is a mixture 
of second rank tensors and pscudotcnsors. In case of pure electromagnetic interactions 
(i = 7) parity is conserved and one has 

= 9l = gJ ^gl^o- (2.2.9) 

In this case Eq. (2.2.7), upon using Eq. (2.2.8) and up to irrelevant terms propor- 
tional to q^ or q'', reproduces Eqs. (2.1.8-10, 16, 18). F3, g^, gi and gr, only contribute 
to parity violating interactions and are often referred to as the parity violating structure 
functions. 

From Eqs. (2.2.1-4) and (2.2.6, 7) one obtains explicit expressions of the cross- 
sections; some of them can be found in [ANS 93]. 

In case of charged current interactions, t^N v{i>)X, the process is, at leading 
order, mediated by the exchange of a meson and it resembles the Zq contribution to 
the neutral current process, with the assignments gy — g^ — 1. In fact Eq. (2.2.1) now 
reads ^ ^ ^ 



with given by Eq. (2.2.7) in which the label i is now W. For a negatively charged 
lepton £~ (which couples to a W~), 

LZ~ = E i^ik', s') 7^(1 - 75) u{k, X)Y [u{k', s') 7.(1 - 75) u{k, A)] 

s' 

^ {l-2\f LI, . (2.2.11) 

Equation (2.2.11) shows that fast ^~ leptons only couple to a IV if they have a negative 
helicity (A = —1/2). Equation (2.2.10) is completed by 

UGM^w \ .2 2 12) 

where is the W mass. For a positively charged lepton one simply changes, as 
usual, the sign of the axial coupling 7^7^, that is one replaces A with —A in Eq. (2.2.11) 
to obtain Lj^^ = (1 + 2A)2 Ll^. 



2.2.1 Weak interaction structure functions and their measurement 

In order to obtain information on the polarized structure functions observed in 
weak interaction deep inelastic scattering we should study, in analogy to Eq. (2.1.22), 
differences of cross-sections with opposite nucleon spin. Rather than considering here the 
most general case, which could be derived from Eqs. (2.2.1-7) or (2.2.10-12) and (2.1.21), 
we look at particular nucleon spin configurations, namely longitudinal or transverse ones, 
like in Eqs. (2.1.24) and (2.1.25) respectively. 

To do so we switch to the more convenient variables x and y = i^/E, using 



E' 



d^ 
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dn dE' yME dx dy dip 



(2.2.13) 



We then notice that in the high energy region one can neglect terms proportional to M/ E 
and that, from Eq. (2.2.6) and for Q'^ values up to ~ 10^ GeV^, one has 



(2.2.14) 



Taking all this into account and adopting the same symbols to denote the nucleon 
spin as in Eqs. (2.1.24) and (2.1.25), one obtains for neutral current reactions: 
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and for charged current reactions initiated by £^ 
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(2.2.18) 
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Equations (2.2.15, 16) show that in our kinematical range the neutral current pro- 
cesses are still dominated by one photon exchange and thus agree with Eqs. (2.1.24, 25) 
respectively. The charged current cross-sections instead, although much smaller than the 
neutral current ones due to the factor rj^ , Eq. (2.2.12), depend also on the parity violating 
polarized structure functions g^^^. (Recall that the ± signs and the labels in Eqs. 
(2.2.17, 18) refer to negatively or positively charged leptons, £^.) 

One could also extract information on the parity violating structure functions in 
neutral current processes by looking at appropriate combinations of cross-sections [ANS 
93] . Ideally, it would be very helpful to perform large experiments with both £~ and 



i'^ leptons with positive (— >) and negative (■<— ) helicities. In fact, always within the Q'^ 
range of validity of Eq. (2.2.14) and noticing, from Eq. (2.2.5), that Qy ~ 0.04 whereas 
gi^ = —0.5, one has: 
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where £~ and ^^ label negative and positive charge leptons respectively. Other interesting 
combinations of cross-sections can be found in [ANS 93]. Notice that in Eqs. (2.2.19 
and 20) a sum over the £~ and helicities is performed, which amounts to considering 
unpolarized leptons. 

In case only i~ unpolarized lepton beams are available one can still obtain infor- 
mation on the parity violating polarized structure functions, by scattering the leptons off 
longitudinally polarized nucleons and measuring [BIL 75]: 

ifd'^a^ d'^at^ \ l/d^cr^ d^at^ \_ d'^a^ d'^a^ 
2\ dx dy dx dy J 2\ dx dy dx dy J dx dy dx dy 



IGttME^ {{1-y) [g,r,^'{gf - gf) + g\rf{gl - gf)] (2.2.22) 



+ [gyv'' gf + g'y gi] + xy{2 - y)g,v''gf} ■ 

Eq. (2.2.22) holds also for unpolarized leptons, provided g^ is replaced by —g^, so 
that on summing and subtracting the £~ and the contributions one recovers respectively 
Eqs. (2.2.19) and (2.2.20). Of course, the above measurements (2.2.19-22) are difficult 
ones in that they single out non leading parts of the cross-sections. However, they may 
be feasable and would provide further information on the structure of the nucleon. 



3 The Naive Parton Model in polarized DIS 

3.1 Projection of Fi,F2,gi and g2 from the hadronic tensor W^^ 

In Section 2 we have written the most general hadronic tensor W^i,, Eqs. (2.1.8-10), 
which describes the unknown coupling of the virtual photon to the composite nucleon 
in terms of the four structure functions Fi = MWi, F2 = gi = A'PuGi and 

g2 = Mu'^G2] in the large Bjorken limit, Eqs. (2.1.16, 18), these arc supposed to 
scale, that is to depend only on the Bjorken variable x = Q"^ /2Mv. By performing DIS 
experiments one learns about the structure functions. 

The hadronic tensor W^^, can also be computed from models of the nucleon, as we 
shall soon illustrate using the Parton Model. In such a case it is convenient to have a 
technique which easily allows one to extract the four different structure functions from a 
knowledge of W^^{N). 

Let us start from the unpolarized case. Defining 
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one can see, from Eqs. (3.1.1-3) and (2.1.8-10), that 



and 



Similarly, in the polarized case one can introduce the projectors 

[{q ■ S)S^ + qx] P.e''^''' 
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such that, from Eqs. (3.1.6-8) and (2.1.8-10), 

PfW^^iN) = g2 (3.1.9) 
PfW^fsiN) = 51 + ^2. (3.1.10) 

Equations (3.1.1-10) can be used in any reference frame and for arbitrary directions 
of the nucleon spin four-vector; however, some care must be taken in particular cases in 
which a hmiting procedure is required. 



3.2 The hadronic tensor and the nucleon structure functions in the Naive 
Parton Model 

In the simplest version of the Parton Model the nucleon is considered to be made 
of coUinear, free constituents, each carrying a fraction x' of the nucleon four-momentum. 
The lepton-nuclcon DIS is then described as the incoherent sum of all lepton-constitucnt 
quark interactions and the hadronic tensor PV^,^(iV), Eq. (2.1.8), is given in terms of the 
elementary quark tensor W/^i, by [LEA 85] 

W,,{q;P,S) = W!£\q;P)+tWl^f{q;P,S) 

= — / — S{x' - x) ng{x',s;S) w^^{x',q,s) , (3.2.1) 

q,s ■ q JO X 

where ng{x', s; S) is the number density of quarks q with charge Cg, four- momentum frac- 
tion x' and spin s inside a nucleon with spin 5" and four-momentum P; the Eg runs 
over quarks and antiquarks; x is the Bjorken variable (2.1.15) and the quark tensor 
Wfj,u{x,q,s) is the same as the leptonic tensor L^^^, Eqs. (2.1.2-7), with the replacements 
f^fj. _^ xP^, k'^ xP^+q^ and a sum over the unobserved final quark spin (s') performed. 
That is: 

w^y{x, q, s) = wlfj {x, q) + tw^^f {x, q, s) (3.2.2) 

with 

wll\x,q) = 2[2x^P^P, + xP^q, + xq^P,-x{P-q)g^''] (3.2.3) 
wl^\x,q,s) = -2m, £^,„^sV (3-2.4) 

and the quark mass must for consistency be taken to be nig — xM, before and after the 
interaction with the virtual photon. Wc will further comment on this point in Section 3.4. 

From Eqs. (3.1.1-5) and (3.2.1-4) one obtains the well known Naive Parton Model 
predictions for the unpolarized nucleon structure functions: 

^i(^) = lT.4<li^) (3-2.5) 

^ q 

F2(x) = xY,elq(x)^2xF,(x), (3.2.6) 
q 

where the unpolarized quark number densities q{x) are defined as 



?(^) = ^ng{x,s]S) 

s 



(3.2.7) 



Similarly, from Eqs. (3.1.6-10) and (3.2.1-4), the polarized nucleon structure func- 
tions are obtained: 



^i(^) = ^Ee?Ag(a:,5) (3.2.8) 

^ q 

g2{x) = (3.2.9) 

where 

Ag(a;, S) = nq{x, S; S) - ng{x, -S; S) (3.2.10) 

is the difference between the number density of quarks with spin parallel to the nucleon 
spin {s = S) and those with spin anti-parallel (s = — -S"). In the Parton Model with free 
non-interacting quarks of mass niq — xM, this quantity cannot depend on the choice of 
S, i.e. Aq{x,S) = Aq{x); this follows because in such ultra simple case one is allowed 
to Lorentz transform to the nucleon rest frame, where the quarks too are at rest and 
where the parton distributions nq{x, ±S, S), with S = (0, S), cannot depend on the spin 
quantization direction S [ANS 92]. 

3.3 Intrinsic p± and g2{x) 

The polarized structure function g2{x), Eq. (3.2.9), is zero in the Naive Parton 
Model where each parton carries a fraction x of the nucleon four-momentum, — xP^. 
However, non zero values of g2 can be obtained by allowing the quarks to have an intrinsic 
Fermi motion inside the nucleon, p'^ — {Eq,px,Py,x'Pz) (in a nucleon moving along the 
5-axis). 

Equation (3.2.1) then modifies to [LEA 85] 

W,,{q; P,S) = J2elJ d'p (^^ j S{2p ■ q - Q') nq{p, s; S) w,,{q;p, s) (3.3.1) 

which holds in the infinite momentum frame {\Pz\ — oo) and where nq{p,s;S) is the 
number density of quarks q with charge Cq, energy Eq, three-momentum p and spin s in 
a nucleon with spin S and momentum P. It is related to the usual quark distribution 
functions q{x) and Aq{x, s; S), Eqs. (3.2.7, 10), by: 

Y: I d^P± riqip, s; S) = g(x') (3.3.2) 

J d^p^ [nq{p, s; S) riqip, -s; S)]Pz = J d^p^ Aq{p, s; S) P, = Aq{x' , s; S) (3.3.3) 

where we have put d'^p = d^pj_ dx'Pz . 

Equations (3.1.9) and (3.3.1) yield 

92{x) =E4l d'p (§) 5{2p .q-Q') nq{p, s; S) Pfwi^^{q-p, s) , (3.3.4) 



q,s 



and insertion of the explicit expression of P^^ , Eqs. (3.1.6, 8) and w^^'*, Eq. (3.2.4), yields 

- E^^I^/'i^pgj^Pp-.-e^) (3.3.5) 
X A5(p, s; S){(P ■ q)[(q ■ S){s ■ S) + q- s\- {P ■ s)[(q ■ Sf + 9^1} 



with b as given in Eq. (3.1.8) and the quark mass m,q = ^E"^ — p^. 

In the nucleon rest frame let S'^ = (0, S), where S is not perpendicular to the z-eods 
(along which the nucleon moves in the infinite momentum frame) ; specify the quark spin 
by choosing s'^ = (0, S) in the quark rest frame (which, for ^ 0, is not the nucleon 
rest frame). Then Eq. (3.3.5), at leading order in M/P^, reads [LEA 88] 

«^W = 5i:'?(l-— JS^A,(...;S) (3.3.6) 

where we have used Eq. (3.3.3). 

As the left hand side cannot depend on s and S, Eq. (3.3.6) appears contradictory. 
However, this is not so since in a fast frame where p± is negligible compared to Pz, both 
and S^, to leading order, are parallel to so that Aq is the hehcity distribution and 
we have 

92{x) = 11:4{^-i)M{x) (3.3.7) 

where Ag(x) is the difference between the number density of quarks with the same helicity 
as the nucleon and those with opposite helicity. When neglecting the intrinsic p± of quarks 
one has ruq = xM, so that from Eqs. (3.3.6) and (3.3.7) one recovers the previous result 
g2{x) = 0. Thus the introduction of intrinsic p± allows a non zero value of g2{x) in the 
Parton Model, but we see that such value shows an extreme sensitivity to the quark mass, 
Eq. (3.3.7). 

We shall argue, as has been emphasized in [JAF 90] , that this throws doubt on all 
purely Parton Model calculations of g2. 

3.3.1 Conflicting Parton Model results for g^^x) 

Before discussing the difficulties of the Parton Model with g2{x), and in order to 
better realize that indeed there arc difficulties, we briefly summarize some contradictory 
statements existing in the literature. 

Some authors [lOF 84] claim that 

9i{x) +g2{x) = , 

whereas Feynman [FEY 72] has 

gi{x)+g2{x) , (3-3.9) 

^ q 

where Aq^ is given by Eq. (3.2.10) with the spin S perpendicular to the nucleon mo- 
mentum. Moreover, on summing Eqs. (3.2.8) and (3.3.7), one obtains, in terms of helicity 
densities [LEA 88], 

1 TTi 

gi{x)+g2{x) = - ^?(^) ■ (3-3.10) 

Finally in [JAC 89] it is argued, using a covariant Parton Model formulation, that 

gi{x)+g2{x)= C—gi{y). (3.3.11) 
Jx y 

These last authors exploit the experimental data on gi{x) to give predictions for 
xg2{x) [see Fig. 3.1]. 



3.3.8) 



In [ANS 92] it is shown that Eq. (3.3.8) is not correct and that Eqs. (3.3.9) and 
(3.3.10) are compatible when neglecting the intrinsic of quarks. In section 5.3, we 
explain why (3.3.11) is incorrect because it neglects twist-3 terms. Nevertheless (3.3.11) 
may not be a bad approximation, so it should be tested experimentally. 

Finally in Fig. 3.2 we show a bag model prediction for g2{x) [JAF 91]. 

3.4 Origin of the difficulties with g^ix) 

As we have seen the Parton Model does not lead to clearcut results for the spin 
dependent structure function g^. Different authors obtain different results and these are 
generally incompatible with each other. 

Perhaps the most general way to understand this difficulty is to recall that the 
Parton Model is fundamentally an impulse approximation in which binding effects (i.e. 
the virtuality) of the struck parton are unimportant for the large transverse momentum 
reactions under consideration. However, in some cases what one is measuring is not just 
a cross-section but an asymmetry, i.e. a difference of cross-sections, and it may happen 
that in this difference the dominant contributions cancel out leaving a result which does 
depend upon the binding energy or virtuality. Just such effects occur in g2. Whenever this 
happens the result is bound to be unreliable. 

We can see this quite directly by reconsidering the antisymmetric part of Eq. (3.3.1), 
valid in the impulse approximation: 

Wl^^\q; P,S)^Y.4j d'P (^) '^(2P ■ Q - Q') n,(p, s; S) w^^{q-p, s) (3.4.1) 

where we ignore the sum over the flavours irrelevant to this issue. 

Now consider the calculation of w'^^'> describing the interaction of the hard photon 
with the quark, as depicted in Fig. 3.3. The final state quark is a 'free' quark and is on 
mass-hell: {p'Y = rrig. In the impulse approximation also the initial quark is considered 
to be free and to have the same mass. But to see the danger of this assertion, let us put 
p2 — m^, where for the moment we allow ^ to represent the fact that the initial 
quark is really a bound quark. Aside from this we treat the incoming quark as free, i.e. 
its wave function is taken as the usual free-particle Dirac spinor u{p, s) for a particle of 
mass m. 

One then finds 

w^,u{q;p, s)^^Tr [(1 + -f5^)(J + m) 7^ (|5 + ^ + m,) 7,] (3.4.2) 
from which one obtains 

■wlt\q;p,s) =2e^^ap{mqs'') 

Equation (3.4.3) is extremely reveahng. We see immediately that for a general 
the result is not gauge invariant {q^wj^^ 7^ 0) unless m = nig, in which case we recover 
the Naive Parton Model, Eq. (3.2.4). Moreover, the offending term, when m 7^ rriq, is not 
small in an infinite momentum frame (where the impulse approximation is supposed to 
be most justifiable) even if (m — niq) is small. 

However, in the special case of longitudinal (L) polarization, if the quark has high 
momentum so that ruq/pz <S 1, the product {mqS^i) — > and the gauge non- invariant 



l-^]p^-^q^ 



rrin 



rrin 



(3.4.3) 



term vanishes because of the antisymmetric Sptvap- For the calculation of gi{x) it is the lon- 
gitudinal polarization that is relevant and there is no unreasonable sensitivity to whether 
or not m equals ruq. Moreover, since the limit limm,^o {^qSi) is finite there is also no 
crucial sensitivity to whether one works with massive or massless quarks. Thus gi [x] can 
be calculated unambigously in the Parton Model. 

Quite the opposite happens for g2{x) where the transverse spin is relevant. There 
is an extreme sensitivity to whether or not m equals rUq and one cannot expect to make 
a reliable calculation of g2{x) in the Parton Model. 

One can put the case even more forcefully. The whole point of quarks is that, in 
their pointlikc interaction with a hard photon, they produce large momentum transfer 
reactions which we are trying to generate for the photon-hadron interactions. But even if 
we define the model by insisting that m — uig, comparing the expression (3.4.3) with the 
general structure of VFj^^ for a spin 1/2 particle (2.1.19), we see that 

g2{x)\quark = . (3.4.4) 

Thus the hard photon-free quark interaction does not possess the cross-section 
asymmetry which we are seeking to explain in the photon-hadron interaction. It is clearly 
unrealistic therefore to try to produce such an asymmetry from quark-partons. 

3.5 WeEik interaction structure functions in the Naive Pcirton Model 

We have seen in Section 2.2 that, when taking into account weak interactions, the 
most general form of the hadronic tensor, Eq. (2.2.7), contains 8 independent structure 
functions, 3 unpolarized and 5 polarized ones. On the other hand, in the Naive Parton 
Model, the hadronic tensor is given by (see Eq. (3.2.1)) 

W;,(q; P, 5) = E 2{p\)x ^' + ^' '^^^ ^^'^■'^^ 

where we have cxplicitely performed the sum over s = ±S. The index i denotes the differ- 
ent kinds of contributions {i = 7, jZ, Z, W) and the J2q runs over quarks and antiquarks. 
The are the (fiavour dependent) quark tensors, given by: 



wl'^ix, q,s) = Yl 4 l.'^iP'^ s'hMP^ Hp'^ s')-f^u(p, s)] 

s' 

wl^'^x, q,s) = Mp'^ «')7/.(^y - 9A'y5)qu{p, s)]* [u{p' , s')'y^u{p, s)] 

s' 

+ Jl^q ["^(P'^ s'hMP^ ["^(P'' s'hi^iQv - 9a15)u(p, s)] (3.5.2) 
w^uix, = Yl ["^(P'^ ^'hi^i9v - 9Al5)qu{p, s)]* [u{p', s')-f^{g^ - g^-f5)qu{p, s)] 

s' 

Q,s) = Y Hp'^ ^')7m(1 - 75)«(P, s)r Hp', ^')7.(1 - 75)«(P, ^)]|(V^),,f 

s',q' 

where p — xP,p' — p + q, s and s' are respectively the momentum and spin four- vectors 
of the initial and final quarks. {gy)q and (g'^)g are the vector and axial couplings of the 
quark of flavour q to the Zq: 

{9v)n,c = ^-^sin^Ow {9a)u,c = ^ 

12 1 
{9v)d,s = -- + -sm'^9w {9A)d,s^-^- (3.5.3) 



In case of charged current, negatively charged leptons couple to u-type quarks (or 
d-tjpe antiquarks) and positively charged leptons couple to c?-type quarks (or -u-type 
antiquarks); one has also to take into account the proper Cabibbo-Kobayashi-Maskawa 
matrix elements {V)qq' occuring in the transition coupling from a flavour g to a flavour 
q'. However, if we consider the contribution of four flavours {u, d, s and c), one always has 
J2q, \ {V)qq'\'^ — cos^ 9c + slu^ 9c — 1, whcrc 9c is the Cabibbo angle. Equations (3.5.2) and 
(3.5.3) hold for quarks; for antiquarks one should only replace 75 with —75 in Eqs. (3.5.2) 
leaving all other expressions unchanged. 

Explicit expressions for Eqs. (3.5.2) can be obtained from the general form 

w^^{x, q,s) = Yl ['^(P'^ s')7/.(vi - ai^5)u{p, s)]* [u(p', s')^^(v2 - a2-f5)u(p, s)] 

s' 

= 2{aia2 + viV2)[2p^Pu-p-q gf,u]-4:aia2mlg^i, (3.5.4) 

- 2via2 mq[2p^Su - s-q g^^] - 2aiV2 mq[2s^pu - s-q g^.^] 

+ 2ie ^ycci3[{via2 + 01^2 )p"g'^ + 2aia2 m.qp'^s'^ + (0102 + ViV2)mq g"s'^] 

by properly fixing the values of Vi^2 and 01^2! in Eq- (3.5.4) we have dropped, as usual, 
terms proportional to q^ or which, when contracted with the leptonic tensor L^^^ give 
negligible contributions proportional to m/ E. 

Equations (3.5.1-4) give the Parton Model prediction of the hadronic tensor PV^,^; 
by comparing it with the general expression (2.2.7) one obtains the naive quark-parton 
model results for the nucleon structure functions. For completeness we list all of them 
here, starting from the electromagnetic case (i = 7), for which we recover Eqs. (3.2.5, 6): 

^1" = ^E<(5/ + 57) f:i = 2xF^, 
2 / 

9l = ^E<(^?/ + ^57) ^1 = (3.5.5) 

where g/ = m, s, c; u stands for the number density of quarks u and so on. The interfer- 
ence contribution {i — 7Z) is: 

f 

i 

gf = Ee.,(5j,,(Ag/ + Ag^) (3.5.6) 
/ 

92 = 94 =0 

9f = 2x^eg^(5^)g^(Agj- Ag/) =2a;^J^ 
/ 

and the purely weak interaction {i — Z) leads to: 

Ff = lT.i9l + 9%A^f + Qf) Fi^2xFf 
^ f 

Fi = ^Ei9v9Mqf - Qf) 
f 

9l = W{gl^g\)q,{^q^^^qi) (3.5.7) 



9i = -lj:i9%fi^<if + ^Qf) 

f 

gl = 0. 

In case of charged current (z = W), on performing explicitely the J2fj one obtains, for 
i~N vX processes: 

" = u + c + d + s F^' = 2xF^~ 
F^~ = 2{u + c-d-s) 

gY~ = (Am + Ac + AJ+ As) = -2^f " (3.5.8) 
gY~ = 2x{Au + Ac - Ad - As) = 22:55^" 

gr = 0. 

i'^N vX processes probe different quark flavours and one obtains the corresponding 
expressions of the structure functions F^"^ and (?]^^ by the flavour interchanges d ^ u 
and s c in the above Eqs. (3.5.8). 

Notice that in the naive quark-parton model the structure functions g\ are always 
zero and one finds g\ = 2xg\ for any i = 7, 7^', Z, W . The functions g^ are nonzero only for 
pure weak interactions (both in neutral and charged current processes). It is interesting to 
note that, in case of neutral current, the integral of is directly proportional to the total 
spin carried by the quarks and antiquarks, as can be seen from Eqs. (3.5.7) and (3.5.3). 
Particular combinations of the structure functions can single out interesting quark and 
quark spin information; some examples are discussed in [ANS 93]. 



4 Phenomenological analysis of the data on g\{x) and its first 
moment Fi 

As explained in Section 2.1.3 the measurement of the cross-section asymmetry using 
a longitudinally polarized lepton beam on a longitudinally polarized hadron target may 
be interpreted as essentially a measurement of the spin dependent structure function 
gi{x^Q'^). Experimental data on g^ were first obtained at SLAC [ALG 78] as early as 1978 
using an electron beam. Further information came from the SLAC- Yale group [BAU 83] in 
1983, with fascinating implications about the internal structure of the proton. And then, 
more recently, very startling results were obtained by the European Muon Collaboration 
(EMC) who scattered a longitudinally polarized muon beam of energy 100-200 GeV on 
a longitudinally polarized hydrogen target at CERN [ASH 88, 89]. The unexpectedly low 
asymmetry found by the EMC led to what was termed a "spin crisis in the Parton Model" 
[LEA 88] and raised serious questions as to how the spin of the proton is built up from 
the spins of its constituents. It also seemed, for a while, to imply the first failure of the 
Parton Model Q- 

The EMC result catalysed a great deal of theoretical research and the "crisis" is 
now beheved to be resolved as a consequence of the discovery of a deep and beautiful 

It is often forgotten that there are in a sense other failures of the NPM, namely the experimental 
observations of large transverse polarization in inclusive A-production at high px first observed in 
1976 and of left-right asymmetry in high pT inclusive pion production on a transversely polarized 
target (see e.g. [HEL 91]). 



connection between the observable measured by the EMC and the axial anomaly [EFR 
88] . Nonetheless the situation is not yet fully settled and there are subtle non-perturbative 
issues and ambiguities in limiting procedures which will be discussed in Section 6. 

Two new experiments are under way to check the EMC measurement of Qi and to 
measure the analogous neutron function g'^. First results have now been presented [ADE 
93, 94; ANT 93] and we shall discuss their implications in Section 4.5. Our discussion in 
this Section will largely rely on the EMC experiment, since this is the experiment at largest 
(Q^) and should thus be least sensitive to higher twist contributions. It is important to 
realize that even if the EMC results turn out to be incorrect, the theoretical discoveries 
catalysed by them remain valid and interesting. In fact at the time of going to press, the 
SMC has just issued results at the same (Q^) as the EMC experiment [ADA 94] and they 
are in substantial agreement with those of EMC. We will comment briefly on the new 
measurement in Section 4.2.2. 

We shall begin by analysing the data in the simple Parton Model. We shall find 
ourselves in difficulty and will therefore invoke the more sophisticated operator product 
expansion with its QCD corrections. This will be seen to offer no significant help and we 
will be forced to conclude that the standard approach is in trouble. A possible resolution 
of these difficulties in terms of an anomalous gluonic contribution will be discussed in 
Section 7, where we deal with the first moment of gi{x) and in Section 8.4 where we 
consider the x-dependence. 

4.1 The SL AC- Yale and EMC data: quark distributions near x = 1 

As derived in Section 3.2, in the Naive Parton Model one has 

g{{x) = [Au{x) + Au{x)] + ^[Ad{x) + Ad{x)] + ^[As{x) + As{x)]^ (4.1.1) 

where we assume that the contribution of heavy quarks is negligible. 

Before discussing the data on gi{x) itself it will be interesting to look at the virtual 
Compton scattering asymmetry Ai{x) introduced in Appendix A. Taking the approximate 
form 

M.) ^ HI (4.1.2) 

and using the standard result 

= [u{x) + u{x)] + ^[d{x) + d{x)] + ^[s{x)+s{x)]^ (4.1.3) 

one has, using (4.1.1), 

^ 4[Am(x) + Au{x)] + Ad{x) + Ad{x) + As{x) + As{x) ^ 

^^^^ ^ 4:[u{x) + u{x)] + d{x) + d{x) + s{x) + s{x) y--) 

The EMC and early SLAC-Yale data [ALG 78; BAU 83] on Ai{x) are shown in Fig. 
4.1. The errors are, of course, large, but there is a discernible trend towards the maximum 
possible value Ai = 1 as x — > 1. 

Given that we know empirically that u{x) ^ d{x), s{x), q{x) as a; — 1 and that 
one must have |Ag(a;)| < q{x) for each quark flavour, the behaviour Ai —>■ 1 as x ^ 1 
suggests that 

Au{x)\^^i ^ u(x) . (4.1.5) 

This implies that those u quarks which carry a large fraction of the longitudinally 
polarized proton momentum are highly polarized along the direction of the proton spin. 



4.2 Analysis of data in the framework of the Operator Product Expansion 

As will be discussed in Section 5, the operator product expansion gives results for 
the moments of gi^2 in terms of hadronic matrix elements of certain operators multiplied 
by perturbativcly calculable coefficient functions. The general result for the moments of 
gi is given in (5.2.1). If, however, the fields involved in the electromagnetic currents in 
(5.1.1) are treated as free quark fields the results simplify via (5.1.10) - in effect the current 
commutator in (5.1.1) can be explicitely evalutated - and one obtains 



r? = £ dx g^^{x) = ]^{«3 + -4 as + ^ aol , (4.2.1) 



where the are [aside from a factor specified below in (4.2.7 and 8)] hadronic matrix 
elements of the octet of quark SU{3)f axial-vector currents J^^ (j = 1,...,8) and the 
flavour singlet axial current J^^ taken between proton states of definite momentum and 
spin direction. (In practice we usually use helicity states). The precise definitions are given 
below. 

In the following we shall explain how the values of 03 and og can be obtained from 
data on hyperon /3-decay. The EMC measurement of can thus, via (4.2.1) be construed 
as the first ever measurement of the singlet matrix element Oq and it will turn out that 
the value of Gq is unexpectedly small. 

We shall show that in the Naive Parton Model 

ao = 2SI (4.2.2) 

where is the component of the total quark spin (carried by all quarks and antiquarks) 
in the direction of motion of a proton of helicity +1/2. That turns out to be almost 
compatible with zero is a great surprise. Naively it might have been expected to be fairly 
close to +1/2 implying ao — 1. When the quark fields are treated as interacting fields 
there are small perturbative corrections to (4.2.1) (see Section 4.4.1) but these have no 
significant effect on the argument below. 
The octet currents are 

Ji, = ^7^75^^)^' (j = 1,2,. ..,8) (4.2.3) 

where the Xj are the usual Gell-Mann matrices and V' is a column vector in flavour space 

(4.2.4) 




and the flavour singlet current is 

= -^7^75^^ ■ (4.2.5) 

(Note that Jg^ is sometimes defined like (4.2.3) with a fiavour matrix Ao/2 = (l/\/6) /.) 
The question of the conservation of the non-singlet axial currents will be discussed in 
Section 6.1. 

The proton states will be labelled by the four-momentum and the covariant spin 
vector S^{X) corresponding to definite helicity A. Recall that 



(4.2.6) 



The forward matrix elements of the J^^ can only be proportional to S'^ and the Uj are 
conventionally defined by 

{P,S\Ji^\P,S) = MajS^ (4.2.7) 
{P,S\Jl\P,S) = 2MaoS^. (4.2.8) 

The relative factor of 2 in (4.2.7 and 8) reflects the fact that the SU{3) currents are 
defined using the generators of the group i.e. \j/2 in (4.2.3). 

Analogous to (4.2.3) one introduces an octet of vector currents 

^^ = ^7^(y)V' (j = 1,-,8) (4.2.9) 

which are conserved currents to the extent that SU{3)f is a symmetry of the strong 
interactions. 

4.2.1 Information from hyperon (3-decay 

Consider now the /9-decays of the spin 1/2 hyperons. If we use the standard S'f/(3) 
labelling [BAI 82] Bj for the hyperons then the hadronic transition involved is controlled 
by matrix elements of the form {Bi\h''WBi^) where U\. is the charged hadronic current that 
couples to the W boson in the electroweak Lagrangian. It is typically of the form of a 
Kobayashi-Maskawa matrix element multiplied by some combination of the currents 
and Jg^. For example AQ = — 1 transitions (A — > p; E ^ n; ^ A) are controlled 
by the current 

u-i^{l - 75)5 = [Jl + iJl - Jt^ - ^4] . (4.2.10) 

If we now assume: 

a) that the 8 spin 1/2 hyperons form an octet under SU{2>)p] 

b) that the currents J^, Jg^ [j = 1,...,8) transform as an octet under SU{3)f with 
J^, J^^ conserved, (the conservation of the axial current will be discussed in Section 
6.1); 

c) that the momentum transfer and mass differences in the hadronic transitions are 
negligible; 

then all the hyperon /3-decays arc described in terms of two constants, F and D which 
occur in the matrix elements of the octet of axial currents [BAI 82]: 

{Bj;P,S\Jl^\Bk;P,S) = 2MBS^{-zf,,kF + d.jkD} ii,j,k = 1,...,8) (4.2.11) 

where the fijk and dijk are the usual SU{3)f group constants. The structure of (4.2.11) 
is dictated by the SU (3) transformation properties of the LHS and is an example of the 
Wigner-Eckart theorem. It should be noted that J^^ does not play any role in the weak 
interactions. 

An analysis in 1983 of hyperon decay data [WA2 83] seemed to be in good agreement 
with the above and led to the values 

0.477 ±0.012, D = 0.756 ±0.011. (4.2.12) 

If now we use the standard SU (3) f assignments for the baryon octect [BAI 82] we 
find from (4.2.11) that 

as = F + D (4.2.13) 
08 = ^ (3F - D) . (4.2.14) 



We see, therefore, that the study of hyperon /9-decay provides us, via (4.2.12), with 
the values of and ag, two of the three parameters occuring in (4.2.1). However there 
is some doubt as to the rehabihty of the values quoted in (4.2.12), and for that reason 
we shall adopt a slightly different strategy. Ultimately, however, it will be seen that the 
essential result is not sensitive to this issue. 

It is straightforward to demonstrate, using only isotopic spin invariance, that 

as ^9^, (4.2.15) 

where is the axial coupling (= Ga/Gv in the language of the old Cabibbo theory) which 
governs neutron /3-decay. We prefer therefore to use the very accurate measurement of 
[PDG 92] 

as^g^^ 1.2573 ± 0.0028 (4.2.16) 
rather than the value F + D^ 1.233 ± 0.016 that follows from (4.2.12). 

4-2.2 The EMC data on and its consequences 

Our main interest will be in the matrix element of the flavour singlet axial current, 
i.e. in qq. Prom (4.2.1) we have 

ao = ^|l2r?- as --^agj . (4.2.17) 

The EMC data on xgi{x) are shown in Fig. 4.2. Of particular interest, and indeed 
the source of all the recent theoretical research in this subject, is the first moment 

r?= l^dxglix). (4.2.18) 
Jo 

Since experimental measurements cannot reach x = or 1 an extrapolation of the 
data is always required in evaluating (4.2.18). The region a; — > 1 is harmless but the 
extrapolation to a; = is potentially tricky. The dashed curve in Fig. 4.2 shows the form 
of the extrapolation used by the EMC for the value of Jl dx' g\{x') as function of x 
for a: —> 0. It is seen that the extrapolation looks perfectly reasonable (but beware the 
logorithmic scale) and it led to the vahic = 0.126 ± 0.010 ± 0.015. 

Recall [see (2.1.31)] that to extract gi{x) from the measurement of the asymme- 
try A\\ one has to utilize information on F2{x,Q'^) and R{x,Q^) which are studied in 
unpolarized DIS. Since the EMC experiment was carried out there have been improved 
determinations of F2 by the NMC group [AMA 92] and of R at SLAC [WHI 90]. As a 
consequence there is a small change and one has [ELL 93] 

[(g2) = 10.7] = 0.128 ± 0.013 ± 0.019 (4.2.19) 

where we have indicated the mean value of for the experiment. 

In Fig. 4.3 we show the new SMC data [ADA 94] on g^{x) at (Q^) = 10 (GeV/c)^ 
compared with the EMC data. There is excellent agreement in general where the measure- 
ments overlap. The new data, which extend to smaller value of x give some indication that 
gi{x) may be increasing as x — 0, though this trend may not be significant. The SMC 
value for the first moment is F? [(Q^) = 10] = 0.136 ± 0.011 ± 0.011, perfectly consistent 
with the EMC result (4.2.19). 



This result is considerably smaller than the value 0.188 ± 0.004 expected on the 
basis of the Ellis- Jaffe sum rule which will be discussed in Section 4.3. Several attempts 
were made to try to explain this discrepancy, initially focusing on the question as to 
whether (4.2.19) can be considered as a reliable number to compare with the theory. If 
we substitute the value of (4.2.19) into (4.2.17) we obtain 

ao = ^|(0.279 ± 0.156 ± 0.228) - ^(3F - D)| . (4.2.20) 

The values (4.2.12) yield (3F - £>)/3 = 0.225 ± 0.013 which implies oq = 0.040 ± 0.117 ± 
0.171. However a more recent analysis of hyperon decays obtains [HSU 88] 

F = 0.46 ±0.01 D = 0.79 ±0.01 -(3F - D) = 0.20 ± 0.01 (4.2.21) 
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which implies 

ao = 0.06 ±0.12 ±0.17. (4.2.22) 

As mentioned earlier one might have expected oq — 1 so that the value in (4.2.22) is 
surprisingly small. 

This unintuitive result and the disagreement with the Ellis- Jaffe sum rule inspired 
several attempts to question the reliability of the EMC value for F^. Because of its fun- 
damental importance we turn to this question. We believe that only two issues require 
serious consideration: the question of the extrapolation to a; = and the question of higher 
twist effects. 



4-2.3 The EMC data: extrapolation to x = 

The extrapolation to small x was queried by Close and Roberts [CLO 88] who 
stressed that the extraction of gi{x) from the data, at small x, is dependent upon the 
assumed behavior as a; of both gi{x) and the usual structure function ^2(0;), and 
that, in particular, the latter may have a more singular behavior than is given by the 
usual Regge analysis of small x behavior. 

We are convinced that the Regge behavior 

gi{x) ~ x-'^-^ (4.2.23) 

where aai(t) is the Regge trajectory of the oi meson (previously referred to as the Ai), 
a = a{0) ~ —0.14 ± 0.20, is correct, and that it is not possible to have a contribution 
to gi (x) from the P <S) P cut [HEX 73] . The reason is that only Regge poles or cuts with 
G{—l)^a = —1 can contribute to those virtual Compton scattering amplitudes that are 
relevant to gi{x). (Here T is the t-channel isospin, G is the G-parity and a the signature). 
However, it is possible to have a three-pomeron cut, but its contribution relative to ai- 
exchange is suppressed by both a factor (m/Q) (lnz/)~^ and a small numerical coefficient. 
Also we are reluctant to accept that the non- Regge singular behavior of F2{x) can be 
relevant at the values of involved in the EMC experiment. It seems generally accepted 
that the extrapolation to small x used by the EMC is reliable and is not the cause of their 
peculiar result. 

It is interesting that there is some support for the validity of the EMC extrapolation 
which comes from a totally different source, namely elastic up — > up. In the latter process 
the momentum transfer dependence of the differential cross-section gives information on 



the axial form factor Ga{Q'^) [AHR 87] whose value at = is directly related to the 
Qj appearing in DIS. One has 



^^(0) = ^|«3 + ^ as - ^ aoj . (4.2.24) 

It is not absolutely straightforward to obtain G'a(O) from the data, since it is some- 
what dependent on the form of the fit assumed for the variation in the measured range 
0.45 < < 1.05 (GeV/c)^. Thus, although no extrapolation in x is required, one does 
have to make an assumption of smooth behaviour in Q'^ as Q"^ ^ 0. The value of Ga{0) 
given in [AHR 87] corresponds to 

ao - V3 as = -0.45 ± 0.27 . (4.2.25) 

Using (4.2.14) and (4.2.21) yields 

ao = 0.14 ±0.27 (4.2.26) 

which is perfectly consistent with the EMC value in (4.2.22). 

The reader is warned that the up up data arc sometimes used as an argument 
in favour of a large polarized strange quark contribution As. This interpretation is only 
vahd in the Naive Parton Model where |[ao — -\/3as] = As. As will be discussed in Section 
8.5 there is also a gluonic contribution arising from the axial anomaly. 



4.2.4 The EMC data : higher twist effects 

The possibility that significant higher twist effects at the values of in the EMC 
experiment might invalidate a comparison of the data with the Parton Model or with the 
leading twist operator product expansion, was studied by Anselmino, loffe and Leader 
[ANS 89] in an approach based upon compatibility with the Gerasimov, Drell, Hearn sum 
rule. A more consistent treatment was later developed by Burkert and loffe [BUR 92, 93]. 

Higher twist effects have also been approached in a quite different way by Bahtsky, 
Braun and Kolesnichenko [BAL 90, 93] based upon the study of higher twist operators 
and QCD sum rules. The latter authors find much smaller corrections than in [ANS 89] 
and in [BUR 92, 93], but the uncertainties in their approach are intrinsically difficult to 
assess. Also there is some question as to the correctness of the whole approach [lOF 92] . 

We shall proceed to explain the two methods. Ultimately we feel that higher twist 
effects are unimportant in the EMC experiment but are surely not negligible in the new 
experiments at = 2 and 4.6 (GeV/c)^. 

a) The Gerasimov, Drell, Hearn sum rule approach. 

The key point is that is connected to the Gerasimov-Drell-Hearn (GDH) sum 
rule [GER 66, DRE 66] for forward scattering of real photons on nuclcons and that this 
indicates rapid Q^-dependences which might still have consequences at the EMC Q^- 
values. 

Define 

^ Vl{Q') . (4.2.27) 

Then at large Q'^ 

2A//"2 

hiQ") ^ -FW n{Q')As (4.2.28) 



where I'i{Q^)as is the asymptotic form given by say the lowest twist contribution to the 
operator product expansion and which is, empirically, positive. But, according to GDH, 



(4.2.29) 



where Kp is the anomalous magnetic moment of the proton {np = 1.79). Figure 4.4 shows 
Ip{Q'^) as given by DIS data and the value at Q'^ — 0. (Also shown is Ip{Q^) — In{Q^) 
relevant to the Bjorken sum rule). Since the Q^-dependence of r^(Q^)As, as given by 
perturbative QCD, is slow and logorithmic, the behaviour of I{Q'^) implies strong higher 
twist effects in T\{Cp). These were parametrized in [ANS 89] in a very simple way by 
putting 



1 



Cp jJL^ 



where 




(4.2.30) 



(4.2.31) 



and is a parameter setting the scale for the Q'^ variation. By vector dominance arguments 
it was suggested that [j? ~ m^. In (4.2.30, 31) F^^^ is taken as a constant, the limit of 
T\{Q^)as as g2 ^ cx). 

Equation (4.2.30) used in (4.2.27) provides a formula for r?((5^) for all 



r?(g^) 



J- l,As 



1 + ^^VQ^ 



1 - 



Cp fjj 



(4.2.32) 



and equating this to the EMC result at (Q^) = io.7 (GeV/c)^ yields T\ a^ ~ 1.2 T^mc^ 
i.e some 20% larger than the EMC result (4.2.19). However the Q^-depcndence implied 
by (4.2.32) is too rapid and was later shown to contradict the Q^-dependence of the data 
by 1.5 standard deviation [ASH 89]. 

More recently Burkert and loffe [BUR 92, 93] have argued that the contribution from 
resonance production has a strong Q^-dependence for small and then drops rapidly 
with and that this should be substractcd out before parametrizing the smooth large 
behaviour. The behaviour of the resonance contribution is shown in Fig. 4.5. Equation 
(4.2.30) is then modified to 



As 



where now 



1 + 



1/^W 



For the A (AT* (1238)) contribution they find [BUR 92] 

7^(0) = -0.78 

and including all resonances up to mass 1.8 GeV [BUR 93] 

7f,,(0) = -1.03. 



(4.2.33) 



(4.2.34) 



(4.2.35) 



(4.2.36) 



Bearing in mind that Kp/A — 0.80 we see that the resonance contribution has an 
important effect on the value of Cp. 



With this revised estimate of Cp, and using // = nip, (4.2.32) suggests that Ff^^ is 
about 8% larger than the EMC result. It appears that the Q^-variation implied by (4.2.32) 
might just be compatible with the EMC data divided into two bins with (Q^) = 4.8 and 
(g2) = 17.2 (GeV/c)2. 

For purposes of comparison with the QCD sum rule approach, let us expand (4.2.32) 
in inverse powers of Q"^ (we ignore IresiQ"^) at large Q^): 



1 + Cp)n'^ ^ 



The leading term in the higher twist correction is then 



(4.2.37) 



In summary the higher twist corrections, estimated on the basis of the GDH sum 
rule, are likely to be less than 8% for the EMC experiment. 

b) The QCD sum rule approach 

In the usual operator product approach explained in Section 5 one keeps only the 
operators of lowest twist r (r = 2 for the study of gi{x,Q^)). The approach of [BAL 
90], based on the QCD sum rule paper of Shuryak and Vainshtein [SHU 82], is a natural 
extension in which one studies also the operators of twist 4 and in which target mass 
corrections are taken into account. The problem is that one has to estimate certain matrix 
elements of gluon and quark condensates and certain correlators of quark and gluon 
operators. The latter is done assuming pseudoscalar meson dominance. 

The leading target mass correction SFit to Ff'" involves the third moment of Qi'^ 
and is given by [BAL 90, 93] 

^F?;^ = ^ldx grix, Q') (4.2.39) 

which yields 

.r;, . ^ (4.2.40) 

which is totally negligible in the EMC experiment. 

The higher twist corrections require an estimate of the nucleon matrix elements 
of certain twist-4 {T4) operators. The contribution to Ff involves flavour singlet (S) and 
non-singlet (NS) operators: 

Srlr, = -^2{l((0')) + (4.2.41) 

where [BAL 90, 93] 

8 r A/f ^ 

((0)) = -[((C/))-— ((V^))J (4.2.42) 

and {{U)), {{V)) are the reduced matrix elements of operator and V^i/,<t, defined by 

(P, S\U^\P, S) = 2MS^{{U)) (4.2.43) 

and 

(P, S\V^,JP, S) = 2M{{V)){(S^P, - S,P,)P,}s^^ (4.2.44) 



where S^a means symmetrization with respect to v^a. 
The actual operators are : 



[/5/Ar5 ^ g[uG,,Yu±dG,,^-d] (4.2.45) 
= g{uG,,^,u±dG^,^,d}s^^ (4.2.46) 



where G^,^ is the matrix 



= (y)- (4.2.47) 



It is a non-trivial matter, using sum-rule techniques, to estimate the nucleon matrix 
elements of these operators and there is some uncertainty about the results. Moreover it 
is claimed [lOF 92] that the handhng of the singlet is incorrect since the axial anomaly 
(to be discussed in Section 6) was not taken into account and certain correlators were 
estimated using pseudoscalar dominance by tt and r/ which is inappropriate to the singlet 
case. In addition note that there are errors in [BAL 90] which are corrected in [BAL 93]. 
In their corrected papers [BAL 93] one has 

_ (0.02 ± 0.013) {GeV/cf 



(0.005 ± 0.003) (GeV/c 
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5Vl^^ = -^^ (4.2.48) 

These corrections are much smaller than those obtained from the GDH sum rule 
(4.2.38) and are quite negligible for the EMC experiment. 

c) Summary on higher twist effects in the EMC experiment 

There is a significant difference between the size of these effects as calculated on 
the basis of the GDH sum rule (4.2.38) and on the basis of QCD sum rules (4.2.48). In 
both cases, however, the effects for the EMC experiment are negligible. 

We shall thus proceed to analyse the theoretical imphcations of the EMC result, 
which gives a value for oq compatible with zero. In the following we shall use the value 
given in (4.2.22). 

4.3 The EMC result: implications in the Naive Parton Model 

In Appendix B we show how to express the protonic matrix element of the ax- 
ial vector current of a quark field of definite flavour in terms of the polarized parton 
distributions, in the framework of the Naive Parton Model. 

We have 

(P,5|^/(0)7m75^/(0)|P,5) = 2MS, f dx [Ag/(x) + Ag/(x)] . (4.3.1) 

J 

It follows from (4.2.3, 5, 7 and 8) that 

ag = dx \l\uix) + /^uix) - Ad(x) - Ad(x)] (4.3.2) 

^0 



1 

08 = ^/ dx [Au(x) + Au(x) + Ad(x) + Ad(x) - 2As(x) - 2As(x)] (4.3.3) 
V3 ^0 

ao = AE= dx AHix) (4.3.4) 

^0 



and 



where 

AE(x) = Au{x) + Au{x) + Ad{x) + Ad{x) + As{x) + As{x) . (4.3.5) 

Notice that (4.3.2-4) inserted into (4.2.1) agree with (4.1.1). 

If, as is at first sight not an unreasonable assumption, one assumes that one can 
neglect the strange quark contribution so that Oq — V^^s one is led to the Ellis- Jaffe sum 
rule [ELL 74] mentioned in Section 4 which gives a much larger value for than found 
by the EMC: 

If 5 ^ 
where we have used (4.2.14, 16 and 21). 



(rf)^_^ = — <! aa + ^ as !> ~ 0.188 ± 0.004 (4.3.6) 



4-3.1 The "spin crisis in the Parton Model" 

Let us consider now the physical significance of AE(x). Since q^{x) count the num- 
ber of quarks of momentum fraction x with spin component ±| along the direction of 
motion of the proton (let us call this the 2;-direction) , the total contribution to coming 
from a given flavour quark is 

cixj (!),-(.) +(-l),-(a;)} 
= \ fdxAq{x). (4.3.7) 

z JO 

It follows that 

ao = 2Sf"'^^' (4.3.8) 

where is the contribution to from the spin of all quarks and antiquarks. 

Now in the Naive Parton Model p^ — Q and all quarks move parallel to the parent 
hadron, i.e. for a quark of momentum p, p — xP. Hence any orbital angular momentum 
carried by the quarks is perpendicular to P and thus does not contribute to J^. In addition 
it is assumed that the gluons are unpolarized, for reasons explained in Section 4.3.2. Hence, 
in the Naive Parton Model, one expects for a proton of helicity -1-1/2 : 

Sl^'^^^' ^ (4.3.9) 

We stress that this ignores effects and assumes only quark and antiquark constituents 
are polarized. 

Prom (4.3.8) and (4.2.22) we get instead of the value 1/2 of (4.3.9) 



(^gquarks\ ^ q ^ q Qg ^ q qq (4.3.10) 

/ Exp 



It was this highly unexpected result which was termed the "spin crisis in the Parton 
Model" [LEA 88] . 

4.3.2 The angular momentum sum rule 

In a more general framework, if we allow for the existence of parton transverse 
momentum and for the possibility that the gluons are polarized, (4.3.9) generalizes to the 
angular momentum sum rule 

= Sf""-"-^' + 5f + = 1/2 (4.3.11) 



and the question then arises as to whether the result (4.3.10) is as surprising as it at first 
seems. On the one hand models of the static properties of the baryons always have S'-wave 
ground states of quarks only, so in these there is no Sf^"'""^ nor /.po-^'o*^^*'. But these models 
deal with constituent quarks whereas the Parton Model involves current quarks. 

Unfortunately nobody knows the precise connection between these so it is perhaps 
unjustified to make assertions about partons on the basis of constituent models. However, 
Lipkin [LIP 90] has given an argument based on partons, which emphasizes how peculiar 
it is that Sf^"""^^^ is so small. We present a simplified version of his treatment. The key 
point is that the SU (6) prediction for (7^, namely 5/3, is really quite close to the measured 
value (4.2.16), so that the proton wave function cannot be too different from the SUiQ) 
wave function. Thus if we write for this proton state 

\p) = cos e \SU{<c>)) + sin 9 (4.3.12) 

where is a state orthogonal to |S't/(6)), the latter being a state in which the valence 
partons are in an SU{Q) configuration and anything else present has zero orbital, and zero 
total spin, angular momentum, then consideration of the contribution to (7^ from various 
possible plausible states IV'), each having different amounts of spin carried by the valence 
quarks and by any other constituents, leads to the bound sin^^ < 3/16. On the other 
hand in the SU{Q) state, = Sf^"''^'^^ = 1/2, so that requiring that 

1 cos^ e + (si'"''-'''] sin^ ^ ~ (4.3.13) 

2 V / ip 

implies the large value 

|^gnarfc.|^> 2. (4.3.14) 

Given the complexity of the state that this implies, it is then hard to understand why 
the proton is the only stable state with J = 1/2, / = 1/2. 

We shall be able to clarify this question a little more in Section 5 when we consider 
the surprising feature that QCD induces a Q^-dependence in Sf""^' and This 
will force us to rethink our interpretation of the physical meaning of the matrix element 

of ^5°;.- 

4-3.3 Trouble with the strange quark 

There is another difficulty associated with the results that follow from analysing 
the EMC and hyperon data in terms of the Naive Parton Model. 

We can use the numerical values (4.2.16, 21 and 22) to solve, via (4.2.14) and (4.3.2, 
3, 4), for the contributions of quarks and antiquarks of a given flavour. One finds: 

Au = t dx \Auix) + Auix)] = 0.79 ± 0.03 ± 0.04 (4.3.15) 

Jo 

Ad = I dx [Adix) + Adix)] = -0.47 ± 0.03 ± 0.04 (4.3.16) 
Jo 

As = [ dx [As(x) + As(x)] = -0.26 ± 0.06 ± 0.09 . (4.3.17) 

Jo 

But Preparata and Soffer [PRE 88] suggested that one can bound the strange quark 
contribution using the fact that, manifestly. 



\As{x)\ < s{x) 



(4.3.18) 



and the knowledge of the behaviour of s{x) as determined from deep inelastic neutrino 
experiments. 

There is an unjustified step in the analysis of [PRE 88] which can be modified [LEA 
88b] , and the argument is not rigorous since a specific form for As is assumed. Nonetheless 
the bound is unlikely to be far from the truth. It yields 



dx[As{x) + As{x)] 



< 0.072 ± 0.030 



(4.3.19) 



which is roughly a factor of three smaller than the value given in (4.3.17). 

We have to conclude that the analysis based on the Naive Parton Model or the 
Operator Product Expansion with free fields is in difficulty. 



4.4 Analysis of Fi using the QCD Improved Pcirton Model 

The analysis in Sections 4.1-3 was based on the Naive Parton Model which emerges 
from the Operator Product Expansion when all fields are treated as free fields. An imme- 
diate question is whether the difficulties we encountered can be removed by the inclusion 
of perturbative QCD corrections. In this Section we shall use the QCD-improved expres- 
sions which will be derived in Section 5 and demonstrate that the corrections do not help 
significantly. We shall also discuss the Bjorken sum rule. 



4-4-^ The operator product expansion for 

As is well known, QCD corrections induce a logarithmic Q^-dependence which 
breaks Bjorken scahng. As is outlined in Section 5.2 the expression (4.2.1) is modified 
to 

r?(Q^) = ^{ («3 + ^ as^NsiQ') + I aoEsiQ')"^ (4.4.1) 

where the coefficient functions Ens and Es have perturbative expansions [KOD 79] re- 
cently extended to order and respectively [LAR 91, 94] 

Ens{Q') = 1 - f - ( ) I ^ ) - I 7oZ ) I ^ ) (4-4.2) 

^^(^') = l-(o:040 -;)(^) (4-4.3) 

where = «s(Q^) and where the upper and lower figures refer to Nf = 3 or 4 respectively. 
Strictly speaking the coefficients of and as are only equal in the case of massless quarks, 
which is assumed in the calculation of (4.4.2). These perturbative corrections will be 
important in analysing the new data at lower Q^, to be discussed in Section 4.5.2. In the 
notation of Section 5.2 we have E^s — E} -^ — E}^; Es — El,^. 

Once interactions are allowed the operators and currents have to be renormalized 
and their matrix elements will, in general, depend upon the subtraction point or renormal- 
ization scale fi. For a conserved current, however, the matrix elements are independent of 
the renormalization scale. These matters were mentioned in Section 4.2 and will be taken 
up again in Section 6. The important point to remember is that a^ and ag are independent 
of II, whereas Oq depends upon /i. 




If we choose = then, in leading logarithmic approximation (LLA), the ex- 
pressions (4.3.2-4) for the in terms of quark densities Ag(x) are simply modified by 
the replacement 

Ag(x) ^ Ag(x;g2) (4.4.4) 

where the Q^-evolution is controlled by the spin dependent Altarelli-Parisi equations 
[ALT 77]. 

For the EMC experiment the mean value of is about 10.7 (GeV/c)^ which 
corresponds to ^^(Q^) ~ 0.24 for 4 flavours. The leading correction term is of order 8% in 
(4.4.2) and even smaller in (4.4.3 and we cannot expect a dramatic change in the values 
quoted in (4.2.22). Indeed one finds 

ao = 0.17 ±0.12 ±0.17 (4.4.5) 

to be compared with the naive expectation ao — 1. The values quoted in (4.3.15-17) 
become, for ~ 10 (GeV/c)^ 

Am = 0.82 ±0.03 ±0.04 

Arf = -0.44 ±0.03 ±0.04 (4.4.6) 
As = -0.21 ±0.06 ±0.09 

which are close to the values given in (4.3.15-17) Q. 

Thus the straightforward perturbative QCD corrections to (4.2.1) in no way alleviate 
the problem, and the value of As is still surprisingly large compared with (4.3.19). 

4.5 The new experiments: neutron data and the Bjorken sum rule 

The operators Jg^ and that give rise to the terms oq and ag in (4.4.1) are 
invariant under isotopic spin rotations. Thus in going from proton matrix elements to 
neutron matrix elements they remain unchanged. J|^, on the other hand, transforms like 
the 3rd component of an isotopic spin triplet and therefore changes sign. It follows that 
the Bjorken sum rule 

dx \gl{x,Q'')-gl{x,Q'')\ = 5b,(Q') (4.5.1) 

^ 

where SBjiQ"^) is the theoretical value of the Bjorken sum rule, given by 

^,o.™a.o.ooo4-^-(-)(^y-(--)(^)'...), 

holds on fundamental grounds. Note that the fact that 03 is independent of is linked 
to the conservation of the non-singlet axial current and is discussed more fully in Section 
6.1. The upper and lower figures in (4.5.2) refer to Nf = 3 or 4 respectively. 

^) Note added in proof: An analysis of all the proton data at = 10 (GeV)^ [ADA 94] alters (4.2.19) 
to = 0.142 ± 0.008 ± 0.011 and consequently (4.4.5) becomes (notice that in [ADA 94] as ~ 0.23 
is used) ao — 0.27 ± 0.08 ± 0.10, still uncomfortably far from oq = 1. 



The Ellis- Jaffe sum rule value (4.3.6), ~ 0.188, would imply, via (4.5.1, 2) at 
g2 = 10.7 (GeV/c)2 

= / dx g^(x, Q^) ~ -0.0003 

Jo 

whereas the measured EMC value of Ff gives hope of a somewhat more sizeable neutron 
result: 

r^* ~ -0.06 at = 10.7 (GeV/cf . 

4-5.1 The new experiments on deuterium and ^He 

The first experiments to measure gi{x) have now begun to report rcsTilts. At CERN 
the Spin Muon Collaboration (SMC) have used a polarized deuterium target at an average 
g2 of 4.6 (GeV/c)2 [ADE 93]. At SLAC the E142 group uses a polarized ^He target at 
(Q2) = 2.0 {GeV/cf [ANT 93]. (For a combined analysis of the data see [ADE 94]). 

The SMC data on the deuteron asymmetry Af{x) at = 4.6 (GeV/c)^ are shown 
in Fig. 4.6, and in Fig. 4.7 the data on xgf{x). The E142 data on the neutron asymmetry 
A'^[x) and on gi{x) at = 2 (GeV/c)^ are shown in Fig. 4.8. The neutron data were 
extracted as explained in Section 2.1.4 using the EMC proton data. 

Soon after the appearance of the new results some rather dramatic claims were 
made that the very fundamental Bjorken sum rule had been violated. In retrospect these 
incorrect assertions arose because of too naive a comparison of the new data with the 
EMC proton data at 10.7 (GeV/c)^. A more careful analysis is necessary, as emphasized 
by Ellis and Karliner [ELL 93] and by Close and Roberts [CLO 93] and AltareUi, Nason, 
Ridolfi [ALT 94]. 

The main issues are the following: 

1) The extraction of gi{x,Q'^) from nuclear data at some (Q^). as explained in 
Section 2.1.4, requires a knowledge of gi{x, Q^) at the same mean Q^. It is unjustified to 
naively combine the EMC data at (Q^) = 10.7 (GeV/c)^ with the deuterium data at 4.6 
and the ^He data at 2.0 {GeV/cf. 

There is no absolutely safe way to overcome this difficulty. It does seem, however, 
that the values of the asymmetries at fixed x, show no significant variation with where 
it has been possible to examine this. On the other hand the variation of F2{x,Q'^) with 
has been much studied. 

In the improved analysis, therefore, gi{x, Q'^) is calculated at the required Q'^ by 

using 

,^(x O^) - ^ ^^'("'Q^) (4 5 3) 

and assuming A^^ independent ofQ^. The resulting dependence of gi{x, Q"^) upon the mean 
is shown in Fig. 4.9 taken from [ELL 93]. In (4.5.3) values of were taken from the 
NMC data [AMA 92] and of RP from the SLAC data [WHI 90]. These g\{x,Q'^) yield, 
[ELL 93], 

T\[2{GeY/cf] = 0.124 ± 0.013 ± 0.019 
rn4.6 (GeV/c)^] = 0.125 ±0.013 ±0.019. (4.5.4) 

Figure 4.10 shows the gi{x, Q'^) extracted from the nuclear data on this basis. 

2) The extrapolation of gi{x, Q"^) to x = 1 and x = in order to compute the first 
moment F" requires some care. 



The SMC experiment at (Q^) = 4.6 (GeV/c)^ covers the range 0.006 < x < 0.6. 
They find 

/ ' gf{x) dx = 0.024 ± 0.020 ± 0.014 . (4.5.5) 

J0.006 

The extrapolation to a; = is taken to be of the reasonable form gf ~ a;~" with 
—0.5 < a < 0, in accord with the discussion in Section 4.2.3 and compatible with Eq. 
(4.2.23). This leads to an estimate of -0.003 ± 0.003 for the integral from to 0.006. 

The extrapolation to a; = 1 is based on the fact that the asymmetry is bounded, 
\A\ < 1, and leads to the estimate 0.002 ± 0.004 for the integral from 0.6 to 1. 

The result is, per nucleon, 



rf [Q2 = 4.6 {GeV/cf] = 0.023 ± 0.020 ± 0.015 . 
Note that using (2.1.55), (4.4.2, 3) and (4.2.14) we have (taking Nf = 4) 



(4.5.6) 
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(4.5.7) 



Taking = 0.28 ± 0.03 at = 4.6 (GeV/c)^ [NAS 93] and using the values in 
(4.2.21), one obtains {ud = 0.058) 



ao = 0.09 ±0.20 ±0.15, 



(4.5.8) 



again consistent with zero and perfectly compatible with the value quoted in (4.4.5). 

For [ELL 93] obtain a result very slightly different from the value quoted by the 
SMC because of the small change in shown in (4.5.4), i.e. 



F^ [4.6 (GeV/c)^] - -0.076 ± 0.037 ± 0.046 . 



(4.5.9) 



The E142 experiment at (Q^) = 2 {GeV/cf covers the range 0.03 < x < 0.6. 
[ELL 93] have re-analysed the extrapolations to x — and x — 1 using the same approach 
as in the SMC case and they estimate contributions to F" of —0.006 ±0.006 for the integral 
from to 0.03 and 0.000 ± 0.003 for the region 0.6 to 1. The result is 



F^ [2.0 (GeV/c) 



-0.028 ±0.006 ±0.009 



(4.5.10) 



which differs somewhat from the value —0.022 ± 0.011 quoted by E142. 



4.5.2 Tests of the Bjorken sum rule 

Using the values of Ff given in (4.5.4) and of F^' given in (4.5.9 and 10) one can 
now test the Bjorken sum rule, bearing in mind, however, that we have assumed that the 
asymmetry is Q^-independent for the range involved. Before comparing with the data it is 
instructive to see to what extent the perturbative QCD corrections influence the expected 
value of the Bjorken sum rule, SsjiQ^) given by (4.5.2). Wc show below in Table 4.1 the 
values of (Q^) various calculated to order as, «^ and starting from the central 
value 0.2096 without corrections. We have used values of as{Q^) from [NAS 93] which 
differ somewhat from the values used in [ELL 93] and also from those used by the SMC 
group. We believe these represent a good compromise between the LEP and unpolarized 
DIS determinations which differ somewhat. 



Table 4.1 Values of q;s(Q^) and of SsjiQ'^) to various orders in perturbation theory 









(0^) = 


= 0.2096 to zero order 


in (GeV/c)2 


Nf 




to 0{as) 


to 0{al) 


to 


2 


3 


0.354 ±0.04 


0.186 


0.176 


0.170 


4.6 


4 


0.283 ±0.03 


0.191 


0.185 


0.183 


10.7 


4 


0.238 ±0.03 


0.194 


0.190 


0.189 



We see that at = 2 the corrections are very important, ~ 20%, and that even at 
— 10.7 they are not neghgible, ~ 10%. 

We come now to the comparison between the experimental and theoretical values of 
the Bjorken sum rule. In allowing for target mass corrections and higher twist corrections 
as discusses in Section 4.2.4 we must recall that the GDH and QCD sum rule approach 
give quite different estimates for the higher twist corrections. 

In general terms one has 

/' dx [g{{x, Q^) - Q2)] = SB,m + SS^i + SS?^ . (4.5.11) 
Jo 

The target mass correction given in (4.2.40) is neghgible for the proton and an 
estimate for the neutron is much smaller [ELL 93] so we can ignore even at Q = 2. 
In the QCD sum rule approach one has, from (4.2.48) 

(4.5.12) 



(0.015 ±0.013) (GeV/c) 



QCD sum rules 

In order to estimate the higher twist corrections in the GDH sum rule approach let 
us rewrite (4.2.37) as 



i 1 [W ) — ^ l,As ~ n2 hAs + ^P,n 



where 



2 1m2 



2 



p,n 



C(o) 



Using the values [BUR 93] 



I^^%0) = -1.028 



7^^^(0) = -0.829 



(4.5.13) 
(4.5.14) 

(4.5.15) 



we find 



GDH sum rule 



= -^|2(rf-r^)^, + Ap-A, 

0.16(GeV/c)^ 



(4.5.16) 



This appears to be a huge correction at small values of Q^. However it is illegitimate 
to use (4.5.16) in a region where the resonance contribution lRes{Q^) in (4.2.33) has not 
yet died out. The (p — n) curve in Fig. 4.5 indicates that Ir^s is not negligible &t Q'^ — 2 
(GeV/c)2. Thus we should add a contribution 

'"-""{Q^) -0.05 (4.5.17) 



2M2 



)2=2 



to the correction (4.5.16) at = 2 (GeV/cf, yielding 



= -0.03 



at 



Q2 = 2 {GeV/cf 



(4.5.18) 



GDH sum rule 



At Q'^ — 4.6 (GeV/c)^ it would appear to be safe to neglect lRes{Q^) and to use just 
(4.5.16). However the corrections do seem surprisingly large. We are unable to estimate 
the errors in (4.5.16 and 18). 

In Table 4.2 we compare the experimental value of the LHS of (4.5.11) with the 
theoretical value of the RHS using the results of Table 4.1 and with the two versions of 
the higher twist corrections. 

Table 4.2 Test of the Bjorken sum rule 



in (GeV/c)2 




Sb,{Q^) + SSt' +6S^i 


QCD sum rule 


GDH sum rule 


2 

4.6 


0.152 ± 0.014 ± 0.021 
0.201 ± 0.039 ± 0.050 


0.162 ±0.023 
0.180 ±0.017 


0.14 
0.15 



It is seen that there is absolutely no significant evidence for a failure of the sum rule. 
Clearly at the low value of = 2 the higher twist effects help to make the agreement 
more impressive, but even without them there is no real contradiction, given the size of 
the experimental errors. 



5 The Operator Product Expansion (OPE) 

The fundamental understanding of the Q^-behaviour of the moments in unpolar- 
ized DIS came originally from a study of the operator product expansion. Later it was 
discovered that the same results could be obtained in the QCD Improved Parton Model. 

When it comes to the polarized case it has been claimed that the two approaches 
yield different results. We shall argue that this is not the case and claim that there is 
perfect agreement between the two methods. 

We shall also clarify the situation with regard to g2{x) itself, to its first moment 
(the Burkhardt-Cottingham sum rule), and to its relationship with gi{x) (the Efremov- 
Leader-Teryaev and the Wandzura-Wilczek sum rules), about which there are many 
misconceptions in the literature. 



5.1 General structure of the OPE 

It is well known, and can easily be deduced from the Feynman diagram Fig. 2.1, 
that the hadronic tensor W^^'^ involved in the expression for the deep inelastic cross- 
section [sec (2.1.1)] is given by the Fourier transform of the nucleon matrix elements of 
the commutator of electromagnetic currents Jfj,{x): 

W,M,P,S) = ^I d'xe''^-'{P,S\[J,{x),MQ)]\P,S) (5.1.1) 

where S'^ is the covariant spin vector specifying the nucleon state of momentum P^. It is 
convenient to introduce an amplitude T^^, which is closely related to the forward T-matrix 



element for Compton scattering of a virtual photon of 4-momentum q and helicities A and 
A'. In the convention for the T-operator 



5 = 7 + i(27r)^ 5\Pf-Pi) T 

one has 

(P, S- q, X'\f\P, S; q, A) = ATra e;{X') T^^" e,{X) (5.1.2) 

where 

T^,{q;P,S)=i J d^x e"i-^{P,S\T{J^{x)MO))\P,S) (5.1.3) 

is given in terms of the time ordered product of the currents. 

Both Wj^i, and Tj^^, may be split into parts symmetric (S) or antisymmetric (A) 
under /i v, and one can show that 

TT 

The symmetric part is independent of the spin vector S^^ and plays a role only in 
unpolarized scattering. We shall therefore concentrate on the antisymmetric part and we 
follow the treatment of Kodaira and his group [KOD 79]. 

The behaviour of T^j, (and therefore Wf^iy) in the deep inelastic limit is controlled 
by the behaviour of the product of currents near the light cone x^ — and can be derived 
from Wilson's operator product expansion. 

It is important to note that the expressions (5.1.1 and 3), for which the operator 
product approach can be utilised, only arise because of the fully inclusive nature of the 
deep inelastic reaction being considered. Indeed the starting point from which (5.1.1) can 
be derived is the expression 

ly^, oc Y,{P,S\J^\X){X\J,\P,S) 

X 

which appears in the formula for the cross-section. Only if the sum is over all final states 
\X) does this reduce to (5.1.1). 

The antisymmetric part of the Fourier transform of the operator product appearing 
in (5.1.3) is expanded in terms of local operators R and coefficient functions E in the 
form [KOD 79]: 

[1 - (-1)"] / 2 \ " 



:jd'xe''^-^T{J,{x)JM) = -^E^ 2^(g^) ^'^--^M„-.E^^{^.^A.?V-i 

X El^{Q\g)K^-'-' (5.1.5) 

r7 - 1 

+ (s^pAa q^q" - £vp\a qp.q" - q e^wXa) 



n 

tl.../J.n-2 1^ 

'2,1 



where g is the QCD coupling constant and where i takes on the values l,2...,8,ip,G, the 
detailed significance of which will be explained below. The operators with i = 1,...,8 
transform like an SU(3) flavour octet. Those with labels ip or G are flavour singlets. 

The factors 5i reflect the charge and isotopic spin structure of the currents. Taking 
for the electromagnetic current (in units of e) 

•^M = E f = u,d,s (5.1.6) 



(a colour sum is, of course, implied), and bearing in mind that we are ultimately only 
interested in the current commutator that occurs in W^,^, one may take for the 5i 

Si = - 6s = S^ = 6g = - (5.1.7) 

3 ^ 3^3 ^ 9 ^ ^ 

and all other 5i — 0. 

Note that, as will be seen below, there does not exist an i?2-type operator with 
n < 2. Nonetheless because of the explicit factor (n — 1) in (5.1.5) we have formally 
written a sum from n = 1. Also there is no operator of the Ri^g type for n = 1. It should 
formally be regarded as zero for n = 1 in the sum over n. 

The leading twist operators occurring are of two types: 

1) The set •••'^"-1 of twist 2. 

For i = 1, 8 these are the flavour non-singlet operators 

^ (.)n-i I^^j ^1 > 1) (5.1.8) 

where S implies complete symmetrization in the indices o", /ii,/i2, ...,/i„_i, the Aj are the 
Gell-Mann SU{3) flavour matrices and D^^ is the usual QCD covariant derivative. For 
i — ip OT G we have the flavour singlet operators 

^ {i)-^{^j,YD'''-D''-'^}s {n>l) 
^aMi...Mn-i ^ {i)^-'Tt{s'"'f^^Gf3^D''K..D''-'G^"-'}s (n > 2) (5.1.9) 

where G^i^ is the usual matrix (in colour space) form of the gluon field tensor [see Eq. 
(6.3.5)]. 

2) The set i^g^'i-'^-^ of twist 3. 

For i = 1, 8 these are the flavour non-singlet operators 

J^^^.^...^r.-2 ^ |^^5^A^a^Mi_ _^Mu-2 l^^j ^1 (^ > 2) (5.1.10) 

where 5" implies antisymmetrization on A, a and symmetrization on the indices /xi, .., //„_2. 
For i — or G we have flavour singlet operators 

^Aji...M„-2 ^ (^)"-i{V;757^D'^D^\..D^'^-=^V}5' {n>2) 

^Aji...,.„_2 ^ (i)"-iTr{£""'^^G;3^D^\..D^"-^G^}5' (n>2). (5.1.11) 

The above operators are often referred to as i?2 i brevity. 

It should be noted that the above list of operators is a complete set of operators in 
the massless quark theory (if rriq ^ there exists a further set of twist 3 operators). We 
shall not deal with these. They are discussed in [KOD 79]. 

The coefficient functions El'^{Q'^,g),E2i{Q'^,g) are the Fourier transforms of the 
singular functions of x that appear in the Operator Product Expansion and have to be 
calculated, using perturbative QCD, as a power series expansion in the coupling g. The 
factors in (5.1.5) are arranged so that for free fields, i.e. g — 0, one has the simple results: 

= 1 E^.{Q'',g^ 0)^1 i = l,...,8 

E^^^iQ^g^O) = 1 E^^^{Q\g^ 0)^1 (5.1.12) 

EIg{Q',9-0) = Ela{Q',9^0)^0. 



Since gluons are electrically neutral it is no surprise to see that the operators in- 
volving the ghion field tensors play no role when g = 0. 

In an interacting theory the i?'s are functions, Ei^{Q^ / fi'^; g) etc., of the renormal- 
ization scale /i^ and they satisfy standard renormalization group equations. The flavour 
octet operators are multiplicatively renormalized, but the G and type singlet operators 
mix under renormalization, except for the case n — 1. This is reflected in the Q^-evolution 
of the coefficient functions (for details see [KOD 79]). 

The proton matrix elements of these operators have the form: 

(P,5|i?i^f •••^"-^|P,5) = ^{S''P''\..P''-'], (5.1.13) 

77/ 

{P,S\E^J'-''^-^\P,S) = Mdi^iS^P^ - S^P'')Pi'K..pi'"-\ (5.1.14) 

where the factors a^, dl^ - essentially reduced matrix elements - reflect the unknown, non- 
perturbative aspect of the dynamics. Depending upon which current is involved, they may 
or may not depend upon the renormalization scale /i^, as is discussed later. 

The connection with the simpler operators used in Section 4 is as follows. For n — 1 
we have 

(i?i,,). = V'TsTafy)^ i = l,...,8 (5.1.15) 



so that via (4.2.7) and (5.1.13) 
Also, note that 



-Jl of eq. (4.2.3) 
Qi = 2a\ . (5.1.16) 



{Rl,^)a = tplbla-ip (5.1.17) 

= -J°, of eq. (4.2.5) 



so that via (4.2.8) and (5.1.13) 



ao = af. (5.1.18) 



5.2 Equations for the moments of gi,2{XjQ^) 

The expressions for the moments are obtained from the OPE in the standard fashion: 

i) 511,2 are related to the absorptive parts of forward scattering amplitudes for virtual 
Compton scattering; 

ii) dispersion relations in the energy 1/ are written for these scattering amplitudes; 

iii) for small enough 1/ the latter can be expanded in a power series in 1/ or, equivalently, 
in 1/x. Only even or odd powers appear because of the crossing symmetry of the scat- 
tering amplitudes. The coefficients in the power series are moments of the absorptive 
parts and thus of the scaling functions gi^2', 

iv) this power series in 1/x is matched to the one which emerges from the OPE when 
the scalar products occurring in (5.1.5) are multiplied out. 

The result is an expression for the moments of gi^2 in terms of the factors a^, 
and the coefficient functions ^'"j and £^2r C)ne finds [KOD 79] for the n-th moment: 



L 



' dx x'^-'g.ix, Q^) = \Y. 5i<KiiQ^^ 9) n = 1, 3, 5 . . . (5.2.1) 
1 ~^ ' 



and 



i 



n = 3, 5, 7 . . . 



(5.2.2) 



Note firstly that the Operator Product Expansion only gives information about the 
odd moments of gi,2{x, Q"^). 

Note also the important feature that (5.2.2) only gives information about the mo- 
ments of g2 for n > 3. In the original literature (5.2.2) is given as holding also for n — 1. 
That this is incorrect has been stressed by Jaffc [JAF 90], and we support his contention. 
Equation (5.2.2) thus gives no information about the first moment 



Let us now concentrate on the first moment, n = 1, of gi in Eq. (5.2.1). Since, from 
(5.1.9), there does not exist a gluon operator for n = 1, the first moment ri(Q^) of gi does 
not receive a direct contribution from gluonic operators. In fact in the free theory, with g = 
0, there is no gluonic contribution at all to gi{x) or g2{x) as explained after Eq. (5.1.12), 
and the constants a^, of (5.1.13, 14) are then independent of //^. Now because as{Q'^) — > 
as (5^ — > oo it often happens that results found in the interacting theory reduce to the 
structure of the free theory results with just the replacements a^(Q^), — >• dl^{Q'^) 

equivalent, in the Parton Model, to the replacements q{x) q{x, Q'^). If this were always 
the case then we would expect no gluon contribution to gi{x,Q'^) as —>■ oo. But 
as already mentioned in Section 4, and as will be discussed in much greater detail in 
Sections 6, 7 and 8 in the framework of the Field Theoretic Parton Model, there is a 
contribution to gi{x,Q'^) from gluons in the nucleon even in the limit Q"^ oo. This 
apparent disagreement between the two approaches has upset many people, leading them 
to question the validity of the Field Theoretic Parton Model treatment. It must be stressed 
that there is no disagreement. The limit g ^ and the limit Q"^ ^ oo are not equivalent. 
The Operator Product Expansion will give a gluonic contribution to ri(Q^) even though 
no gluonic operator exists for n = 1. The point is that the gluonic contribution emerges 
from the quark operators when account is taken of the fact that they are interacting fields 
and not free fields. There is, in fact, complete equivalence between the approaches as is 
explained in detail in Section 8.6. 

When the detailed results [KOD 79] for the coefficient functions calculated to order 
as are fed into (5.2.1) one obtains the fundamental formula (4.4.1) which was the basis 
for the QCD improved analysis of the polarized DIS data and which leads to the Bjorken 
sum rule (4.5.1, 2). 

5.3 Is there a connection between gi{x) and g2{x)l 

There is a much-quoted sum rule, known as the Wandzura-Wilczck relation [WAN 
77] which claims to relate g2{x) to an integral over gi{x). The operator product results 
(5.2.1, 2) allow us to see immediately that the relation cannot be exact. 

From (5.2.1) and (5.2.2) we have that 




(5.2.3) 




(5.3.1) 
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Y^kd^nEUQ^g) (n = 3,5,7...) 



+ 



2n 



where, recall, the contribution of the last term on the RHS is of twist 3. 

Now, if we assume that (5.3.1) is valid for all integer n - and we shall see in a 
moment that this may be a dangerous assumption - we can use the convolution theorem 
for Mellin transforms (which relates the product of the moments of two functions to the 
moment of their convolution) and Eq. (5.2.1) to obtain 

g,(x, Q^) + g^ix, Q^) ^ [' ^ g^{y- Q^) + [twist 3] . (5.3.2) 

Jx y 

In [WAN 77] it is argued, on the basis of a model, that the twist 3 term can be 
neglected, leaving a direct connection between gi and g2- However, this argument is quite 
unreliable. The selfsame model gives nonsensical results for ^1,2(2;) and gi{x)\ 

Thus, as will be further explained in Section 5.4, the neglect of the twist-3 term in 
(5.3.2) is dangerous and the Wandzura-Wilczek relation should not be trusted. 

But even (5.3.2) may be a dangerous expression. For recall that (5.3.1) does not 
hold for n = 1. In extrapolating it to n = 1 we are tacitly assuming [see (5.2.2)] that the 
first moment of 5'2,r2 = lo dx g2{x), vanishes, whereas, as stressed earlier, the Operator 
Product Expansion gives no direct information on r2 (this assumed vanishing of r2 is 
known as the Burkhardt-Cottingham sum rule [BUR 70] and will be discussed in the 
next Section). Thus even (5.3.2) is, a priori, not reliable. 

The correct way to proceed is to simply use (5.2.2) with (5.2.1) for n > 3, i.e. 

£ dx X--' [g,{x) + (-^) g^ix)^ = ^ ^ S,dl,El^{Q', g) (n > 3) (5.3.3) 

where the RHS is a twist 3 contribution. 

This is still a useful equation, for once gi{x) and g2{x) are measured it yields 
information about the twist 3 operator contribution, which, it turns out, depends upon 
correlations amongst the partons in the nucleon - a question of great importance because 
it is just what determines single spin asymmetries in QCD (see Section 10). 

5.4 Does the first moment of g^ix) vanish? 

We have already mentioned that according to the Burkhardt-Cottingham sum rule 
[BUR 70] 

r2(g')= f'dxg2{x,Q^)^0. (5.4.1) 
Jo 

We now wish to ask whether or not we can expect this sum rule to hold. Firstly, 
as stressed in the previous section, the result does not follow from the Operator Product 
Expansion. In fact the OPE gives no information about the first moment of g2 since (5.2.2) 
is not valid for n = 1. 

The result (5.4.1) was derived originally by considering the asymptotic behaviour 
of a particular virtual Compton helicity amplitude whose absorptive part is proportional 
to g2, and rests on the assumption that the asymptotic behaviour is controlled by low- 
lying Regge poles - in this case Regge poles hke the ai with aai{0) ~ —0.14. But a very 
careful analysis by Heimann [HEI 73] showed that there should be contributions from 
multi-pomeron and pomeron-Regge cuts thus invalidating the starting assumption in the 
Burkhardt-Cottingham derivation. 

It should be noted that this same argument invalidates the Wandzura-Wilczek 
result with [twist 3] put equal to zero on the RHS of (5.3.2). For if (5.3.2) were true 
without the twist 3 term one could deduce (5.4.1) from it directly by integrating over x. 



Because of the well known connection between the behaviour of scaling functions 
in the limit x —>■ and asymptotic Regge behaviour, the coupling to multi-pomeron cuts 
implies a highly singular behaviour 

92{x) \ (5.4.2) 

where we have ignored factors of Inx. Thus the integral in (5.4.1) might not even converge! 

That said, we should stress that the couplings of Reggc poles and cuts is a sub- 
tle question and such a remarkable conclusion as (5.4.2) must be tested experimentally. 
Because Regge cuts always arise from non-planar Feynman diagrams, and since these cor- 
respond to twist greater than 2, it might be that the coefficient of in (5.4.2) tends 
to zero as — >■ oo. It is then not clear at what scale of fixed one would expect to 
observe the divergence experimentally. Alternatively, if g2{x) C{Q'^)/x'^ as a; — * with 
C{Q^) ^ 1, it is not guaranteed that the divergence will be seen experimentally if x is 
not small enough. 

Heimann has given a clear illustration of the mechanism which could be behind the 
failure of the Burkhardt-Cottingham sum rule. It involves an illegitimate interchange of 
double integrals in a Fourier transform. 

One can show that g2{x) is given by a Fourier transform of the following type 

/oo p—iXx 
_^c?A^[A/(A)] (5.4.3) 

where the behaviour of /(A) as A —> is such that A/(A) — > as A — > 0. 

Now because of its relation to the absorptive part of a scattering amplitude one 
knows that g2{x) — for \x\ > 1 and g2{—x) — g2{x)- Then 

/ dx g2{x) cos iJ,x = - dx g2{x)e^'^^ 

Jo 2 J-oc 

1 PCX) POO rj\ 

= M dx '^e^^^-^>[\m\. (5.4.4) 

//the integrals can be interchanged, we have 
-■i , , 1 



/•I ^ POO 

I dx g2{x) cos fix — - I dX 6{ijl — A) [A/(A)] 

2 J-OD 



^/^/(/^) (5.4.5) 



so that taking fj, — 



C dx g2(x) = . (5.4.6) 
Jo 

However it is NOT always legitimate to interchange the integrals, and this is obvious 

if the integral on the LHS of (5.4.5) diverges, say for /x = 0. Indeed the Regge behaviour 
of g2{x) as X — > is precisely an example of this type. 

If one takes 

/(A) = A"-^ (0 < a < 1) (5.4.7) 

so that indeed A/(A) — > as A — > 0, one finds upon taking the Fourier transform (5.4.3), 
[HEI 73] 

92{x) oc ^ (5.4.8) 

and the integral of g2{x) diverges. 



5.5 The Efremov— Leader— Teryaev sum rule 

If one assumes that (5.3.3) holds also for n = 2 one obtains 



£ dx x[g,{x) + 2g^{x)] = ^ S^EUQ^ q) 



(5.5.1) 



To lowest order in a^, because of (5.1.12) we need keep only the operators with i = 1, 8 
and i — '^.ln that case, from (5.1.14) and (5.1.10 or 11) the d2 are of the form 



where is a four-vector introduced in the Field Theoretic Parton Model (Section 10) 
such that = and n - P —1. 

It can be shown, in the framework of the Field Theoretic Parton Model, that the 
RHS of (5.5.2) vanishes in the chiral limit. This follows upon integrating cq. (D.20) of 
Appendix D and using eqs. (D.17, 13 and 6). Hence one has the sum rule [EFR 84] 



This is further discussed in Appendix E. 

It should be noted, however, that like the Burkhardt-Cottingham sum rule, this 
result does not follow strictly from the OPE. Moreover if g2{x) diverges hke, or close 
to, 1/x^ as X — > then, as with the Burkhardt-Cottingham case, the sum rule will fail 
because of the divergence of the integral of xg2{x). 

6 The axial anomaly and the axial gluon current K^^ 

As we have seen in the previous Sections the Naive Parton Model interpretation 
of the EMC experiment, based upon treating the flavour singlet axial current 7°^ as 
effectively the quark spin-density operator, leads to a negligible contribution to the proton 
spin from the quark spins, oq = 2S'f""'''^* ^ 0, in contradiction with the Quark Model. The 
fact that the experiment was done at an average ~ 10 (GeV/c)^ and that the Quark 
Model is usually trusted for static properties of the nucleon i.e. for small momentum 
transfers, ~ 0, does not mean that the Naive Parton Model and Quark Model results 
cannot be compared with each other, as is sometimes claimed, so that there is a real 
contradiction. The statement that in the Quark Model we deal with constituent quarks 
and not with partons is somewhat irrelevant, since the states of three constituent quarks 
can be considered as a sum of states involving three valence partons in the same quantum 
state as the constituent quarks, with all the other partons (gluons, qq pairs etc.) having 
zero total angular momentum, as emphasized in Lipkin's argument in Section 4.3.2. 

Secondly, there appears to be no significant Q^-variation in the experimental data 
down to (5^ — 1 (GeV/c)^. The sharp change in T\{Cp) demanded by the Drell, Hearn, 
Gerasimov sum rule occurs for smaller (as discussed in Section 4.1.1), so cannot be 
used, in the region of the deep inelastic data, as an argument for large higher twist effects. 

We must, therefore, return to the foundations of the Parton Model to try to seek a 
resolution to the problem and we shall learn, surprisingly, of subtle effects which force us 
to alter some aspects of our naive interpretation of the Parton Model. 



d2 oc (P, S\il>^S{S ■ D) - ${n ■ D)]il^\P, S) 



(5.5.2) 




(5.5.3) 



6.1 Is there really a "spin crisis"? 



In field-theoretical language, the usual statement that the nucleon predominantly 
consists of three quarks in an SU{6) state (with zero orbital angular momentum) means 
that the 3 quarks —>■ 3 quarks Green function with these quantum numbers has a pole 
at the mass of the nucleon, with J = 1/2, with a significant residue. 

The fact that the nucleon is not just a 3 quark state means that the nucleon pole 
also occurs, albeit with smaller residue, in more complicated Green functions such as 
Sq + g — > 3q + g, perhaps with non-zero orbital angular momentum. The operator J"^, 
containing quark field operators only, is the quark spin density operator for free fields and 
therefore it does measure the total quark spin in the initial parton state. But since, as we 
shall see, J^^ is not a conserved current the expectation value in the nucleon state where 
the partons interact strongly with each other is not the same as the expectation value 
in the initial state made up of free partons. Only in the Naive Parton Model, where we 
ignore these interactions, can we equate the expectation values in the nucleon and free 
parton states. Since everybody knows that only total and not Sz is generally conserved 
one may wonder why it is even necessary to emphasize the non-conservation of J^^. The 
reason is that naively , i.e. from the equations of motion, J^^ is conserved for massless 
quarks. So the non-conservation is anomalous if one works with massless quarks as is often 
done in the Parton Model. 

For and the conservation is exact with massless quarks. It is also exact if 
SU (2) and SU (3) are exact symmetries. Since these symmetries are shghtly broken the 
conservation will be slightly imperfect with quarks of non-zero mass. The currents are 
then partially conserved. We take it that any dependence on the renormalization scale for 
and can be neglected. 

But there is a further consequence of the non-conservation of 7°^ which puts into 
question our interpretation of the expectation value of J^^ as the mean value of the total 
parton spin. To see this consider the protonic matrix element {P\J^^\P) in the Heisenberg 
picture. We may insert a sum over free parton states (assuming, as usual, completeness) 
so that 

{P\JI\P) =^\{kl...kr,\P)\^ {k,...k^\jl\h...kn) . (6.1.1) 

Strictly speaking the parton states should either be 'IN' or 'OUT' type states, but this is 
irrelevant for our discussion. 

Now in an interacting quantum field theory we are forced to renormalize the fields 
and operators, with the result that the value of matrix elements like {ki...kn\0\ki...kn) will, 
in general, depend on the renormalization scale Only if the operator O is conserved 
can one show that the matrix elements are independent of /i^. This is the case for and 
with massless quarks. Since anything of real physical content must be independent of 
/x^ we see that the expectation value of non-conserved operator cannot have any simple 
physical significance. Thus it is misleading to think of the expectation value of J^^ as 
"the physical spin of the parton" - it is not a fixed number, it depends on the value of /x^ 
and it can, in principle, have any value whatsoever. To emphasize this one should always 
indicate the renormalization scale i.e. write {ki...kn\J^^\ki...kn) ^■i- (Sometimes one refers 
to this as the spin at scale //^, but there is no implication that its value should be of order 

1-) 

Bearing this in mind we see that Lipkin's argument (Section 4.3.2) loses its force in 
the context of an interacting quantum field theory. The spin crisis is a crisis only in the 
domain of the Naive Parton Model, as was indeed claimed in the paper which introduced 



the concept of a "crisis" [LEA 88] . 

6.2 On the connection between QCD, the Qucirk Model and the Naive 
Parton Model 

Given that the expectation value of J^^ depends upon the renormahzation scale 
/x^ we would like to know if there is a value of /^^ at which the expectation value agrees 
with the Quark Model result i. e. corresponds to the physical spin of free quarks. We shall 
suggest that this happens as //^ — > 0, but because the Quark Model is an effective theory 
of the non-perturbative regime of QCD we do not have a rigorous argument. However we 
shall adduce some arguments from the perturbative domain to support our contention. 

The Quark Model is characterized by the fact that it does not contain gluons; the 
dynamics is in the quark degrees of freedom. Now neither the quark momentum operator 
nor the gluon momentum operator is conserved in QCD. Only their sum is. Thus the 
momentum fractions of a hadron carried by quarks or by gluons each depend upon the 
renormahzation scale. For large /i^, where ^^(/i^) is small so that perturbation theory can 
be trusted, we know that the momentum fraction carried by the gluons increases to the 
limit 16/25 at /i^ oo. Thus the gluons play a smaller role in the momentum sum rule 
at smaller /x^. Similarly, we shall see in Section 8 that the spin at scale carried by the 
gluons increases without limit as //^ — > oo, so that gluons play a lesser role in the angular 
momentum sum rule at lower /i^. 

Both these examples suggest that one is approaching the Quark Model as /x^ — > 0. 
Of course the above are perturbative arguments but it seems reasonable to assume that 
the trend continues down to non-perturbatively small values of /i^. At the other end of the 
scale we have the partonic picture which was invented to explain Bjorken scaling which 
holds as — > oo, which in the present context corresponds to //^ — > oo. And indeed, 
since q;s(//^) — > as //^ — > oo one does usually obtain the relationships of the Naive Parton 
Model in this limit, but in fact one never recovers Bjorken scaling and one has to utilize 
Q^-dependcnt parton distributions. 

The key exception to this rule is the gluonic contribution to the first moment of 
gi{x, Q^). Here, as will be seen in Section 6.3 and 7.2, we have a QCD correction propor- 
tional to as{Q'^) multiphed by the gluon spin at scale Q^, which increases hke In Q'^. The 
logarithmic decrease of ^^(Q^) as — > oo is just compensated by the increase in the 
gluon spin, leaving a finite, non-zero contribution as — > oo. This is quite anomalous 
and is, indeed, linked directly to the existence of the axial anomaly in QCD. It is thus a 
counterexample to the usual rule that one recovers the relationships of the Naive Parton 
Model in the limit ^ oo. 

This discovery catalyzed by the EMC result is of profound importance and, ironi- 
cally, will stand, even if it turns out that the EMC result is modified by future experiments. 

6.3 The cixial anomaly 

Consider the axial current 



made up of quark operators of definite fiavour /. (An implicit colour sum is always im- 
plied) . From the free Dirac equation of motion one finds that 



7^ 



5/i 
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(6.3.1) 



(6.3.2) 



where niq is the mass of the quark of flavour /. 

In the chiral limit niq — > (6.3.2) appears to imply that j/^ is conserved. If this were 
really true there would be a symmetry between left and right-handed quarks, leading to a 
parity degeneracy of the hadron spectrum e.g. there would exist two protons, of opposite 
parity. However, the formal argument from the free equations of motion is not reliable 
and, as shown originally by Adler, and by Bell and Jackiw [ABJ 69] (in the context of 
QED), there is an anomalous contribution arising from the triangle diagram shown in Fig. 
6.1. As a consequence the axial current is not conserved when niq — 0. One has instead, 
for the QCD case 

d'jL = ^ G::.^r = g Tr [g,^ g'"] (6.3.3) 



where is the dual field tensor 



1 



= 2wGr (6.3.4) 

and where a field vector or tensor without a colour label stands for a matrix. In this case 

G,.^(^y:.. (6.3.5) 

The result (6.3.3), which emerges from a calculation of the triangle diagram (Fig. 
6.1) using rUq — and the gluon virtuality /c^ ^ 0, is really a particular limit of a highly 
non-uniform function. If we take ruq^ 0, k'^ ^ the RHS of (6.3.3) is multiphed by 



2m7fc2 /Jl + 4m2/A;2 + l\ 

T{mi/e) = 1 - , In . (6.3.6) 

' 'l + 4ml/ k^ V^l + 4m2/A;2-i; 

We see that this anomaly corresponds to T — > 1 for {myk'^) — > 0. On the other hand, 
for on-shell gluons, A;2 — Q, and rUq ^ 0, i.e. in the limit {rrig/k^) — > oo the terms cancel, 
T ^ 0, and there is no anomaly. As explained in Section 8.1 the particular case A;2 = Q is 
irrelevant and there is no doubt that the anomaly is relevant. 

The anomaly induces a pointlike interaction between J^^ and gluons. That it is 
pointlike can be seen by taking different gluon momenta ki and /c2 in Fig. 6.1 and noting 
that the amplitude does not depend on the momentum transfer ki — k2 when niq — 
0. Therefore, in computing the matrix element of J^^ in a hadron state, we will get 
a contribution from the gluon components of the hadron as well as the more obvious 
contribution from quarks. From Adler's expression [ABJ 69] for the triangle diagram, 
modified to QCD, one finds for the forward gluonic matrix element of the flavour / 
current (our convention is £0123 = 1) 

{k,X\J^k,X) = '-^e,,,,k''e*^(X)e''(X)T(nil/e) 

= -^Sf,{k,X)T{ml/k') (6.3.7) 



27r 

where A is the gluon helicity and we may take 



Sf^{k,X)^Xk, (6.3.8) 



as the covariant spin vector for almost massless gluons. 



Using the methods of Section 3 we may then compute the gluonic contributions to 
the hadronic expectation value {P, S\J^^\P, S) . In this case the ghions being bound will 
be slightly off-shell i.e. k"^ ^ 0, but small. The full triangle contribution involves a sum 
over all quark flavours. We take rriu^md and to be <^ k"^ whereas mc,mh and rrit are 
^ k"^. The function T^mg/k"^) thus takes the values: 

T — 1 for u, d, s 
T = for c, 6, t 

and the gluon contribution is then given by [see (4.2.8)] 

^gw.(g2) ^ f'dxAg(x,Q') 
2tt jo 

or from (4.2.1) 

rtr'{Q') = -l^A9{Q'). (6.3.11) 

Ag{x^Q'^) is the difference between the number density of gluons with the same helicity 
as the nucleon and those with opposite helicity; its integral Ag(((5^) is the total helicity 
carried by the gluons. Note that if Nf massless quark flavours contribute in the anomalous 
triangle then on the RHS of (6.3.11) 

1/3 ^ Nf{e))/2 

where (e^) is the mean of the squared charges. 

Although (6.3.3) was derived perturbatively to order a^, it is believed to be an exact 
result. It was shown long ago by Adler [ABJ 69] that the anomaly is not influenced by 
higher order corrections at the 2-loop level. More recently Anselm and Johansen [ANS 
89a] showed that the 3-loop diagram of Fig. 6.2 effectively multiplies the anomaly result 
by a cut-off dependent constant. But this constant is the same as the one shown by Adler 
to renormalize Jg^ via Fig. 6.3. Consequently the QED version of (6.3.3) is unchanged. 
These results remain true in QCD for the matrix elements of (6.3.3). Further, it has been 
argued by Jackiw [JAC 85] that (6.3.3) is true even outside the perturbative realm of 
QCD. 

If we consider the anomaly contribution to the hadronic matrix element then it 
might appear possible to regard the contribution (6.3.11) as arising from a QCD correction 
to the quark distribution function. However this is quite incorrect since the important 
region of for the quark lines in the triangle loop integral, must be ^ oo in order 
to produce the pointlike behavior. The triangle really should be regarded as a point 
interaction between the current and the gluons. The result (6.3.10) is of fundamental 
importance. It tells us that the Naive Parton Model formula (4.3.4) for qq (and hence for 
rf in terms of the Aqf) is incorrect. We now have, instead, 

cy 

ao = AE - 3 (6.3.12) 
27r 

This result will be discussed in some detail in Section 7. But we note, immediately, 
that it has the fundamental implication that the small measured value of gq does not 
necessarily imply that AE is small. 



(6.3.9) 



(6.3.10) 



6.4 The cLxial gluon current K^^ and the gluon spin 

It is not difficult to sfiow tfiat if one introduces an axial gluon current 



3 

wfiere, as earlier, the matrix Ap — ^A^, then 



= s^'^p-IV I A, [cp, + \ g[Ap, A,] ) \ (6.4.1) 



d,K^ = ^ G^^Gr = Tr (G,. G"'') . (6.4.2) 



2 '^'^ 

Consequently, if = 0, the modified current 



Jl, ^ 4, -^K, (6.4.3) 



27r 

is conserved, d^Jl^ = 0. 

The matrix elements of the modified singlet axial current 



Jl^Jl-Nf^K, (6.4.4) 



are independent of the renormalization scale and should correspond with the value ob- 
tained in the Quark Model (no gluons; approximately SU{6) quark wave function) i.e. in 
the analogue of (4.2.8), 

{P,S\Jl^\P,S)^2M~aoS^, (6.4.5) 

we expect So independent of and thus oq — 1. 

We remarked in Section 1 that many of the operators corresponding to standard 
dynamical observables are in fact non gauge-invariant in a local gauge theory. In the gauge 
Aq{x) = the gluon spin operator becomes 

Sf = -s.jkAid'Al . (6.4.6) 

But in this gauge the cubic term vanishes for the spatial components of the vector 
and one finds 

Ki = -Sf (gauge A° = 0) . (6.4.7) 

Consider now the hadronic expectation value of K^. By the methods of Section 3, 
and in accord with (6.4.7) we find, in the gauge Aq — 0, 

{P,S\K^'\P,S) = -2MS^{P)Ag. (6.4.8) 

Let us now consider the question of the gauge dependence of this relation. The 
current K'^ is not gauge invariant. In an Abelian theory like QED there is no cubic term 
and the gauge transformation induced in as a consequence of 

A^{x)^A^{x) + d^Aix) (6.4.9) 

can be written in the form 

X^(x) ^ K^{x) - ^ [d''A{x)] e^,p,F^{x) (6.4.10) 



where we have used the fact that F^^ is gauge invariant in QED. Although K^{x) changes, 
its forward matrix elements, or expectation values do not, since the expectation value of 
F^"{x) vanishes. The latter follows because in QED F^" is related to the A^^x) entirely 
via derivatives, and one may utilize [P^^, f{x)] = —idf/dx^. 
In QCD, under 

^ UA^U-' + '-{d^U)U~' (6.4.11) 

one has 

G^, ^ UG^,U~' (6.4.12) 

and one ends up with 

^ K^+'^e^,^fsd''TT{A"U-'df'U) 

+ :^e,,^^Tr{U-\d''U)U-\d''U)U-\d^U)}. (6.4.13) 

The second term in the RHS of (6.4.13) is a total divergence and gives zero con- 
tribution to the expectation value of K^. The third term can also be shown to be a 
divergence [CRO 83] but in this case cannot be discarded because of the non-trivial topo- 
logical structure of QCD as discussed in Section 8.1. The last term may indeed be ignored 
for "small" gauge transformations i.e. those continuously connected to the unit trans- 
formation U — I, but it cannot be ignored for "large" (topologically non-trivial) gauge 
transformations. We shall return to study the consequence of this in Sections 8.1 and 8.4. 



7 Reinterpretation of the measurement of 

7.1 The role of the anomalous gluon contribution 

We showed in Section 4, that the measurement of g{ can be interpreted as effectively 
the measurement of oq which is proportional to the proton expectation value of the flavour 
singlet axial vector current Jg^. We also argued that the value found by the EMC i.e. 
Oo — was most unintuitive in the framework of the Naive Parton Model, since, therein, 
it implies AE ~ 0. 

On account of the anomaly, the connection between and AS is quite different 
from what it is in the Naive Parton Model, and we now have [see (6.3.12)] 

ao(g^) = AE-3^4^A5(g2) 

instead of (4.3.4). 

The anomaly has generated an effective point interaction between the virtual photon 
and the gluons and the small value of can now arise as a result of a cancellation between 
AS and the Q^-dependent gluon spin contribution (3Q;s((5^)/27r) ^g{Q'^). Why these terms 
should roughly cancel will be discusses presently. 

Quantitatively if we take ao — 0.17 at = 10 (GeV/c)^ where ag — 0.24 we can 
let the quarks carry, say 60% of the proton spin i.e. choose AS = 0.6 and (7.1.1) then 
implies 

Ag [g^ = 10 (GeV/cf] ~ 3.8 . (7.1.2) 



(7.1.1) 



This is, at first sight, quite unphysically large Q. But it should not be forgotten that 
according to the discussion in Section 6.1, this is the spin carried by the gluons at scale 
10 (GeV/c)^ and its value could be anything. To get some idea as to whether this value is 
absurd or not, we can use the evolution equation to be derived in Section 7.2 to estimate 
Ag{Q^) at some lower scale, closer to the Quark Model regime. Going as far as we dare, 
we take = AAqqj) where as — I and find that 

Ag{4AlcB)^0.7, (7.1.3) 

a not unreasonable value. 

In (7.1.1) we have shown no Q^-dependence in AS. This follows upon using (6.4.3, 
5 and 8) to get 

Zir Zn 

= AE (7.1.4) 

and the fact that ao, being related to the conserved current Jg^ is independent of the 
renormalization scale /i^ (or Q'^), as we shall explain below. 

Consider the axial charge Qs associated with the conserved current Jg^. We have 

Q, = J d^xJUx,t). (7.1.5) 

Let us ask how it would be possible for to depend upon the renormalization scale /x^, 
in a massless theory. The only way would be via the variable fit and such a dependence 
would then induce a dependence upon t. But we know that for a conserved local current 
the charge is time- independent. Thus must be independent of the renormalization 
scale. It follows, as claimed after Eq. (7.1.4), that oq or AE are independent of Q^. 

7.2 Why the gluon contribution survives as — > C50 

At first sight it seems surprising that the gluonic term in (7.1.1) survives at large 
Q^. It looks like an as correction which usually is expected to disappear as —>■ oo. 
However, this is misleading since the gluon spin at scale behaves just like [a;s((5^)]~^ 
as — > oo. This can be seen either from the spin-dependent Altarelli-Parisi evolution 
equations [ALT 77] or from the anomalous dimensions of the operators involved. We shall 
outline the latter approach. 

Now the axial current J^^ is multiplicatively renormalized which implies that the 
dependence on renormalization scale appears in a multiplicative factor. Consequently one 
can write 
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(7.2.1) 



where 7(0^) is the anomalous dimension of J^g. 

This is the standard form of an anomalous dimension evolution equation. What is 
a little unusual is that the power series expansion for 7(as) begins only in second order 
[KOD 80] I.e. 

7(«.)=72fel + ... (7.2.2) 



Note added in proof: A value of ao ~ 0.27 (see footnote in Section 4.4.1, [ADA 94]) would imply a 
smaller value Ag [Q^ = 10 (GeV/c)^] ~ 3.0. 



with 

72 = 16Nf . (7.2.3) 

The solution to (7.2.1) is then 
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X exp -i^ [a,(g2) - as{Ql)] (7.2.4) 



where, as usual, 

/3o = ll-^iV/. (7.2.5) 
Equation (7.2.4) implies that the quantity 



is independent of Q'^. Let us put it equal to a constant C, where C is the value of (7.2.6) 
at some arbitrary chosen value of Q^. 

We now take the hmit — > oo. Since AE is independent of and ^^(Q^) — > 
we obtain 

^^A^(Q^) = AE-C. (7.2.7) 
Thus it must be that Ag{Q'^) increases hke [^^(Q^)]"^ as increases. 

7.3 The anguleir momentum sum rule revisited 

We have suggested that the results on the first moment in polarized deep inelastic 
lepton-hadron scattering can be explained as a cancellation between a reasonably large 
quark spin contribution AE ~ (0.6 - 0.7), as expected intuitively, and the anomalous 
gluon contribution A^f. The above value of AE is consistent with various relativistic bag 
model calculations and, as will be explained in Section 9.3, with the so-called "Generalized 
Goldberger-Treiman relation" . 

However, in order to accomplish this cancellation, one requires a large gluon spin 
Ag at (Q^) ~ 10 (GeV/c)^ i.e. Ag ~ 4. As we have explained Ag grows indefinitely as 
increases according to the usual QCD evolution equations so that such a large value 
of A^f cannot be ruled out by our expectation that the gluon contribution is small at very 
small Q"^. 

It should be noted that the large Ag has important consequences for the total 
angular momentum sum rule. For a nucleon with helicity +1/2 moving along the z-axis 
we must have 

J. = 5,^ + 5f + L, = i (7.3.1) 

where represents the total orbital angular momentum of all partons. With ~ 0.7 
and independent of Q^, the growing value of has to be compensated by an analogous 
growth in the magnitude of L^. A detailed analysis of how this happens has been given 
by Ratcliffe [RAT 87]. 

In the process of evolution of a quark it radiates a gluon in some preferred helicity 
state and, in so doing, conserves its own helicity. It can then radiate again, with the 
gluon again produced in the same preferred helicity state. Thus gluons of some definite 
helicity are preferentially being radiated and Ag thus increases. But since the total angular 



momentum is conserved at each step of the evolution, something has to compensate for 
the component of spin carried by the radiated gluon. This, it turns out, is provided by the 
orbital angular momentum of the qg-pair in the process q ^ qg. Hence as Ag increases 
with each step in the evolution, so too does with opposite sign. 

The angular momentum sum rule (7.3.1) is thus completely in accord with the QCD 
growth of both and Lz as increases. 

8 The anomaly in the QCD Field Theoretic Model 

We would like to try to understand the previous results concerning the role of the 
anomaly in Ff in more intuitive terms, in particular in the framework of the Feynman 
diagram approach which provides a more rigorous, field-thcorctic approach than the Par- 
ton Model [EFR 88, ALT 88, LEA 88a, CAR 88]. (We refer to this picture as the QCD 
Field Theoretic Model). We also wish to study the gluon contribution to gi{x) itself. 

8.1 The factorization theorem 

According to the factorization theorem, which is the basis for deriving Parton Model 
type results from field-theory, and which has been proved to all orders in perturbative 
QCD, the cross-section asymmetry Aa-y*p 



for the absorption of a virtual photon on a proton of momentum P (the subscripts label 
helicities), is given by 



where the sum is over flavours and the partonic cross-sections refer to photo- absorption 
on quarks or gluons of momentum p = x'P. In the following we consider the 7*-proton 
collision either in the Lab frame with 7* moving along OZ or in an 'infinite momentum' 
frame where the fast proton and 7* move parallel to OZ. 

The partonic sub-process 7*5 is dominated by the Born diagram in Fig. 8.1 which 
physically corresponds to a forward jet with essentially zero kr (momentum perpendicular 
to the photon direction) in the Lab frame or in an 'infinite momentum frame' boosted 
from the Lab frame along the photon direction. (Of course there is another jet, in the 
backward direction in the frame where the 7* energy is zero, arising from the fragments 
of the proton). Integrated over kT the Born contribution leads to the usual expression for 
gi that one obtains in the Naive Parton Model. 

The partonic sub-process 'j*g is dominated by qq pair production as shown in 
Fig. 8.2. It is characterized by the formation of two jets with essentially opposite kr 
and is distinguishable from the process of Fig. 8.1 provided kx is large enough. If one 
calculates the differential cross-section asymmetry dAa^*g/dk'^, where fey and —kr are 
the jet momenta perpendicular to OZ, one finds an interesting dependence on the quark 
mass, as shown in Fig. 8.3, where the fcr-dependence is shown for three choices of quark 
mass: ~ Q^, 7^ but <C and = 0. For kr ^ mq it is negative and for 
kx <S rUg it is positive. For = it is thus always negative. Integrated over kx one 
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Vp = "^1,1/2 ~ 0'-l,l/2 



(8.1.1) 




(8.1.2) 



gets for the contribution to the scahng function of the partonic process £g — > £X, gj ^{x), 
from Fig. 8.2, for each quark flavour f [BAS 93] 

+ 

Here x stands for the Bjorken variable of the partonic process, x — Q^/2p ■ q. For — 
this simphfies to 

gT'ix, g2)U.=o = |g(2a; - l)[ln(gV/) - 2(lna; + 1)] . (8.1.4) 

We now wish to feed this into the analogue of (8.1.2) in order to obtain the gluonic 
contribution to gi{x). In doing so we must remember that in the QCD Field Theoretic 
Model the convolution formula (8.1.2) emerges as the lowest twist approximation to what 
was originally a Feynman diagram with an integration / d'^p. In this integration the point 
= is of zero measure. We saw in Section 6.3 that the anomaly is present if 7^ 0, 
but absent if = 0. From the present point of view the case p^ = is irrelevant and it 
is quite clear that the anomaly does contribute to DIS. 

It should also be remembered that in the QCD Improved Parton Model treating 
the above as a gluonic subprocess assumes that for gluons in the nucleon p^ is small i. e. 

p2 < 

The analogue of (8.1.2) for the contribution of Fig. 8.2 to the scahng function gi{x) 
is, for 3 massless flavours, 

-\f^- A^(^') gt'i^l^') ■ (8.1.5) 

O Jx X 

Taking the first moment of the convolution (8.1.5) we obtain for the contribution of Fig. 
8.2 to rl (for 3 flavours) 

r; = i 

^ 3 

and using (8.1.4) one flnds exactly the contribution (6.3.11) due to the anomaly. 

8.2 Study of the gluonic contribution to gi{x) and the definition of Ag(x) 

However, some caution is necessary about the interpretation of these results because 
they are based on perturbation theory, so may not be reliable as regards contributions 
coming from small kr- In the small kx region, with p^ limited, the momentum p' carried 
by the quark propagator in Fig. 8.2 is such that (p')^< O(Aqcd) so cannot be treated 
perturbatively. Moreover for (p')^ of this scale, the contribution should be considered as 
due to the interaction of the virtual photon with the quark-parton of momenta p'. This 
feature is a common one in the QCD Improved Parton Model. It has nothing specific to 
do with the polarized case. And the same diagram. Fig. 8.2, of course contributes in the 
unpolarized case, where the contribution of small kr is absorbed into a re-definition of 
the quark distribution function q{x), rendering it Q^-dependent i.e. one then deals with 
functions q{x, Q'^) which evolve as varies. 
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(8.1.3) 



dx ^g{x) 



dx g1 ^{x) 



(8.1.6) 



In the unpolarized case the situation is uncontroversial. One finds, for the Born + 
first order terms, a result whose structure symbohcally is of the form 

{1 + asHQ^p^) + asfq{x)} ® q{x) 

+ {asln{Qyp^) + asfg{x)}^g{x) (8.2.1) 

where q{x) and g{x) are unpolarized quark and gluon distribution functions. In LLA we 
neglect the terms asfq{x) and asfg{x), and, correct to order ag, factorize the remaining 
terms by introducing a factorization scale fi^: 

{l + a,ln(gVp2)}?(x) = {l + aMQ''/lJ^'') + OisHii''/p'')}q{x) 

^ {1 + a, ln(gV/i2)}{l + as ln(//Vp')} ?(^) ■ 

We also write 

as ln(gVp2) g{x) = {as HQ^fi^) + a, ln(/xVp2)| g(^^^ 
and then (8.2.1) becomes, to order as, 

{1 + as ln(QV/i2)} q{x, /x^) + a, \n{Qy jj^) g{x) 

where 

q{x, /i^) = [l + as ln(/xVp^)] q{x) + as \n{n^/p^) g{x) . 

In this way the contributions sensitive to low kr and singular in the limit — > gets 
absorbed into new, non-perturbative distribution functions which we cannot calculate, but 
have to measure. 

The leading logarithm In comes specifically from the ultraviolet region of the 
integration which kinematically runs up to 

4 < ^ 1-x _ 

Ax ^ ' 

It is therefore insensitive to different ways of factorizing or separating the low and high 
kx regions and there is no ambiguity about the definition of q{x,fi'^) in LLA. 

The situation is different when working beyond the LLA where one does not neglect 
terms of the form asf{x) in (8.2.1). They do not uniquely arise from the large kr region and 
it is essentially a matter of taste or convention whether they or part of them are included in 
q{x, yU^) or left as a piece of the gluon contribution. For this reason it is important, beyond 
the LLA, to indicate what convention is in use. Hence one uses q{x, /U^)j^, q{x, /i^)Dis etc.. 

In what way is the polarized case different? The key difference is the factor {2x — 1) 
in expression (8.1.4). For it means that the usual large logarithm ln((5^/p^) does not 
appear in the first moment, and the anomaly contribution comes from the second term 
2(lnx + 1) in (8.1.4). But this seems to be like the terms asf{x) in (8.2.1) and therefore 
inherently ambiguous. There are many papers on this question. For access to these see 
e.g. [HAS 91]. We believe that there is a reasonably persuasive argument in favour of a 
particular resolution of the ambiguity. Our approach is similar in spirit to that of Ross 
and Roberts [ROS 90] and Carlitz, Colfins, Mueller [CAR 88] (se also [MAN 90a]). 

To study the role of kx let us introduce a lower cut-off A in the kr integration i.e. 
in evaluating the diagram of Fig. 8.2 we keep kr > A. The result which replaces (8.1.3) 
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The first term in square brackets can be written 
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and we see that the entire term ln{(Q^/p^) [(1— and not just ln((5^/j5^) is insensitive 
to the value of A and therefore can be associated with the region of large kr- This is no 
surprise since the structure Q^{1 —x)/x just refiects the upper limit of the kx integration 
(8.2.2). It thus seems clear that the ln[(l — a:)/x] emerging from this term is a high kx 
effect and surely belongs to the gluon partonic sub-process, in agreement with the anomaly 
analysis. 

Returning to the case nig — 0, A = 0, we thus rewrite (8.1.4), for each flavour /, as 
2 27r 
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introduce the factorization scale //^: 



— Inxfl — x) — 2 



(8.2.5) 



and associate the first term with the gluonic contribution. The rest of the terms are then 
absorbed into the quark distribution A^/ + A^/. 

Finally then, choosing jjF' — Q^, the gluonic contribution from Fig. 8.2 to g{{x) is, 
via (8.1.5) 

9f(.)l....c = i I /; (2I - 1) ,n(i^) ) A,(.0 (8.2.6) 



and the modified quark contribution is 



g\{x^ Q^)|quark singlet — g AE(a;, Q^) 



(8.2.7) 



where the contribution of Fig. 8.2 yields 



AE(a;,g^) = AE(x) + ^ / ^ (2^ - l)[ln(gV/) 
Of course, the first moment 

AJ:^ [ dxJ:(x,Q^) 
Jo 

remains independent of in accordance with the discussion in Section 7.2. Clearly, the 
full definition of AE(x, Q'^) to LLA must also include the contribution from Fig. 8.4. 



(8.2.9) 



8.3 Measuring Ag{x) via 2-jet production 

It is important to note that the definition of Ag{x) adopted above associates it with 
the physical process of the electroproduction of two jets with large transverse momentum 
and —kx respectively. Measurement of this semi-inclusive process thus gives direct 
access to Ag{x). It is determined by a physical measurement and thus, manifestly, must 
be gauge invariant. (The formal question of gauge invariance is discussed in Section 9.3). 

We have for the contribution to gi{x, Q"^) coming from the two-jet events 

7*A^^jet(feT)+jet(-A;T)+X (8.3.1) 

9l-''\x,Q')\kr>K = \ t^9T'ix/x',Q')\kr>A A^(x') (8.3.2) 
where gj ^ was given in (8.2.3) and the lower cut off on kr is given by 

» mj, (8.3.3) 

where is the 'virtuality' of the gluon. 

8.4 The shape of gi{x) and the anomalous gluon contribution 

We have seen that the anomalous gluonic contribution to the first moment Ti per- 
mits an explanation of the EMC experiment in which the total spin carried by the quarks 
can be relatively large provided that A^f = dx Ag{x) is big enough. For example 
^quarks ^ q q jnucieon requires ~ 4 at n [GeV/cf. 

We consider now the x-dependence of gi{x) given by (we ignore perturbative QCD 
corrections) 

9l{x) = ^l^^[Au{x) + Au{x)] + ^[Ad{x) + Ad{x)] 
+ ^[As{x) + As{x)]j+gi{ 

•^) I gluonic 

(8.4.1) 

where 5'i(a;)giuonic is given in (8.2.6). 

We saw in Section 4 that without a gluonic contribution the EMC experiment forced 
us to have AE ~ and this was achieved by having a surprising large negative As 
contribution. 

It is interesting to contemplate the other extreme i.e. to take As{x) ~ 0. Assuming 
a roughly SU (3) invariant description of the sea quarks (this is not a good approximation 
for the unpolarized distributions) gi{x) will in this case depend only upon the valence 
quark distributions Auy{x) and Ady{x) and on Ag{x). The normalization of the valence 
quark distribution is fixed by the values of F and D (see Sections 4.2 and 4.3) and their 
shape is constrained by the requirement 

|AM^(a;)| < u^{x) \Ad^{x)\ < dy{x) . (8.4.2) 

Furthermore at a; ~ 0.3, where xg\{x) has its maximum value (see Fig. 4.2) there 
is essentially no contribution from Ag{x), so that there is a further constraint on A-u^,(x), 
Ady{x) in this region. It then turns out that one requires a significant negative contribution 
near x — 0.1 to fit the data in Fig. 4.2. The question then is: can a large gluonic term 
provide such a contribution? It is found empirically [ROS 90] that this is not possible as 



is demonstrated in Fig. 8.5. One can trace this failure to the feature that the factor in 
parenthesis in (8.2.6) has a double zero at x' = 2x. Given that Ag{x) is bounded by 



\Ag{x)\<g{x) 



(8.4.3) 



and that g{x) dies out rapidly as x increases, the region of overlap in the convolution 
in (8.2.6) where both the kernel and Ag{x) are big is a very small one and in effect the 
integration only samples Ag{x') for a small portion of the region x < x' < 1. 

Of course there is no reason to suppose a total absence of a strange quark contri- 
bution. The problem with the EMC analysis without the anomalous gluon contribution 
was that it demanded a surprisingly large negative As{x). Figure 8.6 taken from [ROS 
90] shows that a reasonable compromise is possible. The fit shown corresponds, for F/D 
= 0.55 (which is approximately the value used in Section 4), to having 
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AS 



0.84 
0.38 



Ad 
Ag 



-0.41 
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-0.05 



3. 



(8.4.4) 



The exact values should not be taken too seriously since they depend upon the sim- 
ple, but reasonable parametrization used for the polarized quark distributions in [ROS 90], 
and because the small x behavior used does conform to the expectations from Regge the- 
ory. However the fit does indicate that the EMC data can be explained with a reasonable 
strange quark contribution, a sensible AS and a moderate Ag ^ 



8.5 The anomaly in i/p-elastic scattering and weak interaction structure 
functions 

We mentioned in Section 6.3 that there is no anomaly in the limit m'^/k'^ 
Analogously, we see that in the limit of Jp^ oo (8.1.3) becomes 
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and the first moment vanishes. In this case it is natural to consider the contribution from 
Fig. 8.2 as a modification of the quark distribution functions. 

There is an interesting application of the role of quark mass in the anomaly in 
elastic neutrino proton scattering, which is mediated by Zq exchange as shown in Fig. 8.7. 
Since the parity- violating Zo-quark coupling involves 7'^75 Fig. 8.7 also depends upon the 
proton matrix element of the axial vector current. In the Naive Parton Model, for forward 
scattering, the contribution is, in principle, proportional to 



5q = (Am + Au- Ad- Ad) + (Ac + Ac - As - As) 



15.2) 



(the sign change for the Is = —1/2 quarks, compared to the electromagnetic case, is due 
to the Zq couplings). 

Experimentally [AHR 87] it seems that 5q is about 12% larger than expected from 
the Am, Ad contribution (see also the discussion in Section 4.2.3). This has been inter- 
preted [ELL 88] as evidence for a large negative As in agreement with (4.3.17), assuming 
a negligible Ac. However, the anomaly also contributes to this process via the diagram of 
Fig. 8.8. 

Note added in proof: Presumably with the new proton data [ADA 94] an even more reasonable AS 
and At; would be possible. 



By the above discussion there is no contribution from the heavy quarks, but for 
each massless quark, in the QCD Field Theoretic Model, 

Aq + Aq ^ Aq + Aq - ^ Ag 

ZTT 

Hence, in (8.5.2) 

Q/Q 

5q^5q+—Ag. (8.5.4) 

In the extreme case As = one then requires a rather large Ag ~ 5 to fit the data. 
As discussed in Section 8.4 such a scenario leads to difficulties with the x-dependence of 
gi{x). But the more reasonable choice As = —0.05, (? ~ 3 is perfectly acceptable. 

The weak interaction structure functions discussed in Section 3.5 also receive an 
anomalous gluonic contribution via a diagram analogous to that of Fig. 8.2, with the 
photon replaced by a Zq (neutral current) or a by (charged current). Such contribution 
leads again, for each massless quark, to the result (8.5.3) [VOG 91]; for the structure 
functions measured in charged current DIS, which involve diagrams with transitions from 
light to heavy quarks [see Eqs. (3.5.8)], there are interesting quark mass effects [VOG 91; 
MAT 92, 92b]. 

8.6 Operator Product Expansion vs. QCD Improved Parton Model 

In Section 5 we stressed that only the fully inclusive electroproduction process could 
be expressed in terms of the matrix element of a product of electromagnetic currents 
(5.1.1) and hence only in this case could one utihze the Operator Product Expansion. 
Our definition of Ag{x) has focused on the semi-inclusive large kr 2-jet process and it 
might thus seem strange that our results for the first moment should agree exactly with 
the anomaly result obtained using the Operator Product Expansion. But this is simply 
a 'red herring'. We have not thrown away the part of the cross-section corresponding to 
small kr- It is simply absorbed into the redefinition of the quark distribution functions. 
So the results of the Operator Product Expansion should agree with the QCD Improved 
Parton Model results. 

Moreover, if we return to the triangle anomaly discussed in Section 6.3 we may 
ask what region of the kr integration involved in evaluating the triangle diagram was 
important for the anomaly. As stressed in [CAR 88] it is indeed the region of large kr 
that is relevant, so there is consistency between the two methods of deriving the anomaly 
contribution to the first moment of gi{x). 

These considerations, as emphasized in [ROS 90], draw our attention to an impor- 
tant issue, namely that whatever cut-ofi' or factorization scheme we use, we must use it 
consistently if we wish to compare results from the Operator Product Expansion with 
those of the QCD Improved Parton Model. We must also be careful when dealing with 
physical quantities that depend on both polarized and unpolarized distribution functions. 

For example, given that g±{x) measures the number density of gluons polarized 
along or opposite to the spin direction of the proton, they must be positive. Thus one 
expects 

\Ag{x)\ = \g4x) - g4x)\ < {g^x) + g4x)\ = g{x) . (8.6.1) 

But g{x) is measured in unpolarized scattering and depends upon the factorization con- 
vention utilized there. Clearly (8.6.1) and analogous relations for quark distributions can 
only be expected to hold if the same cut-off convention is used for Ag{x), g{x), Aq{x) 
and q{x). The relevance to the phenomenology of gi{x) was discussed in Section 8.4. 



(8.5.3) 



Finally, there is another 'red herring' in the literature which needs to be pointed 
out. The anomaly emerges from the triangle diagram only if kx is integrated up to infinity, 
as it must be in the Feynman diagram. It is claimed that, on the contrary, the anomaly 
emerges when integrating up to ~ in the QCD Improved Parton Model. The latter 
statement is misleading. True, at fixed x, /cf,™^ fti Q"^ but to get the anomaly we must 
integrate over x from to 1. The region x — > corresponds, via (8.2.2), precisely to 
k^^ oo. 

In summary there is perfect agreement between the two approaches to the gluon 
contribution in polarized deep inelastic scattering, provided care is taken to use consistent 
factorization or renormalization conventions in both. 



9 Non-perturbative efTects in the interpretation of the measurement of 

rf 

The previous discussion has focused on the role of perturbative aspects of the axial 
anomaly in the theoretical description of polarized DIS. But the axial anomaly also has 
deep and non-trivial non-perturbative consequences in QCD because of the complicated 
topological structure of the non-abelian theory. We first provide a brief survey of the 
topological issues and than turn to their role in DIS. We show that they support the 
conclusions drawn from the reinterpretation of the measurement of including the effect 
of the anomaly. Finally we give a brief survey of the open questions in the understanding 
of the EMC result. 



9.1 Topological effects and the cixial ghost in QCD - a pedagogical reminder 

We follow the approach of Dyakonov and Eides [DYA 81]. A more detailed ped- 
agogical treatment can also be found in [LEA 94]. Consider a purely gluonic version of 
QCD, and choose for simplicity the gauge in which Aq{x) = 0. While satisfying this gauge 
condition the theory is still invariant under arbitrary time-independent gauge tranforma- 
tions 

Aj UAjU-^ + ^ {djU)U-^ (9.1.1) 

where, as usual, Aj without a colour label stands for the matrix {\°'/2)A'j and U = U{x) 
only and where, for a consistent theory, one utilises transformation such that U{x) I 
as |a;| — >• oo. 

The fields Aj play the role of generahzed coordinates in the quantum theory and 
the role of generalized momenta is played by 

djC 

<{x) - ^ = -Ai = E^j (9.1.2) 

where -E" is the colour-electric field. 

The generalized coordinates and momenta satisfy the usual canonical equal time 
commutation relations 

[<(x,i),4(y,i)] = -5ij5ab5\x - y) . (9.1.3) 
The Hamiltonian, quite analogous to QED, turns out to be 



d^x Tr[£;2 + B^] 



(9.1.4) 



where the components of the colour-magnetic field are defined by 



B';^\e,MG'J^. (9.1.5) 
Unlike QED, the fields do change under an arbitrary gauge transformation: 

UEjU-^ Bj UBjU-^ . (9.1.6) 

Bearing in mind that E'^ = A we see that the first term in (9.1.4) plays the role 
of the kinetic energy whereas V = Tr[S^] represents the potential energy. 

The analogue of a wave function in ordinary quantum mechanics is here a functional 
^'[A(x)] of the fields and the stationary states and energies are found from 

H^ = E^. (9.1.7) 

Now the gauge transformations U{x) fall into discrete classes labelled by an integer 
n, the winding number. The class with n = are called small gauge tranformations and can 
be built up from the identity by a series of infinitesimal transformations - by smoothly 
varying the parameters which specify them they can be continously deformed into the 
identity transformation. The classes with n^Q are called large gauge transformations and 
cannot be continuously transformed into the identity transformation (for an elementary 
discussion sec [LEA 94]). 

In Section 6.4 wc introduced the gluon axial current [see (6.4.1)] which is a 
function of the fields A'^. Based on this we can construct a new variable X, a function of 
the fields yl", (often called a collective coordinate): 

2 

X = ^ld'x Koix) (9.1.8) 

and such that it changes by the integer n when the fields undergo a gauge transformation 
of winding number n but does not change under small gauge transformation. In fact 

X^X + ^ Jd^x e,,kTT[{d,U)U-\djU)U-\d,U)U-'] (9.1.9) 

and the change in X can be shown to be just n. Thus 

X^X + n. (9.1.10) 

Note that the non-vanishing of the trace term (9.1.9) is a consequence of the non- 
Abelian nature of QCD. For normal QED, {i.e. in 4-dimensional Minkowski space) on the 
contrary, U = e^^^^^ and the trace term vanishes. 

Now we can imagine expressing the potential energy V in terms of X and any other 
necessary variables. If then we plot V vs. X it must have the same value at all values 
Xn — X + n, since it is, via (9.1.6), invariant under all gauge transformations, large or 
small. In other words if we express the potential in terms of X and the other required 
variables, V = V{X, Y, ...), then V is a periodic function of X with period 1 as shown in 
Fig. 9.1 This periodicity was first noticed by Faddeev [FAD 76], Jackiw and Rebbi [J AC 
76], 't Hooft ['THO 76] and CaUan, Dashen and Gross [CAL 76]. 

Before exploring the consequences it may help to examine a very simple example 
of the above phenomenon. 



Consider electrodynamics in one time and one space dimension, where A" = {Ao, 
and define Ka = SapAf^ , (egi = —^lo = 1) so that Kq = A^. 
Let 

X = ^ dx Ko{x) = / dx A^{x) . 

ZTT J-oo ZTT J-oo 

Then, under the gauge transformation U = e'^^^^\ 

g ^dx ^ g dx 

Consequently 

But the requirement U{x) — > 1 at spatial infinity means that (pioo) and (f){—oo) must both 
be multiples of 27r, but not necessarily the same multiples. Hence X — > X+ (an integer 
n). 

Returning to QCD we see that we are dealing with a quantum problem with periodic 
potential i.e. an analogue of an electron in a crystal lattice. The energy spectrum will thus 
be a band spectrum. Now, near the bottom of the lowest energy band the electron behaves 
like a free particle with E = k"^ /2m* and m* an effective mass which depends on the details 
of the shape of the potential. As a consequence the Green's function (0|T[x(t)x(0)]|0) has 
the property that 

/ dt e^-* (0|T[x(t)x(0)]|0)|_o — ^ (9.1.11) 

like the propagator of a free particle |^ The electron can move freely in the lattice. Note 
that these results are inherently non-perturbative. 

Analogously in QCD this means that the vacuum state wave function is not localized 
in X at one of the minima X = integer, but that the ground state is spreat out in X in 
the form of a Bloch wave 

l^)=E e^"V;...) (9.1.12) 

n 

where the state \n] ..) is approximately localized at the point X = n, and the value of 
6 which specifies the lowest energy states is beyond our power of calculation. The above 
phenomenon occurs because the barriers between the different minima of the potential 
can be penetrated. In QCD the instantons, which are solutions of the Euclidean version 
of the theory, with finite actions, allow one to construct a semi-classical approximation 
to \9), but the structure (9.1.12) is a general result, independent of any consideration of 
instantons. In the following the vacuum state |^) is to be understood as the true QCD 
vacuum. 

Analagously to (9.1.11) the Green's function involving X{t)X{0) will have a pole at 
uj = 0. Bearing in mind how X is constructed from Kq, we expect to find a massless pole 



This is an immediate consequence of the equation for the Green's function of a free particle 



{0\T[x{t)x{0m=z5{t). 



in the momentum space Feynman propagator for K^{x)^ which is the Fourier transform 
of (0|T[ir^(x)iri,(0)]|0). To avoid a plethora of constants, let us introduce 




27r 



(9.1.13) 



where Nf is the number of light flavours - 3 in our case. We thus expect to have for the 
Kfj_ propagator: 



where T and (3 are constants. 

Some comments are necessary about the above expression. The above form will hold 
in any covariant gauge but the value of j3 will be gauge-dependent. The fact that d^j_K^[x) 
is gauge-invariant implies that the longitudinal part of the propagator (9.1.14), i.e. the 
part proportional to q^qu/(f, must be gauge- invariant. Thus the value of T cannot depend 
upon the gauge and any gauge-dependence associated with the propagator can only arise 
from an explicit dependence Tip on f3. 

The massless ghost state thus introduced will be referred to as the 'axial ghost'. The 
'ghost state' exists in an unphysical sector of the Hilbcrt space. (Recall that in quantizing 
a gauge theory the physical states form only a subspace of Hilbert space). Finally note 
that the coupling of the ghost will be such that the apparent 1/g^ double pole never 
actually occurs and, indeed, such that no physical matrix element has a pole at — 0. 

Notice that aside from the normalization factor, the LHS of (9.1.14) corresponds 
exactly to the momentum space Feynman propagator for the 'field' K^lx). We shall in- 
troduce the shorthand notation 



for these momentum space propagators. 

Equation (9.1.14) implies that in QCD there exists a zero mass axial vector ghost 
particle which couples to the axial gluon current K^. The necessity for the ghost pole in 
QCD was recognised as early as 1975 by Kogut and Suskind [KOG 75]. The ghost pole 
was invoked by Veneziano [VEN 79] to resolve the U{1) problem. 

9.2 A reminder about the U{1) problem 

Although we shall not discuss the resolution of the problem we give a brief reminder 
as to what exactly the problem is. 

As mentioned in Section 6.3 the QCD Lagrangian in the chiral limit rriq seems 
to contain an additional symmetry, namely that it appears to be invariant under separate 
flavour transformations on the left and right-handed parts of the quark flelds. Such a 
symmetry would lead to a parity doubling of the hadron spectrum if the vacuum were 
invariant under chiral transformations and therefore unique. For example, it would imply 
a (J)^ = (1/2)^, I = 1/2 partner for the nucleon. Experimentally the nearest resonance 
with these quantum numbers to the nucleon is the A(1535), some 600 MeV heavier than 
the nucleon. Now the quark mass parameters in the Lagrangian are believed to be small 
(m„, nid < 10 MeV/c^ and nis ^ 150 MeV/c^) on the scale of hadron masses, so the 
chiral limit should be a reasonable approximation. It is then difficult to believe that the 



I d^x e"^-^0\T[K^{x)K,{0)]\0) = -iT 



■4 (gt^i^ - (li^(lu/q^)(3 + q^qi^/q^ 
q^ -\- ie 



(9.1.14) 




(9.1.15) 



non-zero quark masses are responsible for the mass gap between the nucleon and the 
A^(1535). 

It is thus assumed that the chiral symmetry is broken principahy as a resuh of 
a spontaneous breaking mechanism which gives rise, for example, to a quark-antiquark 
condensate. There is no longer a unique vacuum, but a continuum of states related to each 
other by chiral rotations and the symmetry is broken by choosing one of these as the true 
vacuum state. As a consequence there should exist a zero mass pseudoscalar Goldstone 
boson. Its mass would be strictly zero in the limit niq = and would be expected to be 
non-zero, but small, in the realistic case of small quark masses. There is one Goldstone 
boson for each broken generator of the original symmetry. 

Since the QCD Lagrangian is, in fact, invariant under U{3) transformations of 
flavour there ought to exist nine light pseudoscalars whose masses would be zero if the 
were all zero. The obvious pseudoscalar nonet is (tt, K, rj, rj'). In the following we shall 
treat the rj as the eigth member of a pure SU (3) octet, and the 77' as a pure SU (3) singlet. 
(In reality there is a small amount of mixing which we ignore since it has no bearing on 
the arguments to follow). 

It turns out, however, that the mass of the ?7'(958) is much larger than expected 
on theoretical grounds. The other 8 pseudoscalars, associated with the breaking of chiral 
SU{3)f, have acceptable masses, so it is only the Goldstone boson associated with U{1)f 
that is in difficulty. It is this that is referred to as the problem'. 

9.3 The axial ghost 

It turns out that simply introducing an axial ghost pole via (9.1.14) leads to a result 
for the 7]' mass which is inconsistent in the chiral limit. The correct procedure seems to 
be the following: 

i) The above ghost is considered as a 'bare' ghost. 

ii) One introduces a quark-antiquark pseudoscalor flavour singlet field rjoix) correspond- 
ing to the 'particle' 

It^o) = — /= \uu + dd + ss) (9.3.1) 
v3 

which would have corresponded to the 77' if there were no axial anomaly. The 770 has 
a mass mo similar to the other pseudoscalor mesons. 

iii) We introduce a direct coupling between the ghost field and r]o(x) which results in 
mixing between the bare ghost and the rjo and then study the consequences of the 
interaction. 

To this end we introduce the bare massless ghost field G^^{x) which couples di- 
rectly to so as to induce the mass zero pole in the propagator (9.1.14). The 
non-perturbative transition amplitude for the - coupling will be indicated graphi- 
cally as in Fig. 9.2. The bare ghost propagator is then given by (9.1.14) without the factor 
T"^. The coupling of the ghost to the 770 will be described by an effective Lagrangian term 
—Am G^{x)d'^r]Q{x) where Am is a coupling constant that will play the role of a mass 
correction. This induces the diagrammatic vertices shown in Fig. 9.3. 

The bare and 770 now mix and the mass eigenstates have to be found by diag- 
onalizing their propagators. The simplest method is by summing the diagrams. Fig. 9.4, 
generated by the interaction. One finds for the masses of the renormalized ghost and 
the renormalized physical 77' 



ma — m^, — ml + (Am)^ . 



(9.3.2) 



The extra term (Am)^ can be adjusted to raise m^/ well above the typical pseudoscalar 
mass mo and this resolves the U(l) problem. 

We can estimate the size of the parameter T involved in the coupling of and G° 
as follows. 

Consider the divergence of the flavour singlet axial current. We have 

d^Jl = Pix) + d^K,{x) (9.3.3) 

where 

P(a:) = 2z^m/g/75g/. (9.3.4) 
/ 

The current Jg^ — A'^ is divergenceless in the chiral limit and is thus analogous to 
the flavour non-singlet currents. We thus define /o by 

{0\J',^-K,\r)o)^iVSfoq, (9.3.5) 

and expect to have /o ~ A = 135 MeV. Now one can show that the parameter T is 
related to /o and Am via 

T2 = x/3/oAm. (9.3.6) 

The value of mo can be obtained from the magnitude of the quark condensates (qq) 
known from QCD sum rules and this leads, via (9.3.2) and (9.3.6) to the estimate T ~ 
460 MeV [DYA 81]. 

For the longitudinal part of the complete propagator for G^, which we shall need 
later, one finds from the diagrams in Fig. 9.4a 

{T[G.G,]),= ^('-^] 4^ (9.3.7) 
q^ \ q^ J q^ — m^, 

and for the non diagonal — > rj' propagator given in Fig. 9.4b 

= (9.3.8) 

q q in^i 

We see, therefore, that the axial ghost plays an essential role in resolving the U{1) 
problem in QCD. A more sophisticated treatment, allowing for some breaking of SU{3)f 
does not affect the above results significantly, but does lead to a successful calculation of 
the mixing between the singlet \rjo) and the eighth member of the octet Irjs)- 

Finally we note that there is no danger of the production of the axial ghost as a 
physical particle. The cross-section for such production would have the form 



daa <x M^Mle''{X)[e''{X)Y6{k'')d''k 



where Ai^ is the Feynman amplitude and e^{\) the polarization vector for a ghost of 
helicity A and momentum k. As long as the axial ghost coupling to physical particles is 
of the derivative type we will always have JA^ oc k^. The cross-section will thus involve 
k^k'^ e^{X)[e^{\)]* 6{k'^). But for a physical particle k^e^{X) — and the cross-section 
vanishes. 



9.4 Role of the ghost in the nucleon spin problem 

We showed in Section 6.4 that the nucleon expectation value of K^^ plays a direct 
role in polarized deep inelastic scattering. In Eq. (6.4.8) we calculated periurbatively the 
value of this matrix element. Now, because of the coupling of to the non-perturbative 
axial ghost, there appears to be the possibility of, in addition to the RHS of (6.4.8), a 
purely non-perturbativc contribution. Similar reasoning applies to matrix elements of ,7°^^ 
[see (6.4.4)] which, perturbatively, was interpreted as the contribution of the quark spin. 
We shall now see how this leads to a new relation for AE. 

In any covariant gauge the nucleonic matrix elements of must have the following 
general form: 



{P',X'\K,mP,X) = Uy{P') jMq') + q,IC2{q') 



75Ma(P) (9.4.1) 



where ^1,2 (?^) are scalar form factors, q^ — P'^ — P^ and A, A' are hehcities. 

We shall now show that the ghost contributes only to /C2(?^)- The relevant Feynman 
diagrams are shown in Fig. 9.5 where the propagators are the complete propagators shown 
in Figs. 9.4a, b. 

For the interaction of the r^o with the nucleon we take a standard pseudoscalar 
coupling term in the Lagrangian: —igrj^NN '4^n^5'4^nVo and for the coupling of a term 
—igooNN G^°^''('0iv75'0iv), chosen so that ghosts cannot be emitted as physical particles. 
These give rise to the Feynman vertices shown in Fig. 9.6. 

The diagrams of Fig. 9.5 yield a contribution to (P', X'\Kn{0)\P, A) of 

^ ^^5-" ^ (^Amg^^Nff - ml c/coivjv) • (9.4.2) 
Hence, comparing with (9.4.1) the ghost contributes only to 

lCt"'\q^) = — — 2" i^rng - mlgG^^N) , (9.4.3) 

i.e. the ghost induces a pole in lC2{q^) at q^ = 0. Since no other contribution can do this, 

(9.4.3) becomes an exact expression for K,2{q^) as q^ — >• 0. 

Note that the gauge-dependent parameter {3 [see (9.1.14)] in the ghost propagator 
does not appear in (9.4.3). Thus the value of K,2{q^) as q'^ — > 0, as given by (9.4.3), must 
be gauge-invariant. 

Now the gluonic contribution that was calculated pertubatively in Section 6 occurs 
only in /Ci(g^). From (6.4.8) one sees that 

^1(0) = -^ A^. (9.4.4) 

We are now able to give the formal argument for the gauge- invariance of A^f. We 
know that d^K^ is gauge-invariant. But its forward matrix element just involves a linear 
combination of /Ci(0) and [q^K^2{q^)]q'^=o i-^. a linear combination of Ag and other physical 
parameters [see (9.4.3)] which do not depend upon the gauge. Thus Ag itself must be 
independent of the gauge. 

Clearly the ghost will also contribute to the modified singlet current J^^ defined in 

(6.4.4) . If we put 

(P', AV°^|P, A) = Uy{P') kGi(g^) + q,G2{q')h,u^(P) (9.4.5) 



then we will have, for the ghost contribution 

and for the perturbative QCD contribution from (6.4.5) and (7.1.4) 

Gi(0) = AE. 



(9.4.6) 



(9.4.7) 



We shall now derive a connection between AE, the r]' decay constant and the rj'NN 
coupling. 

Since in the Parton Model we work with massless quarks, consider the chiral limit, 
• 0. Then in (9.4.3), m^, {Amf, ml ^ and using (9.3.6) we have for (9.4.6) 



m 



G2{q^) = — ^ /o QnoNN + non singular terms . 

Now in this limit Jg^ is conserved. Hence 

= {P',X'\d^Jl\P,X) = iq^{P',X'\J^P,X) . 
Using (9.4.5) we have 

= Uy{P'h,u^{P)[2MG^{q') + q'G2{q') . 

Hence 



2M 



G2{q') 



and via (9.4.8) 



Giiq') = ^^^^ + 0{q') . 



Going to g2 = and using (9.4.7) we obtain [YEN 89, EFR 90] 



AE 



QrioNN 



2M 



(9.4.8) 

(9.4.9) 

(9.4.10) 
(9.4.11) 

(9.4.12) 
(9.4.13) 



known as a Generalized Golberger-Treiman relation (TGV in [VEN 89]). 

Note that the form of (9.4.13) differs from that in [VEN 89]. In [VEN 89] one starts 
with the matrix element of 7°^ expressed hke (9.4.5) in terms of scalar functions Gi^2{q'^)'- 



{P',X'\Jl\P,X)^UyiP') 



l^.G^{q^) + q^G2{q^) 



l5Ux{P) 



(9.4.14) 



[note that Gi(0) is then identical to gq used in Section 4.2 and Section 6] and one finds 



Gi(0) 



2M 



fv'idri'NN ~ ^QNn) ■ 



(9.4.15) 



In this form Fq^^ is the total contribution from gluonic effects, both from the 
ghost and from the gluon-parton, and of course G*i(0) involves not just AS but AE — 
(3«s/27r) A^f. A rcnormalization group analysis then leads [SHO 92] to the identification of 
AE with the first term on the RHS of (9.4.15). However, it seems to us that there is some 
point in separating the ghost and gluon-parton contributions. The former is totally non- 
perturbative and non-Abelian and is absent in QED. The latter is partially perturbative 



and would be present also in an Abelian theory. Further, the two contributions have a 
significantly different kinematical structure in the matrix element of K^^ and so can, at 
least in principle, be measured separately. Obviously one learns about Ag in deep inelastic 
scattering; but also in various high processes, as discussed in Section 11. And it may 
be possible to study ghost effects in NN spin dependent processes in the Regge region, 
as explained below. 

The result (9.4.13) is unchanged if we allow for a non-zero, but flavour symmetric 
quark mass. Moreover it can be generalized [EFR 91] to allow for both SU{3)f breaking 
and isospin breaking. The change is small provided that allowance is made both for 
breaking of the equality of the quark masses and for - 1] - 1]' mixing. If the latter 
is neglected one obtains large and misleading corrections, even the wrong sign for the 
gluonic contribution [CHE 89]! There is a significant cancellation between the effects of 
the two types of symmetry breaking, as is discussed in detail in [EFR 91]. 

Now (9.4.13) is interesting because it suggests, in principle, and independent mea- 
surement of AS. The coupling of 1]' to nucleons can be obtained from studying the role 
of ?7'-exchange in nucleon-nucleon elestic scattering. A phase shift analysis in the frame- 
work of a one boson exchange potential model, leads to g^iNN — 7-3 with, presumably, 
rather large uncertainty [DUM 83] . But in this analysis there was no ghost exchange. So 
the measured coupling does not really correspond to the Qr^NN ™ (9.4.13). The ghost 
would generate a contact NN potential and for small \t\ -C m^,, what is measured is 
actually ^Jg^^J^f^ — iri^' Qa^ff and not Qj^^nn, but, as we shall argue, it seems hkely that 

l"^r?' 9gNn\ ^ 9-qoNN- 

Now Qrj'NN can also be estimated from r]' 27 decay. The decay rate has been 
calculated using a triangle diagram made up of the baryon octet, and assuming, as happens 
for 7r° 27, that this gives the same result as a calculation based on quark triangles [BAG 
90]. This yields ^^,jviv = 6.3 ± 0.4. 

Both the quoted values arc close to the SU{6) value [TOR 84], fif^/TVJV — 6.5. We see 
that there will be consistency between all these values provided that m^/ QgoNN is small 
enough. 

It may be helpful to note that there is an analogous situation suggesting a small 
ghost coupling. In the calculation of the 2-photon matrix element of Jg^, (7(^1)7(^2)! -^5^ |0) 
the conservation is broken by a contribution from the QED anomaly i.e. there is an 
additional term (a — 1/137): 

d'Jl = ■ ■ ■ + ^A^c E 4 W^"^^"'^ (9-4-16) 

where A^^ = number of colours, Cf = quark charge in units of e and F^^ is the electro- 
magnetic field tensor. Of course, in our previous strong interaction considerations, we 
neglected the RHS of (9.4.16) altogether. Here, because we are dealing with photons, it 
is this term, and not the QCD anomaly, that is relevant. Now the above matrix element 
of Jg^ is proportional to g^^i^-y (directly related to the rate r]' — > 27) and earlier studies 
[ROE 90] showed consistency with (9.4.16) at the 10% level. But recently it was realised 
[VEN 92] that the ghost should also contribute. This has the effect of replacing (jf^'-^-y by 
gn'-yy—rrir^i gc-yy and the earlier consistency of (9.4.16) implies that m^/ ga-y-y must be < 5% 

of grj'^-y. 

The smallness of the ghost couphng to quarks in both the A^A^ and 77 case is linked 
to the success of the OZI approximation. It would also be expected to be small in the 



1/Nc expansion [VEN 92], where for any state X, 



rrirf Qgxx 



gri'XX VNc 

but this does not explain why it is as small as it seems to be. 
In summary it thus seems not unreasonable to assume 

^??' 9GoNN 9rioNN 

and thus 

Qr^oNN'^^ - 7. (9.4.17) 

Returning to (9.4.13) we take fo^if^^ 135 MeV and find 

AS -0.8. (9.4.18) 

This supports the interpretation of the EMC result as a cancellation between a 
significant AS and a large gluon contribution as discussed in Section 7.1. It is in agreement 
with the naive expectation that the quark spin contribution is ~ 1/2. At low this is 
the dominant contribution. At high gluon spin and orbital angular momentum far 
outstrip AS. But this should not be a cause for concern. As stressed in Section 6.1 the 
proton probed at large is very different from the static proton. 

It is clearly important to be able to test the validity of the above analysis by 
trying to measure separately g^^^j^fq and QgoNN- O'^g proposal [EFR 91a] is to look for 
the spin-fiip effects caused by the r]' and ghost in a situation where other ffip mechanisms 
are absent or suppressed e.g. in Aatat or D^n in proton-deuteron or deuteron-deuteron 
scattering, in the Regge region \t\ <S s. The I — 1 exchanges of p and tt are forbidden 
and the iP, iP', u and / exchanges are predominantly non-fiip. So 77, ry' and the ghost 
should be the most important exchanges. The ghost exchange presumally contributes a 
fixed pole in the complex-J plane where the 7] and rj' arc associated with moving Regge 
poles. Consequently the energy variation of their contributions will be quite different in 
character and this should become discernable at high energies. 

9.5 Attempts to provide a physical interpretation of the EMC result 

The result of the EMC experiment can be summarized as the statement that [see 
(6.3.12)] 

|ao| = AS-3-^At/ < 0. (9.5.1) 



Now from (6.4.2) we have 

{P', X'ldi'KJP, A) = -3— (P', X'lE"- B^IP, A) (9.5.2) 

TT 

where E"" and are the colour-electric and magnetic fields, and from (9.4.11) the LHS 



IS 



[2M/Ci(g2) + q'K:2iq')] Ux'{P')l^Ux{P) . (9.5.3) 
Consider now the term in square brackets as ^ 0. From (9.4.4) and (9.4.6, 8) 

hm [2M^i(g2)+g2^2(g')] = -2M ^ A^ + V^/o ^,„;viv 

= [gluon-parton] + [ghost] . (9.5.4) 



Use of the Generalized Goldberger-Treiman relation (9.4.13) then imphes, via (9.5.1), 

that 



lim [2M/Ci(?')+gX(g')] 



<0. (9.5.5) 



In other words there seems to be a cancellation between the 'semi-perturbative' 
gluon-parton contribution and the non-perturbative 'ghost' contribution in (9.5.4). It 
remains a question of some interest as to why this compensation occurs. We consider 
briefly some of the physical arguments presented to explain this. 

That flo — was claimed to be quite natural in the Skyrme model approach to the 
'spin crisis' [BRO 88]. The SU{3) Skyrme model deals with an octet oi pseudoscalar Gold- 
stone bosons whose self-interactions are described by the non-linear Skyrme Lagrangian. 
Quark or gluon degrees of freedom do not occur explicitly and the baryons appear as topo- 
logically non-trivial soliton solutions of the equations of motion. Since there are no quark 
fields it is not absolutely obvious what function of the boson fields corresponds to the 
singlet axial current J^^. Nonetheless arguments can be given for the form of J^^ and in 
[BRO 88] it was claimed that ao is of order 1/Nc and, in fact, ao — 0. (In these papers is 
identified with AS, but that is irrelevant). However, this approach has been criticized by 
Ryzak [RYZ 89] who points out that [BRO 88] utihses an effective Lagrangian calculated 
to order 1/Nc in a 1/A^c expansion, yet which does not permit the rj' to decay although 
such decay is precisely a 1/Nc effect. In addition, in this approximation m^' /M oc {Nc)~^^'^ 
which is far from reality. Inclusion of a term to allow rj' decay would destabilize the soliton 
solutions. Ryzak has shown how to overcome this difficulty. He argues that there are many 
uncertainties in the estimation of oq but concludes that one can have oq = 0.2 ± 0.1 in 
Skyrme type models. However it is hard to see why the 77' couples strongly to the nucleon 
in these models. 

An entirely different approach to the compensation is discussed in [FOR 90], where 
it is argued that the cancellation i.e. the validity of (9.5.5) is natural for instanton config- 
urations of the gluon fields when there is just one fiavour. But it is not at all self-evident 
what happens for Nf — 3. 

Yet another proposal [KUH 90] is that there is a hnk between the so-called confor- 
mal anomaly, which involves Gfj^vG^" and the axial anomaly, which involves Gf^^G^^ . By 
making an assumption about analyticity in the mass of the regulator fermion fields it is 
claimed that one has the following exact result in the chiral limit: 

-2Nf 2 
ao = = — (9.5.6) 

for 3 light flavours. Such a connection between the anomalies and the relation (9.5.6) 
seems quite natural in the framework of supersymmetric QCD where both the energy- 
momentum tensor (with its conformal anomaly) and the flavour singlet current (with its 
axial anomaly) occur in the same supermultiplet of currents. But the extension to ordinary 
QGD is far from obvious, and some counter- arguments have been given [ANS 92a]. So it 
is not at all clear that (9.5.6) can be taken at face value. On the other hand attempts 
have been made [DOR 93] to support (9.5.6) on the basis of instanton dominance of the 
^-vacuum. 

It has also been suggested [DOR 91] that instanton dominance of the QCD vac- 
uum could generate large polarized sea quark distributions. In [DOR 91, 93], it is tacitly 
assumed that there is no anomalous contribution to ao, so the EMG experiment is in- 
terpreted as implying AE ~ and the parameters of the model are adjusted so that the po- 



larized sea cancels the standard Naive Parton Model polarized valence quark contribution. 
The instantons are supposed to generate an effective quark-quark interaction which mim- 
ics Pomeron exchange except that the exchanged object has / = 0, G = C = 1, P = +1 
and flips helicity. By the usual connection between high energy behaviour in off-shell 
Compton scattering and the small-x behaviour of distributions [see Section 4.1.1] it is 
argued that for small x, [Aq{x) + Ag(a;)] oc (xln^x)"^. This unusually rapid growth at 
small X is used to produce the large polarized sea quark contribution. 

However, this approach seems to us highly speculative and is governed by the desire 
to obtain AS ~ 0. We have hopefully convinced the reader that the EMC experiment does 
not demand AS ~ and that, moreover, the Generalized Goldbcrger-Treiman relation is 
nicely consistent with the intuitive belief that the value of AS is not far from 1. 

Finally we comment upon an attempt to resolve the issue by means of a lattice 
calculation [MAN 90] of the anomalous gluon contribution to oq. The idea is to try to 
measure Ag by measuring {P, X\Kf^\P, \) , or something related to it, on a lattice. The 
choice of lattice variables to represent is a non-trivial task but seems to have been 
successfully carried out. The simulation is done on a 6^ x 10 lattice at (3 = = 5.7, 
thus corresponding to as — 0.08. 

In [MAN 90] Eq. (6.4.8) was used to extract Ag, which is found to be small compared 
with what is needed to explain the EMC experiment with a reasonably large AS. Namely 
it is claimed that \Ag\< 0.5. 

However, we know from Section 9.4 that (6.4.8) is not the whole story, and the 
relation is modified by the ghost contribution. Indeed from (9.4.1) and (9.4.3) we see that 
the matrix element of in the forward direction is quite ambiguous - the limit g'^ — > 
must be carefully specified. In fact two different matrix elements were studied, namely 

lim^ (0, S'\K ■ S\P, S) and lim^ S'l^^K^lP, S) (9.5.7) 

where P is a unit vector along P and P — \P\. 

Using (9.4.1) and (9.5.4) the values found for these matrix elements lead to the 
bounds 

3a. A(/ 1 \/3/o (y,,,„vjv 
2n 3 2M 
3asAff V3/o gf^^jvjv 
2tt 2M 

We see that the smallness of the numerical results does not imply that Ag is small. 
Moreover the smallness of (9.5.9) just confirms the result (9.5.5) which was obtained from 
the Generalized Goldberger-Treiman relation and is thus in complete harmony with the 
interpretation of the EMC result where a large AS is cancelled by a significant gluonic 
contribution! Similar comments apply to a more recent numerical study [ALT 94a] . 

However, there is another difficulty with the interpretation of the lattice result [EFR 
92; MAN 92]. The calculation is done in the so-called quenched approximation which means 
that closed qq loops are not taken into account. It is easy to see diagramatically that this 
implies that any rj' propagator in the lattice calculation of the matrix elements (9.5.7) is 
effectively a bare propagator i.e. the rj' mass is not shifted from the value ruQ. There will 
thus be a contribution to (9.5.7) proportional to the bare ry' propagator i.e. proportional 
to {rriQ — (f)"^- Bearing in mind that tuq vanishes in the chiral limit rUq — > we see that 



< — (9.5.8) 
~ 20 ^ ' 

< — ■ (9.5.9) 

^20 ^ ' 



there should be a singularity in the quenched lattice calculation of (9.5.7) in the chiral 
limit as — >^ 0. This suggests that the numerical results should be treated with caution 
until the quenched approximation can be overcome. 

It is interesting to note that an analogous problem would arise in a quenched lattice 
calculation of matrix elements of d^J^^, which is central to the Adler-Bardeen relation. 
In this case the order of taking the hmits 5^ — > 0, — > determines whether or not a 
singularity should occur. Thus, in putting (f — Q and then taking the limit — > one 
finds a cancellation of the singularity between the matrix elements of P [see (9.3.3)] and 
the anomaly. This was not taken into account in the lattice calculations. However, taking 
w?q = and then latting — > produces a singularity in q^G2{q^) of Eq. (9.4.14). 

10 ^'1,2(33) in the QCD Field Theoretic Model 

In Section 3.2 and 3.3 we discussed various attempts to calculate g2{x) in the Naive 
Parton Model and in more sophisticated variants thereof. We emphasized that these cal- 
culations do not lead to coherent results. Different authors obtain different results and 
these are usually incompatible with each other. The basic origin of this difficulty was 
exposed in Section 3.4 where we explained that g2{x) simply cannot be calculated within 
the framework of a Parton Model. 

In this Section we attack the problem from the much more solid basis of the twist- 
3 QCD Field Theoretic Model. We calculate both gi{x) and g2{x) so as to be able to 
emphasize the differences between the two calculations. 

The first 'rigorous' derivation of the result for g2{x) is due to Efremov and Teryaev 
[EFR 84]. It is closely connected with the problem of single spin asymmetries in QCD 
[HEL 91]. This work was somewhat extended by Ratchffe [RAT 86] and by Efremov and 
Teryaev in a later paper [EFR 87; QIU 91]. 

10.1 The QCD Field Theoretic Model 

Let us now turn to the QCD Field Theoretic Model. In this case one deals with 
a Feynman diagram one part of which is handled perturbatively and the other, being 
essentially non-perturbative, is parametrized in terms of a small number of unknown 
functions which play the role of the various parton number densities. 

In the simplest diagram only quarks of 4-momentum k from the hadron are involved 
and the diagram is split into a 'hard' part E^'^ and a 'soft' part as shown in Fig. 10.1. 

Both E and $ are absorptive parts for the reactions 7*5 ^ X and hadron — > qX 
respectively. There are no explicit propagators for the quark lines; they are absorbed into 

The mathematical structure which represents the contribution of Fig. 10.1 to W^'^ 

is 

where 

= I d'z e^^-^{P,S\^p{Q)^P^{z)\P,S) . (10.1.2) 

The label a on V'a(-2) is a spinor label, il^ai^) is a column vector in colour space, so 
involves a sum over colour. Colour is irrelevant in E^" . 



Using matrix notation for the spinor indices, (10.1.1) becomes 



_ 1 r 



Tr 



;io.i.3) 



The Parton Model usually emerges upon making the following approximation: 

1) The 7*g interaction is given its simplest form, as shown in Fig. 10.2. 

In the following we calculate only with the uncrossed Born diagram. The result 
for the crossed diagram is obtained at the end by the replacement xbj ^^bj in the 
hadronic matrix elements connected with and is simply to be added to the uncrossed 
result. For a quark of flavour /, charge e/ (in units of e), for the diagram of Fig. 10.2a, 



[IQ.IA) 



where vfiq is the quark mass. Using Fig. 10.2 in Fig. 10.1 gives rise to the famihar 'handbag 
diagram'. In order to isolate the antisymmetric part of W^'^ we make the replacement in 
(10.1.4): 



and recalling (2.1.8) find 



— m, 
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(10.1.6) 
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2) One assumes that the soft matrix element cuts off rapidly for k'^ off the mass-shell 
k^ — rriq, and for k^^ non-coUinear with respect to the hadron momentum P^. 

This is implemented as follows. Consider a reference frame where the hadron is 
moving at high momentum along OZ so that 



P^" = {E, 0, 0, P) with E 
is a 'large' four-vector. We introduce a 'small' null vector 

\P + E' ' ' 

such that 



\P + E' 
n-P-- 



n 



P + E 
= 0. 



One can then write for k'^ 

kf" 



{k ■ n)P^ + - 



{k-n) 



n' 



~\~ knn 



where 



ki^ = (0, fer, 0) . 

In view of the assumption about $ we can say that 

k^^{k- n)P'' . 



(10.1.8) 

(10.1.9) 
(10.1.10) 

(10.1.11) 
(10.1.12) 
(10.1.13) 



It should be noted that some care is necessary in deciding whether the approxima- 
tion (10.1.13) is adequate. We shall see that this depends crucially upon whether we are 
considering a nucleon with longitudinal (L) or with transverse (T) polarization. 



10.2 gi{x): Longitudinal polarization 

For the study of gi we consider a nucleon with helicity A = ±1/2 and it is sufficient 
to approximate (10.1.7) by putting 

{q + ky ^ + {k-n)PP (10.2.1) 

and dropping the term proportional to the quark mass rriq. Then writing 

qP+{k- n)Pf = J dx 6{x - k ■ n)[q + xP^ (10.2.2) 

we can take the integration over d^k in (10.1.7) through to obtain 

Wj^f = I W / dx il + xPYA^ix) (10.2.3) 
where (using Sl to denote longitudinal spin) 

A'-ix) ^ J -0-d'z6{x-k- n)e''-^P, 5i|^(0)7^75^(^)|P, Sl) 
r d\ 

= —e''^{P,SL\mri5H^n)\P,SL) (10.2.4) 

J ZTT 

is a pseudo- vector which can only depend upon the vectors P^, and = SaPpU-^ 
and the pseudo-vector S'^, and which must be linear in S'^. Given that S ■ P = the only 
possibilities are 5*°" and (n ■ S)P'^. Note that with the normalization 

{P\P') = (27r)^ 2E 6%P - P') (10.2.5) 

A'^[x) has dimensions [M]. 

Recall that for a nucleon with 4-momentum given by (10.1.8) 

2A 

S^{X) = —{P,0,0,E) S' = -l, (10.2.6) 
where A is a helicity (A = ±|) so that 

2A 

^/^(A) = —{P" - M'^n") (10.2.7) 

and we may take 

5^(A)^^P^ (10.2.8) 

i.e. MS''^(A) is a 'large' vector. 

In view of (10.2.8) the structures {n-S)P'^ and >S''^(A) are equivalent in leading order 
and the only possibility is then (the factor 4 is for later convenience) 

A"(x, A) = 4MhL{x) 5"(A) (10.2.9) 

where the dimensionless longitudinal distribution function is given by Q 

Ah^ix) = = / ^ e-hLir) (10.2.10) 



' There appears to be a sign error in the definition of hLix) in [COR 92]. Our hL{x) = -5 h1'='^, the 
factor 1/2 being included so as to obtain (10.2.14). Our sign agrees with [EFR 84]. 



where 

^^(^) = M( ^ Q A P,SL\mM{rnmSL) . (10.2.11) 
M{n ■ Sl> 

Substituting into (10.2.3) and adding the contribution from the crossed Born diagram, 
Fig. 10.2b, yields 

W^^\L) = ej^lkLixBi) + /iL(-XBj)] W?''-5'^(A) . (10.2.12) 

Note that the term xP'' in (10.2.3) does not contribute on account of (10.2.8). 
Consequently (10.2.12) is gauge invariant q^Wl-i^^ = 0. Note that (10.2.8), which holds 
only for longitudinal spin, is crucial for the gauge invariance. 

Comparing with (2.1.19) in the approximation oc we obtain, for the contri- 
bution of a quark of flavour /, 

gi{x) ^^e} [h{{x) + hii-x)] . (10.2.13) 

If one treats the quark fields in (10.2.10) as free fields and regards the nucleon as 
an assemblage of free partons one finds 

hi{x) = Aqf{x) hi{-x) = Aqf{x) (10.2.14) 

so that (10.2.13) reproduces the simple Parton Model result for gi{x). Equation (10.2.13) 
provides a field theoretic generalization of the Parton Model result. A more intuitive 
expression is given in Section 10.5. We now consider how the transverse spin case differs 
from the above. 



10.3 g2ix): Transverse poleirization 

In order to see the essential difference between the longitudinal and transverse spin 
cases consider again the result (10.2.12). In the CM of the 7*-nucleon coUision, as far as 
magnitudes are concerned, for the longitudinal case we have 

\q^\^\MS^{X)\ = 0{uy, P-q = Mv (10.3.1) 

so that for the large components of W'^^\L), 

\w[t\L)\ = 0{v/M) (10.3.2) 

assuming that \hL{x)\ = 0{1). 

In the transverse spin case the analogue of (10.2.4) can only be proportional to 
since n • St — 0, and will produce a result like (10.2.12) with S{X) — > St- Given that 
\St\ — 0{1) one has, for the 'large' components, only 

\K^\T)\ = 0{1) . (10.3.3) 

This immediately suggests that care must be exercised in neglecting non-leading terms 
e.g. in (10.1.11). 

Secondly, note that in (10.2.1) the term {k ■ n)P^ of (10.1.11) did not contribute 
because of the fact that, in leading order, P^ oc S^{\). In the transverse case this will not 
happen and the analogue of (10.2.12) will contain a term e^i^p^P^St in W^^^(T) which, 
(analogously to the Parton Model case (3.4.3) when m ^ niq) is not gauge invariant. 



We must therefore return to (10.1.4) and improve upon the approximation (10.1.13). 
However, it will then turn out that the more complicated non Parton Model diagram 
involving gluon exchange, Fig. 10.3, contributes to the same order. This is not entirely 
surprising given that the operator product result (5.2.2) for the moments of g2{x) involves 
A,. 

Amazingly, as was shown by Efremov and Teryaev [EFR 84], this term just cancels 
the unwanted contribution from the {k-n)]f and the mass terms of the handbag diagram 
and the final result is gauge invariant. Essential in this proof is the use of the equations 
of motion for the quark field. 

The analysis to show the cancellation is rather complicated and is carried out in 
Appendix C. We shall only state the result here. It is the exact analogue of (10.2.12), 
namely, including the contribution of the crossed Born diagram Fig. 10.2b, 

K^^^T) = 4 [MxBj) + M-^Bj)] W?'^? (10.3.4) 
where the analogue of (10.2.9) is 

A"(x, T) = AMfrix) (10.3.5) 

withM 

4/t(x) = / - e^™/T(r) (10.3.6) 

where 

Mr) = j^{P,ST\mi5^THrn)\P,ST) . (10.3.7) 

Comparing (10.3.4) with (2.1.19), for the case of transverse polarization, we obtain for 
the contribution of a quark of flavour /, 

9iix) + g^ix) = ^ [fUx) + fU-^)] . (10.3.8) 

Although the surviving contribution comes from the 'handbag' diagram it does not, 
in fact, have any simple parton interpretation as is explained briefly in Section 10.4. 



10.4 The Naive Parton Model for g2ix) revisited 

The result (10.3.8) for gi{x) + g2{x) together with the expression (10.3.6) for frix) 
which only involves quark field operators might suggest that there ought to be a formula 
for gi{x) + g2{x) derivable in the Naive Parton Model. However one can easily see that 
this is not so. 

In the Naive Parton Model we would treat ip{z) as a free field and the proton as 
an assemblage of free partons. The key element in evaluating (10.3.7) would then be the 
matrix element. 

{k,s'mO)j,^Ti^{Tn)\k,s) (10.4.1) 

taken between free quark states labelled by the covariant spin vectors s and s'. This 
matrix element will be multiplied by proton wave functions giving the amplitudes for 
finding quarks with {k, s) and {k, s') respectively. There are essentially two possibilities: 
i) If we choose to use quark states of definite helicity A, A', for massless quarks 
the matrix element (10.4.1) is diagonal i.e. {k, X'\ . . . \k, X) oc 6x>x and the product of 

Our /t(x) = 2/-^°^(x) for reasons explained after Eqn. (10.2.10). 



wave-functions becomes their modulus squared, yielding a probabilistic interpretation. 
But in this case one finds that (10.4.1) is proportional to rnqSr ■ s(A) which is zero in the 
Naive Parton Model with no k±. So any non-zero result is only of order k± and cannot 
be evaluated without a model of the /cT-distribution, an aspect of the nucleon structure 
which is usually considered to be beyond the Naive Parton Model. 

ii) If we choose to use quark states of definite transverse spin e.g. states 

It or i) = l^(\+)±i\-)), (10.4.2) 



v/2 V 

then the matrix element (10.4.1) is zero for massless quarks if s' = s and is largest for 
non-diagonal transitions | T) — >^ | i)- Consequently one does not obtain a probabilistic 
expression from the nucleon wave-functions. Moreover even the non-diagonal transitions 
vanish in the absence of k±. 

Of course the above arguments do not mean that one cannot do model calculations 
of gi{x) + g2{x). They simply imply that the results will not have the traditional Parton 
Model form and will depend upon specific assumptions about the k± dependence of the 
nucleon wave-function. 

10.5 Probabilistic form for gi{x) in the Field Theoretic Model 

We mentioned in Section 10.2 that the hadronic matrix element hL{x), in terms of 
which gi{x) is given in (10.2.13), reduces to the usual partonic /^q{x) when the fields are 
treated as non-interacting free fields. However even in the case of interacting fields one 
can rewrite the formulae (10.2.10, 11) for hL{x) in such a way as to display manifestly 
the probabilistic interpretation. 

From the analysis of the matrix element {P, SL\'4'{0)j^j5'4'{Tn)\P, Sl) given in Eq. 
(E.l) it is easy to see that to leading order the matrix elements with /i — and /i — 3 are 
equal. From the definition of n'* (10.1.9) we can thus replace the factor ^ which occurs in 
the definition of hL{x) by 

;(/ + f) - i^l'^i- (10-5.1) 



P + E" ' ' P + E ' P 
Inside the matrix element (10.2.11) we may then write 

V'(0)^75V'(Tn) -^V'(0)7°75V'(™) = ^V'^(0)75V'(™) 



^^"2p{Tn) -ip\0)^-Y^ip{Tn)^ (10.5.2) 
^ ^ ^-V(rn)) (10.5.3) 




2 

where, in the last step, we have used 

'1±75\2 /1±75 



Finally, inserting a sum over a complete set of states, using translational invariance 
to shift from ■0(Tn) to ■0(0), and carrying out the integration over r in (10.2.10), we 



obtain, taking a nucleon with helicity A = +1/2, 

hL{x) = q+{x) - q4x) (10.5.4) 

where 



X 



(X|i4^V'(0)|P,A = l/2) 



S[n- Px - (l-x)]. (10.5.5) 



This displays manifestly the probabilistic nature of q±{x) even in the fully interactive 
theory. 

The same approach clearly fails for frix) given by (10.3.6, 7), in terms of which 
g2{x) is given by (10.3.8). 

Concerning our expression (10.5.5) for q±{x), one can derive a more accurate result, 
in which one does not make the approximation that the matrix elements with /j, — and 
/J, — 3 are equal, using the formalism of light-cone fields [MAN 91a]. 

Define 

7^ = ^(7° ±7') (10.5.6) 

and projection operator 

P+ = ^7-7+ Pi = P+ (10.5.7) 

and "good" component of the field 

= P+tp . (10.5.8) 
Then a more accurate expression for q±{x) is 



X 



(X|i4^V'+(0)|P,A = l/2) 



5[n -Px-il- x)] . (10.5.9) 



11 Future experiments 

Throughout this paper we have stressed the role of the axial anomaly in providing 
an unexpected gluon contribution to the first moments F^'"^ measured in polarized deep 
inelastic lepton-hadron scattering. As a consequence of this mechanism the data can 
be explained with a large spin contribution AS ~ 0.7 from the quarks, as expected 
intuitively. This is in accord with the Generalized Goldberger-Treiman relation and with 
various relativistic bag model calculations. The smallness of the flavour singlet part of 
is explained as a cancellation between AE and the anomalous gluon contribution. 

However, in order to accomplish this one requires a large gluon spin Ag at (Q^) ~ 10 
(GeV/c)^, i.e. Ag ~ 4. As we have explained, it is inevitable, within QCD-evolution, that 
Ag grows indefinitely with Q^, so even if intuitively we expect Ag very small at small 
we really cannot rule out the above value at (Q^) ~ 10 (GeV/c)^. 

But, clearly, it is of the utmost importance to test this explanation by measuring 
A^r in other, independent reactions. If it is as large as the DIS experiments suggest then 
there will be measurable effects in several other processes. In the following we indicate 
briefly some of the possibilities that have been studied thus far [WIN 92, REY 93]. 

In Section 7.3 we explained how the growth of A^f is compensated for by an analo- 
gous growth in L^, the total orbital angular momentum carried by all partons. The hard 
processes we are about to discuss, being short range interactions, are insensitive to 
and can thus reveal Ag directly. 



11.1 Semi-inclusive deep inelastic lepton-hadron reactions 

Perhaps the most straightforward proposal [CAR 88, MAN 91] is in semi-inclusive 
polarized DIS.The idea is to select events with final state hadrons of large pt relative to 
the 7*-nucleon collision axis. Such events should be enriched by two-jet production via 
the photon-gluon fusion subprocess shown in Fig. 8.2 and will depend directly on ^g. 
The contribution to T\ coming from this subclass of events should be sensitive to Ag( and 
could even be negative for some choice of pr cut-off. 

Very much in the same spirit one could select on events with charm pair production 
[GLU 88] or on J /if) production at large pt [KAL 89; SRI 92]. The mechanism is the same 
as shown in Fig. 8.2. These processes are quite sensitive to A^f. 

Another proposal is to look at the semi- inclusive production of forward kaons, which 
should be sensitive to the strange quarks in the nucleon [CLO 89]. The idea here is to 
measure the polarization of the strange quarks in the nucleon, i.e. As. This is important 
because if the anomalous gluon contribution is NOT the correct mechanism to explain 
the smallness of the flavour singlet part of T\ then one can arrange to have AE ~ by 
invoking a surprisingly large and negative polarized strange quark distribution, as was 
discussed in Section 4.3.3. However the interpretation of the experiment might be difficult 
because some contribution could arise from ss pairs created in the photon-gluon fusion 
diagram. Fig. 8.2. 

Most of the above measurements should be feasible in the present SMC (CERN) 
and SLAC experiments (provided high enough statistics and good enough particle iden- 
tification can be achieved) and in the next generation of experiments now being planned, 
viz. HERMES at DESY and HELP at CERN. 

11.2 Hadron— hadron reactions 

If it is possible to produce high energy collisions of longitudinally polarized pro- 
tons or antiprotons on a longitudinally polarized nucleon target then a whole range of 
experimental possibilities become feasible for studying A^f. 

11.2.1 Hard 'y and Drell-Yan reactions 

Several possibilities have been analysed (see [BER 89]). The classic and most direct 
experiment would be the polarization asymmetry in the production cross-section for hard 
7's at high pt- The obvious asymmetry 

da~^ — da 



da + da 

depends strongly on Ag and Aq. 

It should also be possible to detect the polarization asymmetry in the production 
cross-section for a high pt jet [RAM 88], but this will be less sensitive to Ag than hard 
7 production because quark initiated processes play a relatively larger role. 

Experiments of the above type are planned for the RHIC collider now under con- 
struction at Brookhaven but it might also be feasible to do such measurements using the 
secondary polarized proton beam colliding with a fixed polarized target at the Fermilab 
TEVATRON. 



(11.2.1) 



One of the most interesting possibilities for learning about /S.g in hadron-hadron 
collisions involves the production of Drell-Yan lepton pairs at large pr, i.e. where the px 

of the pair is large with respect to the collision axis. Of course the obvious experiment, 
analogous to the hard 7 case, is the production asymmetry using both polarized beam 
and polarized target. But what is especially interesting about the Drell-Yan case, given 
the difficulty of providing both a polarized beam and a polarized target at high energies, is 
that certain asymmetries sensitive to /S.g exist even when just the beam or just the target 
is longitudinally polarized. 

The problem is that these single spin asymmetries vanish in the usual Born mech- 
anism for large pr production shown in Fig. 11.1. A non-zero asymmetry arises from 
7 — Z*^ interference for large mass pairs [LEA 93, 94a] or from interference between the 
Born diagrams and the one loop QCD corrections to qg — > and qq — > 751 [CAR 92] . 

The most complete information is obtained if the angular distribution of the lepton 
in the pair rest frame can be measured [LEA 93, 94a]. Alternatively, one can measure 
what might be called the handedness [H] of the lepton pair. Let k± be the laboratory 
momentum of the £^ and P be the laboratory momentum of the proton polarized beam 
(or target). Define 

Xi = P-{k+x (11.2.2) 
and let N{Xi) be the number of pairs with given value of Xi. Then 

^ _ N{Xi > 0) - N{Xi < 0) 
~ N{Xi > 0) + N{Xi < 0) 

is sensitive to Ag. 

All such single spin experiments, however, look difficult because of the small cross- 
sections. It has also been pointed out that there is a correlation between the polarization 
of the produced lepton or direct photon and the polarization of the beam [CON 91], which 
depends on Ag. Such an experiment might be feasible with muon pairs, but clearly it is 
very difficult to monitor the final lepton polarization. 

11.2.2 Heavy quark production 

Probably the most sensitive experiment to A^f and, in many ways, perhaps the 
most realistic, is the polarization asymmetry in the production of heavy quarks using 
longitudinally polarized beam and target [COR 88] . In this case the dominant subprocess 
is the gluon fusion reaction g -\- g ^ q -\- q and the asymmetry will be proportional to 
{AgY. Sizable asymmetries are expected and should be measurable at RHIC. 

11.3 Jet handedness 

If it were possible to measure the polarization of the high pt quark produced in 
the partonic reaction 

q + g^q + g (11.3.1) 

where an unpolarized quark collides with a polarized gluon, then one would have a very 
direct measure of the gluon polarization. The advantage of such an approach is that it 
requires only an unpolarized beam colliding with a polarized target. The problem is to 
measure the final quark's polarization in the dominant process where it hadronizes into a 
jet. Thus one has to study the reaction 



(11.2.3) 



p + p^ jet + X 



(11.3.2) 



and somehow monitor the polarization of the quark which initiated the jet [STR 92]. (One 
is effectively trying to measure the longitudinal spin transfer Dn [BOU 80].) 

One can define the handedness of a jet and show that it is proportional to the 
polarization of the initiating parton [NAG 77; EFR 78, 92a]. Namely, analogously to 
(11.2.2) one can define 

X, ^ • jg^ (11.3.3) 

where P j is the jet momentum and "p^ 2 the momenta of two particles in the jet 
chosen according to some definite criteria. For example we could choose the leading and 
the second leading particles in the jet. 

The handedness is then defined in terms of the number of jets with Xj^ > and 
Xl < 0, I.e. 

NjX, > 0) - NjX, < 0) 
~ N{Xl > 0) + NiXL < 0) ■ ^ ^ ^ ^ 

(An analogous quantity can be defined based on using the momenta of the three most 
leading particles in the jet). 

The essential point is that 

H,,t = ajPg (11.3.5) 

where Pq is the degree of longitudinal polarization of the initiating quark and aj is the 
'jet analysing power', the calculation of which lies in the realm of non-perturbative QCD. 
Thus the key issue is the size of aj. 

Miraculously aj can be measured! The point is that in the electroweak process 

e+e" ^ Z° jets (11.3.6) 

the q and q are significantly polarized and their polarization is known with some confidence 
from electroweak theory. Thus by measuring Hjet in (11.3.5) we effectively measure aj. 
One can argue that it is also possible to measure aj independently of electroweak 
theory by studying the correlations between the handeness of the two jets in (11.3.6). 

There is a theoretical estimate that aj ~ (5 - 10)% [RYS 93, EFR 94]. If this 
is confirmed by measurements of (11.3.5) then it will open up the real possibility of 
measuring Ag at presently existing accelerators, i.e at the Fermilab TEVATRON, the 
CERN SPS and the Serpukhov UNK. 



12 Conclusions 

Polarized Deep Inelastic Scattering has proved to be an exciting and controversial 
field in which unexpected experimental results have engendered the discovery of unex- 
pectedly subtle theoretical problems. 

Much new data has become available in the past years from the SMC at CERN 
and from the E142 at SLAC, including, for the first time, results on the neutron. And the 
controversial EMC proton data of 1987 have been confirmed by the SMC. 
The principal points that have emerged are: 
1. The Bjorken sum rule has been tested, and provided sufficient care is exercised in 
dealing with experiments at different mean Q"^ and in incorporating perturbative 
QCD corrections, there is absolutely no evidence for any violation of the sum rule. 
Earlier claims to the contrary were simply too cavalier in handling the data. 

See [NAC 77]: note that aj corresponds to the quantity in that paper. 



2. The measurement of the first moment of gi{x) does not measure the contribution 
of the quark spins to the spin of the proton. There is an anomalous gluon contribution 
which significantly alters the theoretical expression for so that the smallness of 

no longer implies an almost zero contribution from the quark spins. This expla- 
nation requires a significant polarized gluon density and this feature must be tested 
independently. The polarized gluon number density Ag can be measured in polarized 
semi-inclusive deep inelastic lepton-hadron interactions, and in hadron-hadron colli- 
sions involving polarized beam and target. Both hard 7 and Drell-Yan reactions are 
sensitive to A^f, as is the polarization asymmetry in heavy quark production. Another 
possible approach is via "jet handedness". RHIC seems to be an ideal machine for 
exploring the property of A^^. 

3. While Ti itself no longer determines the contribution of the quark spins to the spin 
of the nucleon, there does exist a non-perturbative Generalized Goldberger-Treiman 
relation which connects the quark spin to the coupling of the 1]' to nucleons. This 
coupling can in principle be deduced from a detailed phase-shift analysis of elastic 
nucleon-nucleon scattering or from some other experiments. Estimates already exist 
but a more accurate determination would be of great interest. 

4. It turns out that g2{x) cannot be calculated in the Naive Parton Model nor in the 
QCD Improved Parton Model. Nonetheless an expression can be obtained for it in a 
field theoretic approach. The result does not have a partonic interpretation. 

5. Neither the Burkhardt-Cottingham nor the Efremov-Leader-Teryaev sum rules follow 
from the Operator Product Expansion. Indeed if simple Regge arguments can be 
trusted, the behaviour at small x of g2 should lead to divergent integrals in these 
sum rules. But the argument is not watertight and every effort should be made to 
test the sum rules experimentally Q. 



Appendix A - Some kinematic relations amongst asymmetries and scaling 
functions 

The spin-spin asymmetries A\\ and A±, Eqs. (2.1.27) and (2.1.38) respectively, are 
usually expressed in terms of virtual Compton scattering asymmetries Ai^2, as 



and 
where 



Au=D{Ai+r]A2) (A.l) 



A^ = diA2-^A^) , (A.2) 
^ _ MuGi - Q'G2 _ gi- {AM^x^/Q^)g2 

D-A^. V-^^ (A.5) 

E{l + eR) ' E-eE' ^ ' 



Note added in proof: The first few data on (?2(^) which have just appeared [ADA 94a] have large 
errors and do not allow yet a significant test of the sum rules. 



with 



and, as in Eq. (2.1.32), 




(A.6) 



The analysis of A\\ proceeds through subsequent approximations. One first notices 
that both the coefficient 77 and A2, which is bounded by [LEA 85] 

IA2I < \fR , (A.8) 
are small, so that it is reasonably accurate to write 

^11 fti DAi . (A.9) 
One then neglects the g2 term in Ai, Eq. (A. 3), and replaces Fi{x) by 

FAx) = .^'^""^ (A. 10) 

^ ' 2x[l + i?(x)] ^ ' 

which originates from Eq. (A. 7) for 4M^a;^ <^ Q^. 
Thus one obtains 

, , A\\ F2(x) , 

«'<^> ijwfm] • 

which expresses gi{x) in terms of the measured asymmetry A\\, the unpolarized structure 
function F2{x)/[2x{l + R)] (also measured) and the known coefficient D. It can be shown 
that all the above approximations are actually quite harmless [LEA 88] . 
The analysis of A± is dicussed in Section 2.1.3. 

Appendix B - Current matrix elements in the Pcirton Model 

In the Standard Model all hadronic currents are expressed in terms of quark fields 
ipf{x) where / labels the flavour. We require the forward matrix elements of these currents 
taken between nucleon states. We shall use the quark-parton model to evaluate these in 
the frame where the proton moves along OZ at high speed. 

Consider the generic current 

Jf = ^fT^f (B.l) 

where F = 7" or 7^75, evaluated at the space-time origin — 0. We shall evaluate the 
forward matrix elements of J/ in a proton state with covariant spin vector S^. 

We insert complete sets of quark states, corresponding to partons of four-momentum 
k and helicity A. 

Because the proton is polarized we need to introduce n/(fe. A; S)(Pk as the number 
of flavour / quarks with three-momentum in the range k — > k + (fk and with helicity 
A inside a proton of four-momentum P and covariant spin vector S^. After some algebra 
one flnds for the forward hadronic matrix elements, 

{P,S\Jf\P,S) = J2 [ d^k — nf{k,X;S){k,X\Jf\k,X) (B.2) 



where e is the parton energy and where we have utihsed the fact that for fast quarks 7"^ 
and 7^7^ conserve hcUcity. 

Let us now define more specifically 

Jf = i^f^^f = V5/7"75V'/ ■ (B.3) 

Then treating the quark-partons as 'free' particles in one has 

{k, \\J'}\k, A) = uxi^'ux = 2k'' (B.4) 

and 

{k,X\Jf%\) = uxl^'l^U), = 2m^s"(A) (B.5) 
where s"(A) is the parton's covariant spin vector corresponding to a helicity state A = ±|, 
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s"(A) = — (IfcUfe). (B.6) 
m/ 

On grounds of Lorentz covariance the forward nucleonic matrix elements have a 
similar structure: 

{P,S\Jf\P,S) = 2vfP'' (B.7) 

and 

(P, 5| |P, S) = 2MafS" (B.8) 

where Vf,af are numbers (Lorentz scalars) that measure the strength of the vector and 
axial- vector matrix elements. As explained in [BAI 82] the numbers Vf are exactly known 
when the currents are conserved currents. 

Substituting (B.4 and 5) into (B.2) and using (B.7 and 8) one finds 



A 

and 



MS"af = Y,jd^k — nf{k, A; S) m/s"(A) . (B.IO) 

Now taking in 5"°° 

(a;'P + ^^^^^, kr, x'P^ (B.ll) 

and using it in (B.6), and taking 

2A / 

^"^^^^mY' °' °' 2pj ^^-^^^ 
for a proton of helicity A = we find eventually, from (B.9 and 10) 



and 



dx' [q^{x') + qj{x')] = dx qf{x) (B.13) 



af = dx' [q^{x') - qf{x')] = £ dx Aqf{x) . (B.14) 

where qf{x) and qf{x) are respectively the number densities of quarks / with spin parallel 
or antiparallel to the nucleon spin. 



In the above we considered only quarks. A similar arguments holds for antiquark 
constituents, where the partonic matrix elements are 

{k, \\J'}\k, A) = -vxi^'vx = -2k'' 

{k,\\jf%\) = -^;a7"75^a = 2m/s"(A) . (B.15) 
Finally then one finds 

Vf^ I dx [qf{x) - qf{x)\ (B.16) 

af= [' dx [Aqfix) + Aqfix)] . (B.17) 



The currents occuring in DIS are all just linear combination of the currents Jf and 
so that all results quoted in Section 4.3 follow trivially from (B.7, 8, 16 and 17). 

Appendix C - Transverse spin: the restoration of electromagnetic gauge 
invariance 

We work in a reference system where the nucleon 4-momentum P^^ is given by 
(10.1.8) and where the null- vector is as in (10.1.9). 

We shall choose the gauge A • n — for simplicity. In the general treatment [EFR 
84] where E'^^{...,k; ...) of (10.1.3) does not refer just to the Born ('handbag') diagram, 
the technique is to replace (10.2.1) by 

{q + k)P = qP + {k-n)PP + MP (C.l) 

and to do a Taylor expansion in the small quantity /Sk^, so that (10.1.1) will involve 
{dEi''')/{dkp) evaluated at kP = {q + xP)p. 

The contribution to W^^ from Fig. 10.3 which involves gluon exchange {G) is of the 
form (we show explicitely the colour labels i, j, a) 



X j d''2i(i''z2e"»"e'('»-''''={P,S|V;«(0)^«(22)'/'a,i(^i)l^',5>. 

Wc absorb the colour matrix (A"/2)jj and the factor g into the hadronic matrix 
element to obtain 



where 



$G ^ y dSiciSs e*^i ^i{k2-k,)z2(^p^ ^i^(O) gAp{z2) i^{zi)\P, S) (C.4) 

and where now is a column vector in spin space and in colour space and Ap is the 
matrix Ap — (A"/2)A^. The trace in (C.3) is only over spinor labels. 

In the general treatment a key role is played by the Ward identity 

——(q + xP)^EPP''(q + xP,q + xP) (C.5) 
okp 



which permits a connection between the leading order part of Wq" and the Taylor ex- 
pansion of W^'' from the quark diagram of Fig. 10.2 . The reader is referred to [EFR 84] 
for the general analysis. Here we shall illustrate the phenomenon by taking the simplest 
case where E'^'^ is given by the Born diagram in Fig. 10.2 and E^^p" by the Born diagram 
shown in Fig. C.l. Both E'^" and E'^'"^ are independent of colour. 

For this simplest case it is easier not to use the Taylor expansion (C.l) in (10.1.4) 
but to split ^ + Jji + mq into ^, which yields the result (10.3.4), and the contribution from 
If, -\-mq, which we shall call /^Wj^^\ where 

AW^(T) -'iJ0^, {e'''''^ni.l5%] - m.IVK'^^] } 5[(fc + g)^ - mj (C.6) 
where we have used (10.1.5) and (10.1.7) and where 

$p = y e**="(P, S\i!{0)idpij{z)\P, S) . (C.7) 

To obtain (C.7) we have used the fact that kp can be obtained by differentiating 
e**^^, plus the usual assumption that the hadronic matrix element vanishes for \z\ — > oo. 
Note for future reference that by translational invariance one can show that 

^(0)z dp ^{z) ^ i>m-i dp)ij{z) (C.8) 

are equivalent inside the matrix element in (C.7). 

Introducing / = J dx d{x — k ■ n), and using (C.8), Eq. (C.6) can eventually be 
written as 



2 

AW^('X) = J dx Six - XBj){e^'''^lM-i dp +i dp) - 2m,(j^0 (C.9) 



-,2 

Q 

where we have introduced the shorthand for any X 

dr 



{X) = l£ e'^^{P,ST\mX^irn)\P,ST) . (C.IO) 

Turning now to the gluon diagram, using the Born approximation in Fig. C.l, we 
keep only the leading terms in the quark propagators: 

qP + k^f=iq'' + {ki ■ n)PP = j dxi S{xi - ki ■ n){qP + XiP") 

q" + k^ q" + {k2 ■ n)P'' = j dx2 S{x2 - k2 ■ n){qP + XaP^ ■ (Cll) 
After some manipulation (C.3) becomes 



2 

= ^Jdx,dx2{{Y{4 + ^2f)g4irinW + ^ifhn) 



xDisc 
where, for any X, 



8P-q 

1 / I 



(Xi - XBi){x2 - XBj) \TrP ■ q^ 



(C.12) 



{{X{x,,X2;n))) = I ^^e^-^^^e^^^(^^-^^) 
X (P, ST\i^iO)X{x,, X2, Tin)V(T2n) |P, St) . 



(CIS) 



Note that if X depends on only one of the variables xi or X2, one has the useful 
results 

' dx2 ((X(xi;ti))) - (X(xi;ti)) (C.14) 



dx, ((X(x2;ri))) = (X(X2;0)). 



(C.15) 



Note that although (C.14, 15) look innocuous there is the hidden assumption that 
the Fourier transforms in (C.IO, C.13) can be inverted. We shall return to this question 
in Appendix E. 

The discontinuity in (C.12) is most simply evaluated by sphtting into partial frac- 
tions. It will turn out to be useful to do this in two ways: 



Disc 



1 



_{xi - XBi){x2 - X-Bi) 



Disc 



1 



[X2- Xi 



or 



X2 — Xi 



Xi — XBj X2 — XBj 

[5{Xi - XBj) - 5{X2 - XBj)] 



(C.16) 



Disc 



(Xl - XB-i){x2 - X-B-i) 



Disc 



1 



1 



1 



Xi + X2 — 2XBj \_Xi — XBj X2 — X-B-y 



-ITT 



The 7-matrix identity 



X1+X2- 2xBj 



[5{Xi - XBj) + S{X2 - XBj)] . 



(C.17) 



(C.18) 



is helpful in simplifying the 7-algebra and will be crucial at a later stage of the analysis. 

When using this to simplify Y{i + X2^)^{4 + Xi^)Y, and when the piece antisym- 
metric in III/ is extracted, one can discard all terms of the type Y{]/'oy ^){P-A or q-A)Y, 
since they give a vanishing contribution in leading order. To see this, via (10.1.5), such 
contributions involve terms like 



(pp or q''){P^ or gj£^'^"'^(P, STm0hal,A,{nn)i,ir2n)\P, St) . 



(C.19) 



Now the matrix element is a pscudotcnsor. The leading terms linear in St can then only 
have the forms PaSj and Si Pp. But the latter must be absent since in our gauge A-n = 0. 
Thus (C.19) involves P ■ St or q ■ St which both vanish. 

Finally, the antisymmetric part of the 7-matrix expression in (C.12) becomes 
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(P-q) {~te^"''''-/^-f5Ap{rn)ixi+X2-2xBi) + {xi + X2)i[P''A^{Tn)-P^A^{rn)]} . (C.20) 



Combining this judiciously with (C.16 and 17) we get for the antisymmetric part of 



(C.12) 



lP~q I {g5{xi-XBi){{-e^'""'i.-i^Ap{Tn) 

- %i\P^A^{Tn) - P^A\Tn)\)) -gd{x2 - a;Bj)((e'^''^'^7a75^p(™) 

- i^[P"^'^(rn) - P'^A'iTn)]))} . 



(C.21) 



The trivial dependence on xi,X2 allows us to use (C.14 and 15) to obtain 



= ^ / '^(^ - ^Bj){(-£'^'^'"^7.75^p(0) 

- iii[P''A''{0) - P'^A^(O)]) - {e>''""'^,^5Ap{Tn) 

- iT/tiP^A^irn) - P^A^irn)])} . 



(C.22) 



Now note that the structures of AW^^X) in (C.9) and W^"^^^ in (C.22) are identical, 
so we may add them to obtain 



- m^a^^ + igi\P^A\{S) - P'^A^lO)])} (c.23) 
where we have introduced 

Bp(^) = i^p-gAp{z) 

hp{z) = -idp-9Ap{z). (C.24) 

Consider first the terms involving ^. We can simplify them by analysing the structure 
of the matrix element 

{P^Srli^rA-^ijlP^ST). 
It is a true tensor which can only be built from the vectors P", and 

p.v = n-v^S^ -v^O , (C.25) 

and the pseudovector S"^, and which must be linear in S^. The only possible forms are 
P^^v'^ or P'^v^ of which the latter cannot occur because of the gauge condition A - n — 0. 
It follows that 

if A'') = {^tA") = {t^A") = . (C.26) 
If therefore we write, to leading order, 

^A" = [iP^ + fn'^ - $tS^ - iv^]A'' (C.27) 

then inside (C.23) only the term i/iP^ contributes. Inside the matrix element we thus have 
the equivalence 

^PM^ <^ ^A'' . (C.28) 
We may thus, in (C.23), make the replacement 

iig[P>'A'' - P^A^"] YgA" - YqA^ . (C.29) 

We may also, as a consequence of (C.8) add {dp + dp) so that e.g. 

{i^i'gA''{Q) - i^i^gA^irn)) = {-i^i" d\q) + il^ D\rn)) . (C.30) 

Finally, (C.23) becomes: 

AW^% + Wl^l^^ = ^ / dx 6{x - XBj){(5'^^^'^7.75 Bp(rn) - m,a^-^ 

+ ij^ D\Tn) - iY D^rn)) + (e''"^^7a75 Bp(0) 

- rUga'"' -iYD''{0)+iY d\0))}. (C.31) 



We shall now demonstrate that both matrix elements in (C.31) vanish as a conse- 
quence of the equations of motion 



[IP {z) - m,]ij{z) = Q i,{z)[ip {z) - m,] = Q . (C.32) 
The starting point is the 7- matrix identity which follows from (C.18): 

= -a^'^'Y + ig^f-i^" - ig^P-i" + £^^^^"7,75 . (C.33) 
Taking the scalar product with Dp{z), acting on ipi^)^ using (C.32) yields 

= [-m,(T^'^ + £^'^''"7<x75 Dpiz) + lY d\z) - lY d\z)]^{z) (C.34) 

which is exactly the structure that appears in the first matrix element of (C.31). A similar 
argument shows that the second matrix element also vanishes. 

Thus the miracle is achieved. The unwanted and mass terms have disappeared 
as a consequence of the equations of motion when the two parton handbag diagram is 
combined with the three parton gluon diagram. The result is (10.3.4). 

Appendix D - Distribution and two-parton correlation functions for transverse 
spin 

It is of some interest to consider the distribution and correlation functions associated 
with the matrix elements of the quark and gluon operators appearing in Appendix C and 
to derive certain relationships between them which follow from the equations of motion. 
The structure of these correlation functions will be further studied in Appendix E. 

Consider first (a'^'^) in the notation of (C.IO). Its leading order tensorial structure 
can only be e'^^°''^ PaSj . Thus we put 



and find that Q 



where 



{a^'^) = 4/iT(x)e'^""^P«sJ (D.l) 



4/it(x) = / — e^™ hrir) (D.2) 



hrir) = {P,STmO)MTHrn)\P,ST) . (D.3) 



It is hxix) that is the true analogue of h^lx) of (10.2.10) and which has a simple 
parton interpretation in terms of probabilities to find quarks with transverse spin in 
a transversely polarized nucleon. As we have seen, however, in DIS this term appears 
multiplied by the quark mass [Eq. (10.1.7)] and in fact eventually cancels out. Nonetheless 
it is a fundamental distribution reflecting a different aspect of the structure of the nucleon, 
and, as has been emphasized by [COR 92] , can be measured directly in Drell-Yan reactions. 

Consider next the matrix element ((7^^ D (rn))). In leading order one can only have 
the tensor structures P'^f ^ or P'^v'^ but one can show that the latter term is not permitted 
by time-reversal invariance. So one can put 

((7'' D\rn))) = MB^ixi, X2)P^v'' . (D.4) 



^) Our hrix) = l/2h^"^{x) for reasons explained after Eq. (10.2.11). 



It follows that 
where 

B''in,T2) = -j^^{P,ST\m4iv-Dinn)i,{T2n)\P,ST). (D.5) 
Note that it can be shown [EFR 84] that as a consequence of time reversal invariance 

B''{x,,X2)^-B''{x2,Xi) . (D.6) 

Now using (C.14) one has 

{iY D\rn) - D^irn)) ^ iM J dx2 B^{xi, X2) e^"'''^PaS^ (D.7) 

where we have used 

= e'^^'^fp^S^. (D.8) 
Finally, consider {{j'^j^ D {rn))). In leading order, one can only have 

((7^75 BVi^)» = MB^{xi, X2)P''S't + MB^{xi,X2)P'S^ (D.9) 



with 



B^{xi, X2)= fP'^ e^"^"=^ e'^'^"'^B^{Ti, T2) (D.IO) 

J ZTT ZTT 

B\t^.T2) = -^^P,ST\milAS^- D{T^n)\i>{T2n)\P,ST) (D.ll) 

B^{x^,X2) = J e'^'^' e^^^(^^-^^)S^(Ti,T2) (D.12) 

B^{rur2) = -^(P,5T|VS(0)M5ri[n • L'MlV'MlP, 5^) ■ (D.13) 
From time reversal invariance one has [EFR 84] 

B^(xi,X2)^B^(x2,xi). (D.14) 
Note that in (D.12), because of the gauge condition, 

n- B(Tin) = in^d^ , (D.15) 
so that the ti integration can be done yielding 

B^{x^, X2) = -^5{x, - ^2) / ^ e-^-^ 

x{P,ST\'ipiO)^Tl5n''daipiT2n)\P,ST) ■ (D-16) 

It is fairly obvious, and is shown in Appendix E, that the derivative n ■ d acting 
on ip is equivalent to d/dT2 which may then be switched so as to act on the exponential, 
yielding 

B^(xi, X2) = 45(xi - X2)x2Mxi) (D.17) 



where 



and 



where 



where /r is the function appearing in (10.3.6). 
Using (C.14) and (D.9) we have 

(£'^'^'"^7.75 Dpirm)) = M jdx2 {B^{xi, X2) - B^{x,, X2)} s^'^^'^PaSl • (D.18) 
The vanishing of the matrix element of (C.34) thus impUes, using (D.l, 7 and 17) 

dx2 {B^ixi, X2) - B^{xi,X2) + iB^'ixi, X2)} = j^hrixi) . (D.19) 
Taking account of (D.l 7) this becomes: 
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dx2 {^5^(xl, X2) - B\x^, X2)] = j^hrixi) - 4xJt{xi) • (D.20) 

Use of the equations of motion acting to the left upon ^(0) eventually leads to an analogous 
relationship: 

dxi {iB^ixi, X2) + B^{xi, X2)} = J^hT{x2) + Ax2fT{x2) ■ (D.21) 

Appendix E - The Burkhcirdt— Cottingham and the Efremov— Leader— Teryaev 
Sum Rules in the QCD Field Theoretic Model 

We shall show how the above sum rules can be derived from a careful study of the 
structure of the matrix elements involved in Section 10 and in Appendix D. The other 
ingredient is the equation of motion for ip{z). Also one has to make assumptions about 
the invertability of certain Fourier transforms. One sees very clearly where the weak point 
is in the 'proof of these sum rules (see Sections 5.4 and 5.5). In the case of the Efremov- 
Leader-Teryaev sum rule, the original derivation [EFR 84] was a little unsatisfactory, 
because it appealed to a particular reaction to derive properties inherent to the nucleon. 
The present derivation deals only with nucleon matrix elements. 

Consider first the forward matrix element of the bilocal operator 

on the hght cone z"^ — 0. Its most general form is 

^(V^(0)7'^75V^(^))p,5 = A^S>' + {z ■ S)A2P^ + {z ■ S)A^z'' (E.l) 

where ( ... )p^s is short for (P, 5*1 ... \P-iS). The scalar functions Ai^2,z are functions only 
of z • P. ' 

Prom (E.l) we deduce 

-^{miV^{z))p,s = A'^S^P^ + A2P^S'' + A^z^S'' 

+ {z- S)[A^P^'P'' ^A'^z^'P'' ^A^g^'''] (E.2) 



where 



(E.3) 



We now put — rn^ as required in the matrix elements involved in Section 10 and 
Appendix D. Then 

^{mi'l^i^{rn))p,s = A,S^ + r(n • S)[A^P^ + rA^n^] (E.4) 
where now Ai — Ai(r), and 

^(V;(0)7'^759''V'M)p,5 = A'^S^^P^' + A^P^S^ + riA^n^'S'' 

+ {n- S)[A'^Pi'P'' ^TA'^ni'P'' + A^gi''']} . (E.5) 

We assume that all scalar functions arc such that tA{t) ^ as r — >■ for all terms 
occurring in (E.4 and 5). Then at r = we have the simple structures 

^(^(0)7''75^(0))P,5 = ^1(0)5'' (E.6) 

and 

^{miVm)p,s = ^'i(0)^'^P^ + AMP'S"" . (E.7) 
We shall also require, from (E.5) 

^(Vi(0)75#(Tn))p,5 = -T(n • S)[M^A;^ + 5^ + tA^] (E.8) 

so that at r = 

^(^(0)75#(0))p,5 = 0. (E.9) 
Finally note, from (E.5) that 

^(^(0)7^75n-a^(Tn))p,5 = j^,S^ + {n ■ S)[{A2 + tA'^)P^ + t{2A, + TA'^)n^ 

= ^^(V5(0)7''75V'M)p,5 (E.IO) 

a result used in Appendix D. 

Consider now the gluonic matrix element 

1/M {i:{Q)Yi,gA^{z)i^{z))p,s 

with z = TTi. Its most general form is 

T{S-n) [BiP^P" + rB^P^n" + tBsU^P" + r'^B^n'^rf] 

+B5S^'P'' + B'^P'^S'' + TB'iS^'rf + tBsu'^S" . (E.ll) 

The gauge condition rinA^ = implies that 

^5 = 0, rBi^-Be, tB^^-B^ (E.12) 

so that 

^{mi'l^9A''{Tn)ij{Tn))p,s = TB^[{S-n)P^P^-P^S^] 

+t{S ■ n)[B2P^''n!' + rB^n^'n''] (E.13) 
WB2,[{S ■n)n^P'' - n>'S'']+TB^S>'rf . 



Notice the crucial feature, that the imposition of the gauge condition, together with 
the assumptions about the vanishing of products hke tB{t) as t — > 0, leads to the 
vanishing of (E.13) at r = ie. 

(VS(0)7^75^^'^(0)V(0))p,5 = . (E.14) 

This result will be crucial for deriving the Efremov-Leader-Teryaev sum rule. 

Let us now relate some of the above coefficients to the functions occurring in the 
discussion of gi and g2. 

Prom (10.2.11) and (E.4) we have 

/iL(r) = Ai(r)H-ryl2(r). (E.15) 

Prom (10.3.7) and (E.4) 

/T(r)=A(r). (E.16) 
Then from (10.2.13) and (10.2.10), if the Pourier transforms can be inverted, 

-,2 



Jo o 



p2 



= -^Ai(O) by (E.15). 



(E.17) 



Similarly, from (10.3.8 and 6) 



e2 

dx [gi{x) + g2{x)] = -^/t(O) 



= -jA,(0) by (E.16). 
Equations (E.17 and 18) imply the Burkhardt-Cottingham sum rule 



(E.18) 



Cdxg2ix)^Q. (E.19) 
Jo 

We have already discussed in Section 5.4 why the above derivation may fail because 
of the non-invertability of the Pourier transforms. We turn now to the Efremov-Leader- 
Teryaev sum rule. 

Consider first Eq. (D.20) which followed from the equations of motion. Integrating 
over xi, using (D.6), (D.IO) and (10.3.6) there results: 



B^(0,0) '"^^^ 



dr 



(E.20) 

r=0 



where we have taken the quark mass to be zero for simplicity and where we have taken, 
on the basis of (10.3.6), 

4x/T(x) = z/^e-^^(r). (E.21) 
Now because of (E.14), from (D.ll) 

^^(0,0) = -^(v;(o)^75(^T • d)m)p,s 

so that via (E.7) 

5^(0,0)^^2(0). (E.22) 



Use of this and (E.16) in (E.20) yields 

^ = ^'1(0) . (E.23) 



r=0 



Now by arguments similar to those that lead to (E.21), we have 



8 /-i .dhL^T 
ax xgAx) — i 

e} Jo 



T=0 



dr 

= t[A[{0) + A2m by(E.15) (E.24) 
= 2a; (0) by (E.23). 



Similarly we have 



T=0 



(E.25) 



^2/ dx x[gi{x) + g2{x)] ^ 2i—fT{T) 
ej Jo dr 

= 2a; (0) by (E.16). 
Subtracting (E.24) from (E.25) yields the Efremov-Leader-Teryaev sum rule 

r dx x[gi(x) + 2g2(x)] = . (E.26) 
^0 

The same caveats apply to this 'proof as do to the Burkhardt-Cottingham sum rule. 
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Figure captions 



Fig. 2.1 Feynman diagram for inelastic Icpton-nucleon scattering: IN — >■ IX. 
Fig. 2.2 Definition of the angles a, (3, 9, (p, and ©. 
Fig. 2.3 Further definition of the angles a, 9 and (f). 

Fig. 2.4 Spin configuration relevant to Eq. (2.1.39) for the measurement of A± with 

= 0. 

Fig. 2.5 a) Lowest order Feynman diagrams for weak interaction contributions to IN — > 
IX and b) to C^N u(P)X. 

Fig. 3.1 Predictions for xg2{x) and x[gi{x) + g2{x)] based on Eq. (3.3.11) [JAC 89]. 

Fig. 3.2 Prediction for g2{x) in the bag model [JAF 91]. Dotted hne = twist-2 contri- 
bution; dashed = twist-3; solid = sum. 

Fig. 3.3 Interaction of hard photon with quark. 

Fig. 4.1 EMC (•) [ASH 89] and SLAC-Yale (o) [ALG 78; BAU 83] data on Ai{x). 

Fig. 4.2 EMC data on xgi{x). Also shown is the extrapolation to a; = of the integral 

of gi{x) (see text). 

Fig. 4.3 Comparison of the new SMC data on g^{x) at {Q^) = 10 {GeV/cf [ADA 94] 
with EMC data. 

Fig. 4.4 Behaviour of Ip{Q'^) (curve 1) and /p(Q^) — IniQ"^) (curve 2) at large and 
as given by the GDH sum rule at = (cross and triangle respectively, from [ANS 
89]). 

Fig. 4.5 The resonance contribution to I{Q'^) for proton, neutron and their difference 
(from [BUR 93]). 

Fig. 4.6 SMC data on the deuteron asymmetry Af at = 4.6 (GeV/c)^ (from [ADE 
93]). 

Fig. 4.7 SMC data on xgf{x) at = 4.6 (GeV/c)^ (from [ADE 93]). 

Fig. 4.8 E142 data on the neutron asymmetry ^4^(0;) and on gi{x) at = 2 (GeV/c)^. 
See text (from [ANT 93]). 

Fig. 4.9 g{{x,Q^) at various as extracted from the deuterium and helium-3 data 
by [ELL 93] (see text). Continuous curve, = 2 (GeV/c)^; dotted, = 4.6; dot- 
dashed, ^ 10.7. 

Fig. 4.10 gi{x,Q'^) as extracted by [ELL 93] (see text). Continuous curve, Q'^ — 2 
(GeV/c)2; dotted, ^ 4.5, dot-dashed, ^ 10.7. 

Fig. 6.1 Triangle diagram giving rise to the axial anomaly. 

Fig. 6.2 Diagram giving rise to renormalization of axial anomaly. 

Fig. 6.3 Diagram giving rise to renormalization of the axial current. 

Fig. 8.1 Born diagram for 7*g — > q. 

Fig. 8.2 Lowest order diagram for ^*g — > qq. 



Fig. 8.3 Dependence of dAa'^*'^/dk^ on kr for several choices of the quark mass. 
Fig. 8.4 Lowest order diagram for 7*g g q. 

Fig. 8.5 Result of fitting gi{x) with As = 0, and Ag constrained by the unpolarized 
glue. The sohd line corresponds to F/D = 0.63 and the dashed hue to F/D = 0.55. 

Fig. 8.6 Result of fitting gi{x) with As ^ 0, and /S.g constrained by the unpolarized 
glue. The sohd line corresponds to F/D — 0.63 and the dashed hue to F/D — 0.55. 

Fig. 8.7 Zq exchange contribution to elastic up up. 

Fig. 8.8 Anomaly contribution to elastic up up. 

Fig. 9.1 The QCD potential V{X, Y, ...) as a function of X. 

Fig. 9.2 Graphical representation of the non-perturbative coupling of the bare axial 
ghost to the current X^. 

Fig. 9.3 Graphical representation of the coupling between the axial ghost and the field 

rjo{x). 

Fig. 9.4 Graphical representation of matrix elements a) {T[G^Gi,]) and b) (TfG^Ty']). 

Fig. 9.5 Graphical representation of matrix element (P', A'|ir^(0)|P, A). 

Fig. 9.6 Graphical representation of interaction of rj^ and ghost G° with nucleon. 

Fig. 10.1 Separation of 7*+ hadron — > X into soft and hard parts. 

Fig. 10.2 a) Born diagram for 7*5 — > 7*5 and b) the crossed version. 

Fig. 10.3 DIS interaction involving quark-gluon correlation. 

Fig. 11.1 Conventional Born mechanism for large pt Drell-Yan pair production. 

Fig. C.l Born diagram for E^^f"^. There is also a crossed diagram analogous to Fig. 
10.2b. 
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